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ABSTRACT 


Parameter  plane  techniques  were  first  introduced  in  an  IEEE  paper 
dated  November  4,  1964.  The  paper  dealt  mainly  with  the  theory  of  the 
parameter  plane  whereby  the  roots  of  a  polynomial  could  be  determined 
graphically  in  terms  of  two  parameters  which  may  appear  linearly  in 
any  of  the  coefficients. 

In  this  text,  parameter  plane  techniques  are  applied  to  the 
compensation  of  linear  feedback  control  systems  by  both  graphical  and 
analytical  means.  Parameter  plane  equations  are  extended  to  include 
parameter  products  and  three  parameters.  An  attempt  is  made  to  show 
the  complementary  roles  of  the  parameter  plane  and  root  locus. 

The  writer  wishes  to  express  his  appreciation  for  the  assistance 
and  guidance  given  him  by  Dr.  G.  J.  Thaler  of  the  U.  S.  Naval  Post¬ 
graduate  School  in  this  investigation. 
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Introduction . 


1  . 

The  analysis  and  synthesis  of  feedback  control  systems,  or  the 
compensation  of  same,  can  be  effected  by  three  general  methods.  The 
first  of  these  can  be  called  the  integral  criteria.  Here  a  cost  function, 
in  which  is  inherent  the  system  design  specifications,  is  minimized  with 
respect  to  certain  variable  system  parameters.  This  method  is  mainly 
applied  to  the  statistical  properties  of  feedback  control  systems.  The 
second  method  is  the  Bode  frequency  response  method  whereby  the  systems 
open  loop  transfer  function  is  manipulated  to  obtain  the  desired  system 
response.  This  method  has  its  inherent  weaknesses,  such  as  difficulty  of 
application  to  non-unity  feedback  control  systems,  difficulty  in  inter¬ 
preting  closed  loop  transient  response  in  terms  of  open  loop  frequency 
response,  difficulty  of  varying  more  than  one  parameter,  and  the  short¬ 
comings  inherent  in  the  approximations  which  are  required  when  applying 
the  frequency  response  method  to  compensation  problems.  Third  are  the 
algebraic  methods.  Under  this  heading  can  be  listed  the  root  locus  method. 
The  shortcomings  of  this  widely  known  and  widely  used  method  are  familiar 
to  its  users.  Its  greatest  disadvantages  lie  in  plotting  the  actual  locus 
of  roots\  and  in  the  fact  that  only  one  parameter  can  be  varied  conven¬ 
iently.  In  references  (1)  and  (2),  Ross-Warren  and  Poliak  respectively 
developed  algebraic  methods  of  cascade  compensation  using  root  locus 
techniques,  but  the  inherent  disadvantages  of  the  root  locus  were  still 
present.  In  reference  (3)  Mitrovic  developed  an  algebraic  and  graphical 


computer  program  is  introduced  in  section  (5)  which  makes  this 
disadvantage  less  significant. 


1 


method  of  obtaining  the  roots  of  a  polynomial  in  terms  of  two  variable 
parameters.  Later  in  reference  (4)  Ohta  developed  some  additional 
sketching  techniques  which  greatly  facilitated  the  plotting  of  the 
Mitrovic  curves.  In  references  (5)  and  (6)  Choe  and  Hyon  respectively 
applied  and  extended  the  Mitrovic  method  to  the  compensation  of  linear 
continuous  feedback  control  systems.  The  inherent  disadvantage  of  the 
Mitrovic  method  is  that  the  variable  parameters  may  appear  in  no  more  than 
two  coefficients  of  the  characteristic  equation,  which  reduces  the  flex¬ 
ibility  of  the  method.  In  reference  (7)  Siljak  introduced  a  method  of 
obtaining  the  roots  of  a  polynomial  in  terms  of  two  variable  parameters 
which  can  appear  in  any  and  all  the  coefficients  of  the  polynomial.  In 
this  text,  Siljak 's  method  is  applied  and  extended  to  the  compensation 
of  linear  continuous  feedback  control  systems.  General  methods  of  com¬ 
pensation  will  be  developed  and  an  attempt  will  be  made  to  relate  the 
root  locus,  and  the  parameter  plane  as  a  set  of  complementary  techniques 
which  when  applied  in  conjunction  with  one  another  represent  the  most 
adequate  tool  to  date  for  solving  the  problem  of  compensation  of  linear 
feedback  control  systems. 

The  relationship  between  being  able  to  place  the  roots  of  a  polynom¬ 
ial  at  specified  locations  in  the  S-plane  and  the  compensation  of  feed¬ 
back  control  systems  is  as  follows.  The  basic  idea  is  that  any  feedback 
control  system,  including  any  added  compensators  which  may  contain  vari¬ 
ables,  can  be  reduced  to  or  can  be  represented  by  a  ratio  of  two  poly¬ 
nomials  which  is  the  closed  loop  transfer  function.  Well  known  methods 
are  available  whereby  a  specified  system  response,  in  terms  of  over¬ 
shoot,  bandwidth,  settling  time,  steady  state  accuracy,  etc.,  can  be 
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obtained  by  placing  a  pair  of  complex  conjugate  roots  of  the  character¬ 
istic  equation  at  a  specified  location  in  the  S-plane,  while  ensuring  that 

the  real  part  of  this  complex  root  pair  (called  the  dominant  roots)  is 
smaller  in  magnitude  than  the  real  parts  of  the  remaining  roots  of  the 
characteristic  equation.  The  problem  of  compensation  therefore  reduces 
itself  to  one  of  moving  the  roots  of  the  characteristic  equation  to  the 
desirable  locations.  The  usefulness  of  the  parameter  plane  to  achieve 
this  will  soon  become  apparent. 
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2.  Derivation  of  dio  basic  parameter  plane  equations. 

A  feedback  control  system's  characteristic  equation  can  be  re¬ 
presented  as  a  polynomial  of  the  following  form: 

f(S)  =  £  V  *  0  (2-1) 

K-O 

Where  the  coefficients  a,(k  =  0,  l,....,m)  are  real,  and  S  is  the 

I - 21 

complex  variable  S  =  -G+  jw  =  -  J  w  +  jw  ||  1-  ^  .  (2-2) 

w  is  the  undamped  natural  frequency,  and  J w  is  the  relative  damping 
coefficient.  It  is  noted  in  reference  (7)  that  S  may  be  represented 
by  the  following: 

Sk  =  wk(Tk(-  7  )  ±  j  /  1-  72  Uk(-  7  ))  (2-3) 

where  TR(-  J  )  «  (-l)kTk(  J  )  and  UR(-  J  )  =  (-l)k+1Uk(  J).  J  ) 

and  Uk(  J  )  are  Chebishev  functions  of  the  first  and  second  kind  respect¬ 
ively.  Values  of  zeta  will  be  considered  such  that  0  —  -  1  and 

values  of  w  such  that  0  s  w  s  .  The  values  of  Tk  and  Uk  are  tabu¬ 
lated  in  Appendix  I  for  various  values  of  zeta.  But  more  useful  to  digi¬ 
tal  computer  employment,  they  can  be  obtained  from  the  following  recursion 
relations: 

Tk  +  i<  7  >  -2  V  7  >  -  W  }  >  - 0  <2-5> 

uk  + 1<  t  >  -2  1  >  -  \.i<  y  >  - 0 

Here  To<  J  )  -  X ,  T,  (  J  )  -  J  ,  U0(  J  )  -  0,  Uj(  f)  -  1 

Substituting  equation  (2-3)  into  (2-1)  and  setting  the  real  and  imagin¬ 

ary  parts  to  zero  independently  one  obtains: 
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(2-6) 


&0W  V  •  7  >  '  0 

HvV-  ?  >  ■ 0 


Employing  equations  (2-5)  one  obtains  from  equations  (2-6): 


I.(-d  kvVi<  ?  >  ■ 0 

|o(-l)kV\(  ?  )  -  0  (2-7) 

Consider  the  coefficients  of  the  characteristic  equation  (2-1)  as 
linear  functions  of  the  variable  system  parameters  as  follows: 

\  =  bk^  +  ckp  +  dk  (28) 

Employing  the  above  relation  for  a^,  equations  (2-7)  give  the  following 
relations : 


Where: 


B 


1 


oCB][  +  jB  Cx  +  Dj  ■  0 

oL  B2  +  (2  C2  +  D2  =  0 


w  k  k 

£  (-l)VX  . 

K  =  o  k  k-1 


^  k  k 
B  =  £  (-l)kbkw  U. 
2  k?o  K  k 


(2-9) 


ct  *  Io(-1)kckuK-i  v  I^-dVX  <2-io> 
D1  *  j^^VVl  D2  -  Jo  <-1>  V\ 

Since  equations  (2-9)  are  two  linear  equations  in  the  two  unknowns  alpha 
and  beta,  Cramer's  rule  can  be  applied  to  obtain: 


oL  = 


ClVC2Dl 

B1C2'B2C1 


B2DrBiD2 

B1C2'B2C1 


(2-11) 


Equations  (2-11)  are  now  functions  of  zeta  and  w.  Hence  by  fixing  w  and 
varying  zeta  or  by  fixing  zeta  and  varying  w,  the  constant  w  or  constant 
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zeta  5  piano  contours  respectively  can  be  mapped  into  the  real  domain 
of  the  alpha  beta  plane  or  parameter  plane. 

In  reference  (7)  the  following  relationships  are  utilized: 


Sk-P.+Jw/l-  J2  Q. 


P.  . .  +  2w  P,  +  w  P  =0 
k  +1  k  k-1 


(2-12) 


Qk  +  1+  2w  Qk  +  w  Qk  1  =  0 


Pq  =  1,  Pi  =  -w J  ,  0Q  =  0,  Qx  =  1 

pk  -  -w  i  vVi 

Here  P^  and  are  related  to  the  Chebishev  functions  by 

Pk  =  7  }  =  (’1>kwkTk(  7  )  (2-13) 

_  k- 1__  f  —p  v  f  ~  1  k-l„  /  ~7  \ 

\  -  ”  v-  ?  >  ■  <-p>  ”  v  i  > 

By  using  equations  (2-12),  and  (2-13),  one  obtains  proceeding  as  before: 


m  w 

£oWi  =  0  K#oVk  =  0  <2‘14> 

Employing  equations  (2-8),  (2-14),  along  with  Cramer's  rule  one  again 
obtains  equations  (2-11)  where  the  following  expressions  now  apply: 


Jo  VW 

B2  = 

Jo  Vk 

w\ 

^  °k^k-l 

K :0  K  K  1 

C2  = 

^ L  ck^k 

k*o  K  K 

jj? 

£  Wi 

k-o  K  K  1 

°2  = 

£  dkQk 
k-o  k  K 

(2-15) 
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Equations  (2-11)  and  (2-15)  arc  useful  for  mapping  constant  zeta-omega 
curves  from  the  3~planc  into  the  parameter  plane.  As  will  be  seen  later 
these  curves  play  an  important  role  in  dominance  considerations. 

If  the  complex  variable  S  is  substituted  in  equation  (2-1)  by  lett¬ 
ing 

S  =  -  (2-16) 


where  sigma  corresponds  to  values  of  S  along  the  real  axis,  then  in 
accordance  with  equations  (2-8),  the  characteristic  equation  (2-1)  be¬ 
comes  : 


w  k  k 
oL  £  + 


e  k  - 


£  (-l)kd,^k  =  o  (2-17) 

K=0  k 


K-0  K  1  K-0 

The  above  represents  a  straight  line  in  the  alpha-beta  plane  for  a  given 
value  of  sigma.  Hence  a  point  on  the  real  axis  in  the  S-plane  maps  into 
a  straight  line  in  the  alpha-beta  plane.  Also  for  a  given  value  of  alpha, 
beta,  and  sigma  which  satisfies  equation  (2-17),  then  the  characteristic 
equation  (2-1)  must  have  a  real  root  at  minus  sigma.  For  the  constant 
zeta  and  omega  curves  as  defined  previously,  if  for  certain  values  of 
alpha  and  beta,  say  for  a  value  obtained  from  equations  (2-11)  with  a 
certain  value  of  zeta  and  omega,  then  the  characteristic  equation  (2-1) 

I - T 

has  a  pair  of  complex  roots  at  S  =  -  J  w  +  jw ^  1-  J  . 

It  is  important  to  note  that  by  applying  equations  (2-11)  and  (2-17) 
one  can,  for  a  specified  value  of  zeta,  omega,  and  sigma,  compute  the 
value  of  alpha  and  beta  so  that  the  characteristic  equation  will  have  a 

I - P 

pair  of  complex  roots  at  say  S  =  -  J ^w^  +  jw^  J  1-  ,  and  a  real 

root  at  S  =  •  The  m-3  remaining  roots  of  the  characteristic  equa¬ 

tion  can  then  be  determined  by  dividing  out  the  three  known  or  specified 
roots.  This  method  where  zeta,  omega,  and  sigma,  or  just  zeta  and  omega 
are  specified,  and  the  computations  for  alpha  and  beta  are  done  algebrai¬ 
cally,  will  be  referred  to  as  the  albegraic  parameter  plane  solution. 
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To  solve  the  problem  in  general,  for  all  values  of  zeta,  omega,  and  sigma, 
it  is  necessary  to  plot  a  family  of  parameter  plane  curves  for  various 
values  of  zeta,  omega,  sigma,  and  if  desired,  zeta-omega.  On  the  result¬ 
ing  parameter  plane  plot  one  can,  by  picking  an  M  point  or  operating 
point,  graphically  read  from  the  curves  the  values  of  the  m  roots  corres- 
ponding  to  the  m  order  characteristic  equation.  This  latter  method  will 
be  called  the  graphical  parameter  plane  solution. 

The  algebraic  solution  has  the  advantage  that  the  labor  of  plotting 

2 

the  curves  can  be  avoided,  but  it  has  the  disadvantage  that  without  the 
curves  it  is  sometimes  difficult  to  pick  the  most  optimum  value  of  zeta 
and  omega  so  as  to  ensure  dominance  and  still  meet  the  system  specifica¬ 
tions.  The  graphical  solution  has  the  advantage  that  one  has  a  picture 
of  the  way  the  roots  of  the  characteristic  equation  move  around  in  the 
S -plane  as  alpha  and  beta  are  varied.  This  enables  one  to  pick  the  values 
of  alpha  and  beta  corresponding  to  the  best  values  of  zeta,  omega,  sigma, 
and  zeta-omega  for  all  the  roots  of  the  characteristic  equation.  This 
latter  feature  of  the  parameter  plane  points  out  a  strong  argument  for 
trying  to  obtain  the  parameter  plane  curves.  If  a  digital  computer  is 
not  available  then  by  using  the  relationships  derived  in  section  (3-3) 
under  sketching  techniques,  along  with  a  desk  calculator  or  slide  rule, 
the  curves  can  be  plotted  with  some  labor.  Under  these  circumstances  it 
is  questionable  whether  the  algebraic  or  the  graphical  solution  would  be 
better.  Which  one  is  used  is  a  matter  of  personal  preference. 

2 

The  computer  program  presented  in  section  (6-1)  is  most  helpful  in 
reducing  this  labor. 
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3.  Application  of  the  parameter  plane  to  the  compensation  of  linear 

continuous  systems . 

3-1.  Algebraic  solution. 

In  this  section  it  will  be  assumed  that  the  system  performance  specifi¬ 
cations  have  been  given  in  terms  of  placing  a  pair  of  complex  roots  at  a 
specific  value  of  zeta  and  omega,  say  J  ^  and  w^,  with  the  error  coef¬ 
ficient  Ke  being  greater  than  or  equal  to  a  specified  value.  If  the  specifi¬ 
ed  location  of  the  roots  is  such  that  after  computation  of  the  necessary 
value  of  alpha  and  beta,  the  remaining  roots  of  the  characteristic  equa¬ 
tion  are  located  so  that  the  specified  roots  are  not  a  dominant  pair,  then 
either  a  different  value  of  zeta  and  omega  will  have  to  be  used  (possibly 
at  the  sacrifice  of  some  measure  of  the  system  performance),  or  a  differ¬ 
ent  method  of  compensation  will  have  to  be  used. 

In  section  (3-2)  a  method  is  presented  whereby  the  dominancy  specifi¬ 
cation  can  be  met  by  introducing  a  third  parameter. 

3-1-1.  Feedback  compensation. 

Figure  (3-1)  represents  a  unity  feedback  control  system.  In  order 
to  meet  the  system  specifications  a  feedback  compensator  H  will  have  to 
be  used.  Let 


G  =  K _  = 

e(S) 


K 


Sm  +  e  S 
m-1 


m- 1 , 


.+e,^ 


(3-1) 


Where  K  is  the  forward  path  gain  which  can  be  varied  and  e(S)  is  a  poly¬ 
nomial  in  S  representing  the  poles  of  the  open  loop  transfer  function  of 
the  uncompensated  system.  The  letter  L  in  equation  (3-1)  corresponds  to 
the  system  type.  For  a  type  0  system,  L  =  0,  for  type  1,  L  =  1,  for  type 
2,  L  =  2,  etc.  By  definition,  the  error  coefficient  is  given  as  follows: 
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Figure  3-1 


(3-2) 


K 

e 


lim  SLG 
S-*0 


cc 


Here  G cc  is  the  open  loop  transfer  function  of  the  compensated  system. 
Sometimes,  as  for  example  in  reference  (8)  the  error  coefficient  is  design 
ated  as  for  a  type  zero  system,  for  a  type  one  system,  and  for  a 
type  two  system. 

Tachometer  plus  acceleration  feedback. 

Here  one  lets 

H  =  K  S  +  K  S2  (3-3) 

t  a  v  ' 

The  resulting  compensated  system's  characteristic  equation  becomes: 

e(S)+K+(K  S+K  S2)  =  0  (3-4) 

t  a 

or  by  expanding  e(S),  equation  (3-4)  becomes: 

Sm  +  e  1Sm“1  +. ..+(e0+KK  )S2+(e.+KK  )S+e  +K  =  0  (3-5) 

m-i  z  a  l  t  o 

where  L  is  taken  to  be  zero  for  a  type  zero  system  which  one  can  consider 
as  the  most  general  case.  The  following  results  also  apply  to  a  type 
one  system  if  e^  is  set  to  zero,  or  to  a  type  two  system  if  both  e^  and 
e^  are  set  to  zero,  etc.  From  equations  (3-2)  and  (3-4)  the  error  co¬ 
efficient  becomes: 

K  =  lira  5°  K  =  k  (3-6) 

e  S-*o  e(S)+K(KtS+KaS2)  eQ 

or  if  the  uncompensated  system  is  type  one: 


K 

e 


K 


ei+KKt 


and  if  the  uncompensated  system  is  type  two: 


(3-7) 


K  = 
e 


(3-8) 


Note:  If  the  uncompensated  system  is  type  two,  the  compensated  system 
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would  be  typo  one  if  tachometer  feedback  or  tachometer  plus  acceleration 
feedback  is  used. 

In  the  compensated  system* s  characteristic  equation  (3-5)  let  alpha 

=  KK  and  beta  =  KK  .  Equation  (3-5)  then  becomes: 

a  t  * 


Sm  +  e  .sm_1+.  .  ,+(e„+  oOS2  +  (e  +  B  )S  +  e  +K  =  0  (3-9) 

m-i  Z  l*o 

Recalling  from  equation  (2-8)  that  in  general  the  coefficients  of  the 
characteristic  equation  are  of  the  form: 

\  =  +  ck  P  +  ^ k 9  am*  corrcsPondingIy-  *h 

equation  (3-9),  m  =  k,  one  finds  that: 

d  =  e  -r  k ,  b  =  c  =  0,  d,  =  e,  ,  b„  =  0,  c.  =  1,  d.  =  e  ,  b.  =  1, 
o  o  *  o  o  *1  1’  1  *1  2  2  2 

cn=  0,  e  .  =  d.  -  ,  b.  -  =0,0.  .  =  0 ,  etc, 

2  m-1  k-1  k-1  k-1  9 


Then  from  equations  (2-10)  one  obtains: 


2  2  2 
Bj  =  (-1)  w  U1  =  w 


B2  -  w  u2 


C,  = 


D,  = 


-wU  =  0 
o 

m 


j <=0  K  K  1 


C2  =  *wUl  =  "w  (3-10) 

^  ^  ,  ixk  k¥T 

D2  ■  io  <‘1>  dk“  °k 


Use  was  made  of  the  fact  that  U  ^  -  -1,  Uq  =  0,  and  =  1  (see  appen¬ 
dix  I-B).  Using  the  expressions  for  alpha  and  beta  as  given  in  equations 
(2-11)  along  with  the  above  information  the  following  relations  evolve: 


cL  -  CiVC2D1  = 


B1C2-B2C1 


-  (-l)Vw1 

K-o  k 


k-2. 


U 


k-1 


(3-11) 


-w 


$c  <-«kdk”k'1<Uk-U2Uk-l) 


At  this  point  alpha  and  beta  may  be  linear  functions  of  K,  the  for¬ 
ward  path  gain,  and  one  can  use  the  steady  state  error  specification  to 
put  K  in  terms  of  alpha  and  or  beta.  Since  zeta  and  omega  were  assumed 
to  be  specified,  then  equations  (3-11)  can  be  used  to  solve  for  alpha  and 
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0 


\.  , 


* 


* 

r 


rr 

t,  i 


CJ 

I 


pt«' 


r  0  . 


=■•  .  I 


•>: 


>•  i 

f 

? 


12 


beta.  Then  K  and  K  can  readily  be  determined. 

a  u. 

Example  3-1 

The  system  given  in  figure  (3-2)  is  to  be  compensated  using  tacho¬ 
meter  plus  acceleration  feedback.  The  system  specifications  are  as 
follows: 

1.  Complex  roots  at  =  .7,  w  =  10. 

2.  K  —  6,  not  to  be  reduced. 

e  * 

Solution: 

From  equation  (3-2): 

K  =  lira  SLG  =  K  ^6 
6  S*o  CC  2  +  KKt 

or  K  — 12  -f  6KK^ .  The  compensated  system’s  characteristic  equation  is: 

S3  +  s2(3  +  KK  )  +  S(2  +  KK  )  +  K  =  0  (3-12) 

a  t 

Letting  alpha  =  KK  and  beta  =  KK  equation  (3-12)  becomes: 

a  t 

s3  +  S?'(3  4-  oC  )  +  S(2  +  £  )  +  K  =  0  (3-12a) 

Employing  equations  (3-10)  the  following  can  be  obtained: 

bl  =  ioo  cx  =  0 

Dx  =  -1100-K  B2  =  140 

C2  =  -10  D2  =  -1120 

Using  equations(2-ll)  one  obtains: 

ot  =  10 (-1100-K)  B  =  140 (- 1100-K)  +  56000  (3-13) 

-1000  -1000 

Note:  The  preceding  expressions  could  have  been  arrived  at  directly  by 
employing  equations  (3-11).  From  the  steady  state  accuracy  specifica¬ 
tions  it  is  necessary  that: 

K  -  12  +6^3  ,  hence  let  K  =  12+6  {9  (3-14) 


13 


Using  equation  (3-14)  and  (3-13)  beta  is  found  to  be: 

p  =  140 (-1100 -12 -6  B  )  +  56000  =  625 
-1000 

Therefore  K  =  12  +  6(625)  =  3762  and  oL=  .01(1100  +  3762)  =  48.6 

Since  ol  =  KK  then  K  =  48.6  =  .0129 
a  3  3762 

Also  from  3  =  KK  it  is  seen  that  K  =  .166. 

'  t  t 

The  compensated  system's  characteristic  equation  becomes: 

S3  +  51.62S2  +  627S  +  3762  =  0  (3-15) 

2 

Now  zeta  =  .7  and  omega  =  10  corresponds  to  S  +  14S  +  100.  When  this 
quadratic  is  divided  out  of  equation  (3-14),  the  remainder  is  S  +  37.62. 
Hence  zeta-omega  of  the  desired  roots  =  7  <<37.62,  and  the  complex  roots 
are  dominant  so  the  problem  is  solved. 

Example  3-2 

The  problem  presented  in  example  (3-1)  will  now  be  solved  by  intro¬ 
ducing  the  error  specifications  at  the  beginning  of  the  solution  instead 
of  at  the  end.  Equation  (3-12)  is  as  follows: 

S3  +  S 2 (3  +  KKa)  +  S(2  +  KK  )  +  K  =  0 
Again  let  alpha  =  KKfl  and  beta  =  KKt.  Then  from  the  steady  state  error 
requirement  one  obtains: 

K  =  12  +  6(S  (3-16) 

Substituting  this  expression  for  K  into  equation  3-12  results  in: 

S3  +  S2(3  +  d(_)  +  S(2+.(3)  +  12  +  6  jS  =  0 

Therefore  b  =  0,  c  =  6,  d  =  12,  and  all  other  coefficients  are  as 
o  ’  o  o 

before.  Hence  from  equations  (3-10)  it  is  seen  that: 

B  =  100,  Cx  =  -6,  D1  =  -1112,  B2  =  140,  C2  =  -10,  D2  =  -560 

and  cC  =  -7760  =  48.6  £  =  625  K  =  3762 

-160  ' 
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The  above  result  agrees  with  example  (3-1). 

Note:  Equations  (3-11)  could  not  be  used  here  since  they  are  based  on 

applying  the  accuracy  specifications  at  the  end. 

Tachometer  feedback  only. 

Let  H  =  KtS  (3-17) 

The  characteristic  equation  of  the  compensated  system  becomes: 

Sm  +  em-1  +  o2S2  +(e1+ot)S  +  e  +  |3  =  0  (3-20) 

o 

From  equations  (2-10)  one  obtains: 


B 


1 


=  0 


c 


1 


D 


1 


M)Wi 


B2  =  -W 


c2  =  0 

D„  = 


(3-21) 


K=o 


(-l)kd,  wV 


From  equations  (2-11)  it  is  seen  that: 

W\  |  If  1  i  i 

=  $.0  (_1)  V  ’  uk  and  f5  =  tio  (-D  dkwVi  <3‘22> 

If  a  given  zeta  and  omega  are  specified,  the  alpha  and  beta  can  be 
computed  from  equations  (3-22).  The  error  coefficient  is  then  determined 
directly  from  equations  (3-6),  (3-7),  or  (3-8).  Hence  the  error  coeffi¬ 
cient  is  fixed  once  zeta  and  omega  have  been  chosen,  so  if  a  certain 
error  specification  is  to  be  met,  the  specified  values  of  zeta  and  omega 
may  have  to  be  adjusted  to  meet  it. 

If  the  error  specification  was  the  overriding  specification  to  be 
met,  then  zeta  could  be  fixed  at  some  reasonable  value.  By  means  of  the 
given  Kg,  alpha  could  be  computed  from  equations  (3-6),  (3-7),  or  (3-8). 
Equations  (3-22)  could  then  be  used  to  solve  for  first  omega  and  then 
beta.  The  calculations  would  be  more  tedious  however. 

Example  3-3. 

Figure  (3-3)  shows  the  same  system  as  used  in  the  previous  two  ex¬ 
amples  only  now  tachometer  feedback  alone  will  be  tried.  The  same  system 
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specifications  are  to  be  met,  i.e.,  21  6,  zeta  =  .7,  and  omega  =  10. 

The  compensated  system's  characteristic  equation  then  becomes: 

s3  +  3S2  +  (2  +  KKt)S  +  K  =  0  (3-23) 

Using  equations  (3-19)  and  (3-23)  one  obtains: 

S3  +  3S2  +  (2+o^l  )S  +jS  =  0  (3-24) 

From  equations  (3-22)  alpha  is  found  to  be: 


cL  =  -2  +  30(1.4) -100(.96)  =  -56 

Since  alpha  is  negative  it  is  seen  that  positive  tachometer  feedback  is 
required.  Since  the  coefficient  of  the  first  power  of  S  in  the  character¬ 
istic  equation  would  then  be  negative,  the  system  would  be  unstable. 

Hence  the  desired  system  specifications  can  not  be  met  with  tachometer 
feedback  alone. 


Acceleration  feedback  only. 

Let  H  =  K  S2 
a 


(3-25) 


The  characteristic  equation  of  the  compensated  system  then  becomes: 


Sn  +  e  .Sm_1+...+  (e.  +  KK  )S2  +  e.S  +  e  +  K  =  0 
m-1  2  a  1  o 

Let  alpha  =  KKa  and  beta  =  K. 

Then  from  equations  (3-26)  one  obtains: 

s"  -  VIs"'1  +---+  <e2  ^  >S2  +  elS  +  eo  +p  -  0 

Using  equations  (2-10)  and  (3-28)  results  in: 


(3-26) 

(3-27) 

(3-28) 


n  2„  2 

B1  =  W  Ui  =  w 


ci  ■  u-i  ’  -1 

M  k  k 
Dt  =  ^  (-1)  d  w  u 
1  <-o  K 


B2  '  "  U2 


c2  =  o 


(3-29) 


k-1 


D,  =  |  (-l)kd,  w\ 


K-O 


Solving  for  alpha  and  beta  results  in: 
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(3-30) 


Calculations  for  alpha,  beta,  and  the  error  coefficient  are  performed 
in  the  same  manner  as  with  the  preceding  tachometer  feedback  calcula¬ 
tions. 

Example  3-4 

The  system  of  examples  (1)  and  (2)  will  now  be  compensated  using 
acceleration  feedback  as  indicated  in  figure  (3-4).  As  before  -ET  6, 
zeta  =  .7,  omega  =  10.  Therefore  Kg  -  K/2  and  the  error  coefficient  is 
unaffected  by  the  acceleration  feedback.  Hence  one  can  choose  K  =  12  to 
meet  the  specifications.  The  compensated  system's  characteristic  equa¬ 
tion  then  becomes: 


S3  +  (3  +  KK  )S2  +  2S  +  K  «  0 
a 


(3-31) 


If  K  in  equation  (3-31)  is  set  to  its  prescribed  value  of  12,  only 
one  parameter  remains  and  the  parameter  plane  equations  produce  an  indeter¬ 
minate  solution.  Therefore  K  will  be  left  as  the  variable  beta.  Since 
beta  is  fixed  by  the  chosen  values  of  zeta  and  omega,  then  so  is  the  error 
coefficient,  and  it  will  most  likely  not  agree  with  the  error  specifica¬ 
tion.  In  view  of  this,  the  solution  proceeds  as  follows.  Making  the 
usual  change  of  variables  in  equation  (3-31)  results  in: 


(3-32) 


By  employing  equations  (3-30)  one  can  solve  for  alpha  and  beta. 


cL  =  -1  [  -2/10  +  3(1.4)  -10(.96) ]  =  4 


1.4 

(3  =  3(10) 2  - ( 10) 3 ( 1 . 4)  -U4  [-2(10)  +  3(10)2(1.4)-(10)3( .96) ] 
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From  the  negative  value  of  beta  it  is  concluded  that  the  desired 
values  of  zeta  and  omega  cannot  be  obtained  using  acceleration  feedback, 
and  of  course  neither  can  the  desired  error  specification  be  obtained. 

One  would  therefore  choose  another  method  of  compensation. 

If  one  chooses  to  use  feedback  compensation  then  perhaps  tachometer 
plus  acceleration  feedback  should  be  tried  first  using  equations  (3-11) 
and  the  appropriate  steady  state  error  specifications.  If  the  specifica¬ 
tions  cannot  be  met  in  this  manner  then  it  is  obvious  that  they  cannot  be 
met  by  either  tachometer  or  acceleration  feedback  separately.  In  this 
case  either  the  system's  specifications  must  be  modified  or  another 
scheme  of  compensation  must  be  employed.  If  it  is  found  that  the  specifi¬ 
cations  car.  be  met  by  combined  tachometer  and  acceleration  feedback,  then 
if  desired,  equations  (3-22)  or  (3-30)  can  be  used  to  see  if  tachometer  or 
acceleration  feedback  alone  can  be  used. 

Cases  where  feedback  is  not  around  the  forward  path  amplifier. 

Figure  (3-5)  illustrates  a  compensation  situation  that  sometimes  occurs 
in  practice.  This  is  the  situation  where  it  is  not  possible  or  practical 
to  get  at  the  input  terminals  of  the  error  detector  and  the  feedback  has 
to  be  inserted  at  the  output  terminals  of  the  amplifier.  The  solution  to 
this  problem  is  solved  by  means  of  an  example. 

Example  3-5 

Figure  (3-5)  shows  the  system  that  was  used  in  example  (3-1)  only 
now  the  feedback  is  inserted  at  the  output  terminals  of  the  amplifier 
represented  by  the  gain  K.  The  same  system  specifications  are  to  be  met 
i.e.,  —  6,  zeta  =  .7,  and  omega  =  10.  The  characteristic  equation  now 

becomes : 

S3  +  (3  +  K  )S2  +  (2  +  K  )S  +  K  =  0  (3-33) 

C 
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Letting  alpha  =  K  and  beta  =  K  equation  (3-33)  becomes: 
a  t 

S  +  (3  +cC  )S^  +  (2  +  ^  )S  +  K  =  0 


(3-35) 


Comparing  equation  (3-35)  to  equation  (3-12a)  it  is  seen  that  they  are 
identical,  so  the  solution  obtained  for  alpha  and  beta  in  example  (3-1) 
applies.  This  points  out  an  important  advantage  of  the  parameter  plane 
method.  This  is  that  the  solutions  depend  only  on  the  characteristic 
equation  and  not  on  the.  system  that  the  characteristic  equation  was 
formed  from. 

From  example  (3-1)  it  was  found  that  alpha  =  48.6,  beta  =  625,  and  K 
=  3762.  In  this  example  there  are  no  additional  computations  necessary 
to  find  K  and  K  ,  since  these  are  now  the  system's  parameters  alpha  and 

3.  t 

beta.  So  K.  =  625  and  K  =48.6. 
t  a 

This  same  general  principal  can  be  applied  to  control  problems  in¬ 
volving  tachometer  feedback  alone  or  acceleration  feedback  alone. 

3-1-2  Cascade  compensation. 

Figure  (3-6)  represents  a  unity  feedback  control  system  which  in 
order  to  meet  the  system's  specifications  a  cascade  compensator  Gc  is  re¬ 
quired.  Let  G  have  the  form  of  equation  (3-1). 

i.e. ,  G  =  K _  =  _ K _  (3-1) 


e(S) 


Sm  +  e  1Sm_1+...+e  SL 
m-1  l 


K  is  the  forward  path  gain  which  can  be  varied  and  e(S)  is  a  polynomial 
is  S  representing  the  poles  of  the  open  loop  transfer  function  of  the  un¬ 
compensated  system.  The  letter  L  again  corresponds  to  the  system  type. 


Le  t : 


G  =  P(S  4-  Z)  (3-36) 

C  Z(S  +  P) 

This  compensator  has  a  D.C.  gain  of  unity  so  its  placement  in  the  forward 
path  will  not  affect  the  steady  state  accuracy.  It  is  assumed  that  the 
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uncompensated  system's  forward  path  gain  had  previously  been  adjusted  to 
give  the  correct  steady  state  accuracy.  Using  the  form  of  G as  indicated, 
the  values  of  Z  and  P  are  computed  to  give  the  desired  system  response. 

If  P  is  less  than  Z  a  lag  network  is  needed  and  the  factor  P/Z  of  the 
compensator  is  inherently  present  due  to  the  physical  nature  of  the  compen¬ 
sator,  which  is  assumed  to  be  an  R-C  network.  See  appendix  II  concerning 
the  nature  of  lag-lead  R-C  networks.  In  this  case  all  forward  path  ampli¬ 
fier  gains  can  remain  unchanged  to  meet  the  stipulated  accuracy  specifica¬ 
tions.  If  however,  the  computed  values  of  Z  and  P  are  such  that  P  is 
greater  than  Z,  a  lead  network  is  required  and  the  compensated  system's 
forward  path  gain  will  have  to  be  raised  by  a  factor  of  P/Z.  The  physical 
nature  of  the  lead  network  is  such  that  the  factor  P/Z  is  not  inherently 
present,  so  to  maintain  steady  state  accuracy  this  factor  will  have  to 
be  provided  either  by  adding  an  amplifier  in  cascade  with  the  lead  net¬ 
work  or  by  raising  the  gain  K  of  the  existing  amplifier  by  this  factor. 

Continuing  then,  the  procedure  is  as  follows.  The  compensated  system's 
forward  path  transfer  function  is: 

G  =  G  G  =  K  P  S  +  Z  =  TSfs  +  PA)  _ K  (3-37) 

CC  C  e(S)  *  Z  *  S  +  P  S  +  P  '  e(S) 

Applying  the  definition  of  the  error  coefficient  to  the  compensated  system 

one  obtains: 

K  =  lim  SLf  K  ,  ^(S  +  P/Tf  )  j  =  K  (3-38) 

e  S-X)  [_e(  S)  (S  +  P)  J  eL 

Again  assuming  a  type  zero  system  where  L  =  0,  the  compensated  system's 
characteristic  equation  becomes: 

e(S)  (S  +  P)  +  K  (S  +  P/T  )  =  0  (3-39) 

or  in  general  after  expanding  equation  (3-39): 
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(3-40) 


Sm+1  +  (P  +  e  )Sm  +  (Pe  +  e  )Sm"1+...+  (Pe  +  e  )S2  + 
m- 1  m-1  m-z  Z  1 

(KX  +  Pe.  +  e  )S  +  P(e  +  K)  -  0 
l  o  o 

Letting  alpha  =  P  and  beta  =  K  ,  equation  (3-40)  becomes: 

Sm+1  +  ((/-+  e  )Sm  (e  ,oL+  e  +. .  .+(e,  oC  +  e.)S2  (3-41) 

m-i  m-i  m-Z  Z  1 

(K  {3  +  ex  oL  +  eQ)S  +  °^-(e0  +  K)  =0 

Comparing  equation  (3-41)  to  the  general  form  of  the  characteristic  equa¬ 
tion  as  specified  in  equations  (2-1)  and  (2-8),  it  is  apparent  that  K  = 
m  +  1  (the  order  of  the  equation),  bQ  =  eQ  +  K,  cq  =  dQ  =  0,  b^  =  e^, 
c^  =  K,  d1  =  eQ,  b2  =  e2,  c^  =  0,  d2  =  ej^,  etc. 

It  is  important  to  note  that  the  parameter  plane  variable  beta  re-  . 
presents  the  pole  to  zero  ratio  of  the  cascade  compensator.  In  references 
(1)  and  (2),  Ross-Warren  and  Poliak  respectively,  utilized  the  concept  of 
root  relocation  zones  to  divide  the  S-plane  into  regions  where  lag  compen¬ 
sation  or  lead  compensation  is  required.  By  assigning  variables  in  the 
above  manner,  the  parameter  plane  is  effectively  divided  into  corresponding 
regions  above  and  below  the  straight  line  beta  =  1.  For  values  of  beta 
less  than  one,  a  lag  network  is  required  and  for  values  of  beta  greater  than 
one,  a  lead  network  is  required.  In  addition  if  beta  is  greater  than  10 

f 

or  less  than  .1,  a  multiple  le^dor  multiple  1®.^  respectively  is  required. 

A  multiple  section  compensator  will  also  be  required  if  the  computed  value 
of  either  P  or  Z  turns  out  to  be  negative.  A  method  is  given  in  section 
(4-6)  for  the  design  of  double  section  compensators.  If  more  than  two 
sections  are  required,  the  compensation  has  to  be  done  in  steps. 

In  this  case  complex  roots  are  placed  at  some  intermediate  value  of 
zeta  and  omega  and  a  new  characteristic  equation  is  then  computed.  An¬ 
other  compensator  can  then  be  designed  on  the  basis  of  this  new  character¬ 
istic  equation  to  place  the  roots  at  the  desired  value.  Thus  by  compensa- 
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tion  in  steps  in  conjunction  with  the  double  section  design  of  section 
(4-0),  a  four  section  compensator  could  theoretically  be  designed.  The 
use  of  more  than  four  sections  is  questionable  and  in  this  event  it  would 
be  better  to  employ  combined  cascade  and  feedback  compensation,  or  per¬ 
haps  feedback  compensation  alone. 

On  the  basis  of  equations  (3-41)  and  (2-10)  it  is  found  that: 

B1  =  ~(eo  +  K)  +  w  e2  +•••+("!)  w  X-3  +  W  Uk-2 

C1  =  0 

p>  2  ,  ,/•  ,\k-2  k-2_t  ,  t  ..k-1  k-l„  -l 

D1  =  W  el  •  •+(’1)  em-2W  Uk-3  +  <’1>  em-lW  Uk-2  + 

(-1)^*0^  (3-42) 

B2  =  -wei  +  w2e2U2  +...+(-l)k_2Uk_2  +  (-l)k'1wk'1Uk_1 

C2  =  -wK 

_  -  2  „  .  .  .  k-2  k-2,. 

D2  =  -weo  +  «  ejU2  +. .  .+(-1)  em_2w  \_2 

and  from  equations  (2-11)  it  is  found  that: 


B. 


B2DrBiD2 

-wKB. 


(3-43) 


In  equations  (3-42),  one  sets  eQ  =  0  for  a  type  one  system,  eQ  =  e^  =  0 
for  a  type  two  system,  etc.  On  the  basis  of  equations  (3-42)  and  (3-43) 
a  cascade  compensator  can  be  designed. 

Example  3-6 

Problem: 

Design  a  cascade  compensator  for  the  system  shown  in  figure  (3-7) 

to  place  a  pair  of  characteristic  roots  at  zeta  =  .5  and  omega  =  1.  The 

error  coefficient  K  should  be  50. 

e 
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Solution: 


From  figure  (3-7)  it  is  apparent  that  =  K/10  =  50  or  K  =  500.  The 
characteristic  equation  is: 

S4  +  S3(8  +  P)  +  S2 (17  +  8P)  +  S(10  +  17P  +  KY )  +  P(10  +  K)  -  0  (3-44) 

Letting  alpha  =  P  and  beta  =  If, equation  (3-44)  becomes: 

S4  +  S3(8  +o C)  +  s2(17  +  8o4)  +s(10  +  17oC+  500.0  )  +  (10  +  500)c4  =  0(3-45) 
Applying  equations  (3-42)  and  (3-43)  one  obtains: 


=  -503 


=  .0179  =  P 

But  =  P/Z  so  Z  =  1.529. 

G  then  becomes:  .0117  (S  -u  1.529) 

C  (S  +  .0179) 


C2  =  -500 
D2  =  6 

|S  =  .0117  =  % 


This  is  a  lag  network  where  the  factor  .0117  is  inherent  in  the  R-C  filter 
design.  Due  to  the  small  value  for  gamma,  the  size  of  the  filter  compon¬ 
ents  may  be  unreasonable,  and  a  double  lag  network  could  be  designed  us¬ 
ing  the  method  of  section  (4-6-2).  Instead  the  problem  will  be  solved  in 
section  (3-1-3)  using  combination  cascade  plus  tachometer  feedback  compen¬ 
sation  which  will  permit  the  use  of  a  single  section  lag  network. 

3-1-3  Combination  cascade  and  feedback  compensation. 

Derivative  signal  not  enclosing  a  cascade  compensator. 

In  figure  (3-8),  G  =  K/e(S)  and  G  =  ^  S  +  P  (3-46) 

c  sTT 

First  let  H(S)  =  KfcS.  Then  in  view  of  equations  (3-46)  the  compensated 
system's  forward  path  transfer  function  becomes: 


G  15  S  +  P  .  K/e(S) _  =  y  S  +  P  .  _ K _  (3-47) 

CC  S  +  P  1  +  KH(S)/e(S)  S  +  P  e(S)  +  KH(S) 
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i 


Figure  3-8 


(3-48) 


and  K  =  lira  SL  G  = _ _ K _ 

r;*0  CC  e  (!J)  +'  lira  KH(S)/SX 

S^O 

Note:  The  system  type  will  change  if  the  lowest  order  of  the  deriva¬ 
tive  signal  fed  back  is  less  than  the  system  type  number.  A  few  error  co¬ 
efficients  are  given  in  table  (3-1)  for  different  types  of  uncompensated 
systems.  In  the  table  let  L  be  the  type  number  of  the  uncompensated 


system. 

Table  of  K 

e 

L  — > 

0 

1 

2 

H(S)^ 

KtS 

K/e 

o 

K/(e1+KK1) 

1/KX 

K  S2 
a 

K/e 

o 

K/ej 

K/(e2+KK2) 

Table  (3-1) 


The  compensated  system's  characteristic  equation  is: 

Se(S)  +  Pe(S)  +SKH(S)  +  PKH(S)  +  K  Y  S  +  PK  =  0  (3-49) 

Using  equation  (3-1)  for  e(S)  and  letting  alpha  =  P,  beta  =  "ft  ,  and 
k  -  m+1,  equation  (3-49)  becomes: 

+  e  1)S^C^’  +  (e  .c4+e  _)S^  ^  H-. .  ,+(e.  +  KK  4-  e-  cC.  )S 

m- 1  m- 1  m-Z  l  z  z 

+  (e,  oC.  +  KKoC+  K/3  +  e  )S  +o/  (K  +  e  )  =  0  (3-50) 

1  z  «  o  o 

The  parameter  plane  variables  are  then: 

»1  ‘  -%'K+  +-"+<-1>k'2Vl”k'\-3+  <-l)k-1wk-\.2 


Cl.‘° 


2- . k"2  ,„k-V 


Dj  ■  +  “t)»  e2w  »2+- ••+<-!>  em-2"  -k-3 

+  (-1)k’lem-l”k’luk-2  +  <-1)k”k°k-l 


(3-51) 
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B2  13  '(G1  '!‘  KJ<t)w  '■  e2V/2lJ2  +---+(-Dk"2Vlwk’\-2 


,  .  ,  .k-1  k-1 

-I-  (-1.)  w  U 


k-1 


C2  =  -Kw 

D2  =  -eQw  (e],  +  KKt)w2U2  +. .  ,+(-l)k"2ein_2wk-2Uk_2 

,  >  -vk-1  k-l„  ,  .  ,.k  k„ 

+  <-l>  e„-i”  Vi  +  (-1)  “  uk 

In  tenns  of  equations  (3-51)  one  can  solve  for  alpha  and  beta: 

cf—  ~  -D,  /B  @  =  B2D1  "  B1D2  (3-52) 

-wKB1 

One  can  now  compare  the  above  expressions  for  alpha  and  beta  with 
equations  (3-42)  and  (3-43)  to  see  the  effect  of  tachometer  feedback. 

Let  Bj,  C|,  Dj,  B2,  C2>  and  D2  represent  the  quantities  given  by  equations 
(3-42)  where  only  cascade  compensation  was  used.  Then  in  terms  of  the 
primed  quantities,  equations  (3-51)  become: 


B1  =  Bi  B2  -  B2  -  “V 

Cj_  =  CJ  =  0  c2  =  c2  (3-53) 

Dx  »  DJ  +  KKtw2  D2  =  +  KKfcw2U2 

Let  the  expressions  for  alpha  and  beta  as  given  by  equations  (3-43) 
be  designated  JL  '  and  Q  Then  in  terms  of  these  quantities,  equa¬ 
tions  (3-52)  can  be  expressed  as: 

o£-  =  (-DJ  -  KKtw2)/Bj  =o4  '  -  KKtw2/Bj  (3-54) 

(3  =  (-B^wKt  +  KtDj  +  K2Kw2)/B|  +  KtwU2  +  g  * 

Now  it  can  be  seen  from  equations  (3-54)  how  tachometer  feedback  modifies 
the  values  of  cA  '  and  p,  1  as  computed  for  cascade  compensation  alone. 

For  instance  if  the  pole  to  zero  ratio  6  '  is  too  small,  then  tachometer 
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feedback  can  be  used  to  increase  this  ratio,  if  for  the  specified  values 

of  zeta  and  omega,  Bj  is  negative,  etc. 

2 

Now  letting  H(S)  =  K  S  in  figure  (3-8),  the  characteristic  equa- 

a 

tion  becomes: 


Se(S)  +  Pe(S)  +  KK  S3  +  Pkk  S2  +  K  ''S  S  +  PK  »  0 

cl  3 


(3-52) 


By  making  the  same  substitutions  as  in  equation  (3-49)  one  obtains: 

S^+(c4+e  -)S  *+(e  0)S^  ^  +• • .+(e0  +  e~oC  +  kk  )S^ 

m- 1  m-i  m-z  z  J  a 

(ex  +  e2  o4  +kkaoC)S2  +  (el0^  +  K  ^  +  eQ)S  +  o£  (K  +  eQ)  =  0  (3-53) 

It  then  follows  readily  that: 


B,  =  -(K  -r  e  )  +  (e  +  KK  )w2  -e„w3U0  +.  ,.+(-l)k  2e_  ,wk-2U. 


.  ,  ,xk-l  k-1 
+  (-1)  w  U. 


k-2 


m-1 


k-3 


ci-° 


2  . 3 . k-2 


D1  =  elW  “  (e2+KKa^W  em_2 

+<-1>k\.lwk'\-2  + 


k-3 


B,  =  -ew+  (e.fKK  )w2U0  -  e,w3U,+. .  .+(-l)k_2e  1wk_2U. 

2  '  2  a  2  3  3  m-1  1 

t  , .k-1  k-1 
(-1)  w  0k.! 


k-2 


(3-53a) 


C2  =  -Kw 


.k-2 


D2  -  -V  +  elW‘U2  -  (e2+KKa)»-U3  +...+(-1)- * 
<-1)k‘lem-l”k‘luk-l  + 

After  writing  equations  (3-53a)  in  terms  of  the  primed  quantities,  which 
correspond  to  the  situation  of  cascade  compensation  only,  it  can  be  seen 
that  the  following  expressions  result: 


B,  =  Bi  +  KK  w 
11a 


B2  =  B*  +  KKw  U2 
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(3-54) 


D.  =  D'  -  KK  w 
11a 


3 


Solving  for  alpha  and  beta  one  obtains: 


oC=  -(D'  -  KK  w3)/ (B '  +  KK  w3) 

1  a  jl  3, 

P  =  (3^  +  KKaw2U2) (Dj  -  KKaw3)  -  (B  j  +  KK,w2) (D^  -  KKgw3U3) 


(3-55) 


-wK(Bj  +  KKflw2) 


There  is  no  straight  forward  way  of  observing  from  equations  (3-55) 
how  c4  1  and  '  are  modified  by  acceleration  feedback.  However,  assum¬ 
ing  that  one  would  logically  try  cascade  compensation  alone  before  attempt 
ing  combination,  one  could  easily  compute  the  alpha  and  beta  with  accelera 
tion  feedback  added  since  the  quantities  B|,  B2,  Dj,  and  D2  in  equations 
(3-55)  would  already  have  been  computed. 

Example  3-7 

Problem: 

Design  a  cascade  compensator  with  tachometer  feedback  as  shown  in 
figure  (3-9),  so  as  to  place  a  pair  of  characteristic  roots  at  zeta  =.5 
and  omega  =1.  Kg  should  be  50. 

Solution: 

Referring  to  example  (3-6)  it  is  seen  that  this  is  the  same  problem 
except  that  tachometer  feedback  has  been  added.  As  was  found  in  example 
(3-6)  the  value  of  beta  using  a  single  section  compensator  was  too  small. 
The  problem  now  is  to  see  if  tachometer  feedback  will  increase  this  value. 

As  was  determined  previously,  K  =  500,  Bj  =  -503,  Cj  =  0,  D|  =  9, 

B2  =  -9,  C2  =  -500,  D2  =  6,  '  =  .0179,  and  ^ 1  =  .0117,  where  the  primes 

have  been  added  to  indicate  cascade  compensation  only.  Using  equations 
(3-55)  one  finds: 


33 


u 


A5 


Figure  3-9 


(3-55a) 


|3  =  (9Kt  +  9Kt  +  500K2)/(-503)  +  Kfc  +  .0117 

It  follows  from  the  above  equation  that  for  beta  to  increase  the  following 
inequality  must  hold: 

Kfc  >  (18Kfc  +  K2)/503 

or  (3-56) 

K2  <  (485Kt)/500 

Since  one  excludes  values  of  less  than  or  equal  to  zero,  equation  (3-56) 
reduces  to: 

0  c  Kt  <  .97 

can  arbitrarily  be  taken  as  .4.  Beta  can  then  be  obtained  from  equa¬ 
tion  (3-55a)  and  is  found  to  be  .238  which  is  in  the  acceptable  region  of 
.1  2  £=*  5  10. 

Calculating  alpha  from  equations  (3-54)  one  finds  that: 

c>4  =  o4  '  -  KKtw2/Bj  =  .414  =  P 

The  cascade  compensator  becomes: 

.238  (S  -f  1.744)  ,  which  can  be  synthesized  as  an 
(S  +  .414) 

R-C  lag  network.  The  compensated  system's  characteristic  equation  is  as 
follows: 

S4  +  8.415S3  +  220. 3S2  +  219. 06S  +  211.69  =  0  (3-57) 

The  roots  of  equation  (3-57)  are: 

r^  =  -.499  -  j.866 

r  «  -.499  +  j.866 
2 

r3  =  -3.708  -  j 14. 077 
r4  =  -3.708  +  jl4.077 

The  roots  r^  and  r ^  are  the  desired  ones,  and  since  .499  3.708,  they 

are  also  dominant. 
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Derivative  signal  enclosing  a  cascade  compensator . 


In  figure  (3-10),  G  =  K/e(S)  and  G  =  If  S  +  P  ,  where  these 

C  S  +  P 

quantities  have  been  defined  previously.  Let  H(S)  =  K^S.  The  forward 
path  transfer  function  is  then  seen  to  be: 

G  =  GcG _  -  _ K  Y  (S  +  P/'S  ) _  (3-58) 

CC  1  +  H(S)GcG  (S  +  P)e(S)  +  K  "8  (S  +  P/*  )H(S) 

Since  the  D.  C.  gain  of  G£  is  unity,  the  error  coefficients  for  figure 

(3-10)  are  the  same  as  those  given  in  table  (3-1). 

By  expanding  equation  (3-58)  into  the  characteristic  equation  one 
gets  after  letting  alpha  =  P  and  beta  =  Y  : 

2 

Se(S)  +  Pe(S)  +  KY  KtS  +  KPKfcS  +  K  ft  S  +  KP  =  0,  or  after  expanding: 

Sk  +  (cJ-  +  e  )S  +  (e  .0^.+  e  „)Sk  2  +.  ,.+(e.  +  KK  8  +  Pe„o£-  )S2 

m-1  m-1  m-z  1  t  r  2 

+(el°^  1  +  KK  +  K  @  +  eQ)S  +  oY  (K  +  e  )  =0  (3-59) 


Then  as  before  one  obtains: 

B.  =  -(e  +  K)  +  w2e.  +. . .+(-l)k_2e  ,wk‘2U,  ,  +  (-1) 

1  o  2  m-1  k-3 

Cx  =  KKtw2 

_  2  .  . Nk-2  k-2  .  ,  ..k-1  k-1 

D1  "  elW  +...+(-!)  em_2w  Uk.3  +  (-l)  e^w  U 

+(-l)kwkUk_1 

B2  =  -(el  +  KKt)  w  +  e2w2u2+-*-+<-1)k’2em.iwk’2uk.2  + 
C2  =  -Kw  +  KKfcw2U2 

D2  .  -eo„  +  +  <-l)k‘2V2"k'2uk-2  +  <’I>k‘lV 

+  (-Dkw\ 


k-1  k-llt 
w  U, 


k-2 


k-2 


,  i .k-1  k-1 
(-1)  w 


k-1 

1W  \-l 


U 


k-1 


With  primed  quantities  corresponding  to  the  case  of  cascade  compensation 
only  it  follows  that: 
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Figure 


C,  =  C'  +  KK  w2  =  KK  w2 

lit  t 

di  =  d; 

B2  -  B2  -  KKtw 

C2  =  C‘  +  KKtw2U2 

D  =  D' 

2  2 

and 

d-  =  KKfcw2D^  -  <c2  +  KKtw2U2)Dj 

B|(C2  +  KKtw2U2)  -  (B2  -  KKtw)KKtw2 


(3-60) 


(3-61) 


P  ■  <B2  -  KKtw)Di  '  B1D2 


B|(C^  +  KKfcw2U2)  -  (B2  -  KKtw)KKtw2 


In  figure  (3-10)  letting  H  =  K  S  ,  the  characteristic  equation  becomes: 

cl 

sk  +  (oL+  em_1)Sk'1  +  (em_1oC  +  em_2)sk’2  +**‘+(e2  +  e3^  +  ^a 
+  (ex  +  e2o(.  +  KKa«l  )S2  +  (e^  +  K  0  +  eQ)S  +  o£  (K  +  e  )  =0  (3-62) 

Then  in  terms  of  the  primed  quantities  it  is  found  that: 

B.  =  B.'  +  KK  w2 
11a 

C.  =  C'  -  KK  w3U0  =  -KK  w3U0 

11  a  2  a  2 


di  =  d; 


(3-63) 


B2  =  B2  +  KKaW  U2 


C2  =  C2  -  «aw'  U3 


D  =  D ' 
2  2 

and 


c/.  =  -KK  w3U0D1  -  (Ci  -  KK  w  u0)D ' 

a  2  2  2  a  3  i 


(B[  +  KKaW2)(C»  -  KKaw3U3)  +  (B^  +  KKw2U2)KKw3U2 


(3-64) 
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(3  =  +  KKaw2U2)Dj  -  (Bj  +  KKaw2)D^ 

(B|  +  KKaW2)(C^  -  KKaw3U3)  +  (B^  +  KKaw2U2)KKaw3U2 

Comparing  the  expressions  for  alpha  and  beta  for  the  case  of  the  deriva¬ 
tive  signal  not  enclosing  a  cascade  compensator  to  the  case  of  the  de¬ 
rivative  signal  enclosing  a  cascade  compensator  one  finds  that  the  latter 
are  considerably  more  complex.  So  if  one  had  a  choice  between  the  two 
methods,  the  former  could  be  tried  first  since  it  is  easier  to  analyze. 

If  the  values  for  alpha  and  beta  thus  obtained  were  still  not  acceptable 
then  the  latter  method  could  be  tried. 

3-2  Dominancy  of  the  specified  roots. 

In  the  preceding  examples  nothing  was  done  in  the  calculations  to 
make  the  specified  roots  a  dominant  pair.  As  was  mentioned  in  section  (1), 
being  able  to  predict  a  system's  response  on  the  basis  of  the  location  of  a 
pair  of  complex  roots  was  based  on  the  assumption  that  the  magnitude  of 
the  real  part  of  the  primary  or  specified  roots  was  much  less  than  the 
magnitude  of  the  real  parts  of  all  the  other  roots  of  the  characteristic 
equation.  In  most  cases,  if  the  real  part  of  the  primary  roots  is  one 
half  to  one  fifth  or  less  of  the  real  parts  of  all  the  secondary  roots, 
the  system  is  said  to  be  dominant  in  the  primary  roots.  In  many  cases 
the  system  will  still  meet  the  specifications  even  if  two  pairs  of  complex 
roots  have  the  same  real  part,  providing  the  zetas  for  both  pairs  of  roots 
meet  the  specifications,  and  the  undamped  natural  frequencies  are  such 
that  the  component  time  responses  are  not  highly  additive.  The  presence 
of  closed  loop  zeros  will  also  greatly  affect  the  dominancy  factor  needed. 
For  instance  even  if  there  exists  a  characteristic  root  whose  real  part  is 
closer  to  the  origin  than  the  real  part  of  the  primary  root,  the  presence 
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of  a  closed  loop  zero  could  make  the  residue  of  the  close  in  root  negli¬ 
gible  as  compared  to  the  residue  of  the  parimary  roots.  However,  if  possi 
ble,  one  tries  to  make  the  real  parts  of  all  secondary  roots  as  large  in 
magnitude  as  possible. 

In  the  preceding  examples  it  should  be  noted  that  in  many  cases 
there  were  actually  three  and  sometimes  four  variable  parameters.  For 
instance,  the  forward  path  gain  was  usually  set  at  a  fixed  value  in  the 
computations  so  as  to  meet  the  minimum  steady  state  accuracy  requirements. 
There  is,  however,  usually  no  reason  why  the  gain  cannot  be  raised  above 
the  minimum  value,  thus  permitting  a  third  degree  of  freedom.  When  cas¬ 
cade  and  feedback  compensation  are  employed  simultaneously,  the  forward 
path  gain  and  tachometer  gain  become  the  third  and  fourth  parameters. 

3-2-1  A  method  of  employing  a  third  parameter. 

Recall  that  the  system  characteristic  equation  has  the  following  form 


(2-8)).  In  order  to  meet  the  system  specifications,  one  places  a  com- 


2 

1-^j  ,  which  implies  that 


(3-65) 


Since  the  coefficients  of  equation  (3-65)  are  known,  this  quadratic  can 
be  divided  out  of  the  characteristic  equation,  leaving  a  polynomial  which 
contains  all  the  secondary  roots  of  the  characteristic  equation.  Since 
only  two  of  the  degrees  of  freedom  or  variable  parameters  were  used  in 
fixing  the  roots  of  equation  (3-65),  the  remaining  variable  parameters 
will  appear  in  the  coefficients  of  the  quotient  polynomial,  and  it  is 
these  coefficients  that  can  be  varied  to  achieve  dominance.  Instead  of 
division  to  find  the  quotient  polynomial,  coefficients  of  like  powers  will 
be  equated  to  achieve  a  set  of  equations.  Let  the  quotient  polynomial  be 
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given  by: 


fx(S)  = 


v“  -  • 


(3-66) 


where  n  =  m-2,  i.e.,  equation  (3-66)  is  order  two  less  than  the  character¬ 
istic  equation.  Using  equations  (2-1),  (3-65),  and  (3-66)  it  is  seen 
that: 


(3-67) 


Equating  coefficients  of  like  powers  and  taking  a^  =  1,  results  in: 

a  =  f  =1 
k  n 

\-l  =  fn-l  +  2  7  1W1 

2 

a,  0  =  f  0  +  f  ,2 
k-2  n-2  n-1  7 11 


2  (3-68) 

a_  =  f  +2  ^  i w  f  +  f-w. 

2  o  H  1  1  21 

®1  =  2  }  lVo  +  fl6l 

a  =  f  w2 
o  o  1 

The  formulas  (3-67)  can  be  solved  for  the  coefficients  f  in  terms  of  the 
coefficients  a.  The  solution  will  be  made  for  the  following  cases: 

Case  of  k  =  3,  n  =  1 
Equation  (3-67)  becomes: 


(S2  +  2  j  +  wJ)(f1S  +  fQ)  =  S3  +  a2S2  +  aLS  +  aQ 
Equating  coefficients  of  like  powers  one  obtains: 
a3  =  1  “  fl 
a2  =  fo  +  2  7  l6lfl 

ai =  fiwi  +  2  7  iVo 

Solving  for  the  coefficients  f  results  in: 
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1 


£1  - 


£o  ■  ao/wl  =  <al  ?lV  =  a2-2  }  l“l 

Case  of  k  =  4f  n  =  2 

Proceeding  as  before  the  a  coefficients  become: 
a4  =  1  =  f2 
a3  =  2  f  1W1  +  fl 
a2 
al 


*7  + 

2  ? 

,w,  f  + 

1 

? 

111 

,  2 
f„  w. 

+  2 

?  _w,  f 

1  1 

f  w* 
o  1 

(1  1  o 

When  solved  for  the  coefficients  f,  equations  (3-70)  yield: 

£2-  1 

£1  ”  S3  ‘  2  f  l”l  ”  a2/2  ^  l”l  '  ”l/2  /  1  '  ao/2  /  l”l 
£1  =  1/wl  <al  -  2  ?la„/ui> 

£„  ’  V“l  *  al/2  ?  l“l  -  Vl/2  ?  1  +  W1  *  a2  •  2  ?  l"la3 

w2  +  4  ^i»2-a2  -  2  ^  lwla3  +  wl*^3 
Case  of  k  =  5,  n  =  3 

Proceeding  as  before  the  coefficients  a  are: 


a5  -  1  -  £3 


°4  =  2  ?1W1  +  f; 


a3  =  W1  +  2  J  Iwif2  +  fl 


2 

~f  Iwifl  +  wlf2 


an  =  f  +  2  ?  ,w,  f,  +  w,  f 
2  o 


a  =  2  ?  ,e. f  +  f  ,w 
1  /  1  1  o  1  i 


a  =  f  w, 
o  o  1 


(3-69) 


(3-70) 


(3-71) 
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Solving  for  the  coefficients  f  one  obtains: 


f2  '  a4  •  2  ?1»1  *  S2/Wl  "  *o(4  ?l/ul  -  1/wl>  -  2  /  lal/Ml 

f2  =  a3/2  ?  l”l  '  V2  ?  l“l  +  ao/ut  •  ul/2  ?  1 

fl  al/wl  -  2  }  lao/Ml  *  a3  -  2  ?  l”la4  +  U3”l  (3‘73) 

£o  ‘  ao/ul  "  al/2  ?  1U1  ‘  a3wl /2  ?  1  +  Vl  +  s  ?  l"l 
£o  ‘  a2  -  a32  }  l”l 

In  fomulas  (3-69),  (3-71),  and  (3-73),  the  coefficients  a  are  of  the 
form  a^  =  +  c^  ^  +  d^.  In  section  (3-1),  formulas  for  alpha  and 

beta  were  developed  for  various  forms  of  compensation  techniques.  Since 
the  formulas  for  alpha  and  beta  will  in  general  contain  other  variable 
parameters,  then  the  coefficients  f  as  derived  above  will  be  functions  of 
these  parameters.  In  most  cases  the  coefficients  f  will  be  functions  of 
only  one  parameter  (or  they  can  be  made  to  be).  Hence  the  roots  of  f^(S) 
can  readily  be  placed  algebraically  for  the  cases  of  n  =  1,  and  n  =  2. 

For  the  case  of  n  =  3,  a  root  locus  sketch  will  usually  suffice.  Although 
the  coefficients  f  have  been  derived  for  only  up  to  a  fifth  order  case, 
they  could  easily  be  obtained  for  higher  order  cases  if  necessary. 

3-2-2  Applications  of  the  dominancy  technique. 

Example  3-8  (Third  order  characteristic  equation) 

Problem: 

Compensate  the  system  of  figure  (3-11)  with  a  cascade  compensator  to 
obtain: 

1.  Characteristic  roots  at  zeta  =  .5  and  omega  =  40. 

2.  K  -  250. 
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3.  The  specified  roots  are  to  be  made  dominant. 

Solution: 

The  characteristic  equation  of  figure  (3-11)  is: 

S3  +  (4  +  P)S2  +  (4P  +  K  *  )S  +  KP  =  0 
3  2 

or  S  +  (4  +o4  )S  +  (4o£  +  K  @  )S  +  K<^-  =  0 

where  alpha  =  P  and  beta  = 

Now  G  =  K/e(S)  =  K/(S2  +  4S) ,  so  e  =  0,  e,  =  4,  e.  -  1,  U  =  1,  U  =  0. 

o  1  A  2  J 

Equations  (3-42)  are  now  applied  to  obtain: 

Bj  =  -K  +  w2  =  -K  +  1600 

Cl-° 

Dj  =  4w2  -  w3U2  =  -5.7xl04 
$2  =  -4w  +  w2U2  =  1440 
C2  =  -wK  =  -40K 
D2  =  4w2U2  -  w3U3  =  6400 

From  equation  (3-43)  are  obtained: 

0^  «  5.76xl04 

-K+1600 

8  =  (1440) (-5. 76xlQ4)  -  (-K  +  1600)(6400)  (3-74) 

-40K(-K  +  1600) 

Now  for  any  value  of  K,  equations  (3-74)  will  provide  a  value  of  alpha 
and  beta  to  provide  characteristic  roots  at  zeta  =  .5  and  omega  =  40. 

The  value  of  K  will  now  be  chosen  on  the  basis  of  dominance  and  steady 
state  error  considerations.  To  satisfy  the  error  specifications  it  is 
necessary  that  K  2  1000.  Since  f^(S)  is  of  order  one,  the  appropriate 
equations  are  (3-69),  hence: 
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1 


fo  =  «„/»*  -  (aj  -  wj)/2  ?  -  »2  -  2  ?  ^ 

From  Che  characteristic  equation  it  is  seen  that 
a2  =  4  +c>C 

a  =  4  oC  +  K  (3 

a  =  K  cL 
o 

The  real  part  of  the  specified  roots  is  J ^w^  =  20.  Arbitrarily  choos¬ 
ing  a  dominance  factor  of  five,  the  dominance  criteria  becomes: 
fQ  >  5  ^  ^w^  =  100.  To  satisfy  this  requirement,  the  simplest  form  of 

f  will  be  chosen,  namely  f  =  a  /w, .  Therefore  it  is  seen  that  f  = 
o  o  o  1  o 

K  oL  /1600  =  (5.76x10^) _  =  36K  >  100,  where  equation  (3-74) 

1600 (-K  +  1600)  (1600  -  K) 

was  employed. 

The  above  inequality  implies  that  K  >  1180.  Since  K  >1180  also 
satisfies  the  error  specification,  a  value  of  K  =  1200  is  chosen  arbi¬ 
trarily.  Using  this  value  of  k,  the  following  quantities  are  computed: 
alpha  =  144,  beta  =  4.2,  and  fQ  =  108.  As  a  check,  the  expression  fQ  =  a^ 

-2  J  jW^  can  be  employed.  Hence: 

f  =  (4  +  144)  -  2 (.5) (40)  =  108. 

Example  3-9  (Fourth  order  characteristic  equation) 

Problem: 

Compensate  the  system  shown  in  figure  (3-12)  employing  tachometer 
plus  acceleration  feedback  to  obtain: 

1.  Characteristic  roots  at  zeta  =  .5  and  omega  =  2. 

2.  K  -  12. 

e 

3.  The  specified  roots  should  be  dominant. 
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Figure  3-12 


Solution: 


The  characteristic  equation  from  figure  (3-12)  is: 

S4  +  16. 5S3  +  (73  +oC  )S2  +  (82.5  +  (3  )S+25  +  K  -  0 

Here  alpha  =  KK.  and  beta  =  KK  . 

a  t 

G  =  K/e(S)  =  K/(S4  +  16. 5S3  +  73S2  +  82. 5S  +  25) 

By  inspection  it  is  seen  that: 

eQ  =  25,  =  82.5,  e 2  =  73,  e3  =  16.5,  e^  ■»  1,  U2  **  1,  U3  =  0, 

and  =  -1.  By  employing  equations  (3-11)  one  can  find: 
alpha  =  (185  +  K)/4 

and 

beta  =  (K  -  23)/2 

From  the  characteristic  equation  it  is  seen  that: 

a^  =  1,  a3  =  16.5,  a^  =  73  +  c/-  ,  =  82.5  +  (S  ,  aQ  =  25  +  K. 

2 

Since  the  quotient  equation  f^(S)  is  a  quadratic,  i.e.,  S  +  f^S  + 
fQ  =  0,  it  would  be  desirable  from  the  dominancy  standpoint  if  f^  >  5  J  ^w^ 
=  5.  However,  looking  at  the  dominancy  equations  for  this  case  (equations 
(3-71)),  it  is  seen  that  one  of  the  several  expressions  for  f^  is  f^  =  a3  - 
2  f  1W1*  Since  all  the  expressions  for  f^  have  to  be  simultaneously 
satisfied,  it  is  seen  in  this  problem  that  f^  is  a  fixed  constant  since 
a3  and  J  ^w^  are  fixed.  This  is  true  even  though  the  remaining  express¬ 
ions  for  f^  involve  one  or  more  variable  coefficients.  Therefore  f^  is 
found  to  be  14.5.  Observing  the  pertinant  expressions  for  f  it  is  seen 

that  none  of  them  contain  only  constants,  so  the  simpler  expression  fQ  = 

2 

a  /w.  is  chosen.  Now  since  f,  =  14.5  7  5,  a  dominant  situation  already 
o  1  1 

exists,  but  the  system's  performance  can  be  further  improved  by  choosing 
a  reasonable  value  of  zeta  and  omega  for  the  secondary  roots.  From  the 
error  specification  it  is  necessary  that  K/25  -  12,  or  K  S300.  Now 


48 


f  (S)  =  s2  +  14. 5S  I-  aQ/w2  or  f][(S)  =  S2  +  14. 5S  +  6.25  +  .25K.  For 

2 

K  =  300,  f  (S)  becomes:  S  +  14. 5S  +  81.25.  Therefore,  2  ?  »w  = 
i  '  2  2 

2 

14.5,  w2  =  81.25  or  w2  =  9.  Then  zeta  =  .806.  These  are  reasonable 

values  for  2  and  w2  since  the  secondary  roots  taken  by  themselves  would 

produce  much  less  overshoot  and  a  much  smaller  settling  time  than  the 

primary  roots.  Using  this  smaller  value  of  K  one  can  compute  alpha  and  beta. 

alpha  =  121.2  beta  =  138.5 

Since  alpha  =  KK  +  300K  and  beta  =  KK  +  300K  then  K  =  .405  and 
t  t  a  a  t 

K  =  .462. 
a 

As  an  added  bonus  of  the  method,  all  roots  of  the  characteristic  equa¬ 
tion  are  now  known  and  the  time  response  could  be  computed  if  desired. 

Example  3-10  (Fifth  order  characteristic  equation) 

Problem: 

Compensate  the  system  shown  in  figure  (3-13)  to  obtain: 

1.  Characteristic  roots  at  zeta  *  .5,  and  omega  =  4. 

2.  K  i  8. 

e 

3.  The  specified  roots  should  be  made  dominant. 

Solution: 

The  characteristic  equation  is  after  letting  alpha  =  P  and  beta  =  X 

S5  +  (  +  17)S4  +  (17o4+  84)S3  +  (148  +  KKfc  +  84c4)S2  +  (148c/  + 

KK  oi  +  K  (2  +  80)S  +  ©(  (K  +  80)  =  0. 

G  =  K/e(S)  =  K/(S4  +  17S3  +  84S2  +  148S  +  80) 

From  the  above  it  follows  that:  e  =  80,  e.  =  148,  e  =  84,  e,  =  17, 

o  ’1  2  3 

e^  =  1,  U2  =  1,  =  0,  =  -1,  and  =  -1.  Employing  equations  (3-42) 

one  obtains: 
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=  K  +  1196 


Dx  =  -1986 

Therefore: 


B2  =  492 
C2  =  -4K 
D2  =  -1276 


alpha  =  1986/(K  +  1196)  beta  «  (1276K  -  555000) /[4K(K  +  1196)] 

For  beta  to  be  positive  it  is  necessary  that  K  be  greater  than  435.  The 
accuracy  specification  implies  that  K  be  greater  than  640.  For  K  =  640, 
beta  =  .0555  which  is  out  of  the  desired  range  of  . 1  —  (S  -  10. 

Of  course,  increasing  K  makes  beta  even  smaller.  At  this  point  one  could 
design  a  multiple  section  compensator  in  accordance  with  section  (4-6-2) 
or  feedback  compensation  could  be  added.  The  latter  course  of  action 
is  chosen  and  tachometer  feedback  not  enclosing  the  cascade  compensator 
is  chosen.  The  applicable  equations  are  (3-54).  From  them  are  obtained 
the  following,  where  the  above  alpha  and  beta  now  become  o^.'  and  ^  '. 


alpha 


1 _  (1986  -  16KK  ) 

K  +  1196  C 


beta  = 


-492(4)Kt  -  1986Kt  +  16KK*  +  4k  +  ^  , 
K  +  1196  t 


To  simplify  the  analysis  let  K  =  640  which  meets  the  accuracy  specifica¬ 
tion.  Then  alpha  and  beta  become 


.  .  1986  -  10210K 

alpha  =  _ t 


1836 


Here  alpha  is  positive  for  less  than  .193. 

beta  =  5.64K?  +  1.83K  +  .0555 

t  t 


(3-75) 


(3-76) 


The  appropriate  dominancy  equations  are  (3-73)  where  it  is  seen  that 

neither  f  ,  f, ,  or  f „  are  restricted  to  constant  values.  From  the  charact- 
o*l'  2 

eristic  equation  it  is  seen  that: 
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4 
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1 


a.  88 +17 

4 

a3  =  17d  +  84 

a2  =  148  +  KKt  +  84oC 

a =  148°4  +  KKt<4+  K  (3  +  80 

a  =  oL  (K  +  80) 
o 

On  Che  basis  of  these  equations  the  following  expressions  are  chosen 
from  equations  (3-73)  for  the  f^(S)  coefficients: 
f2  =  -  2  ^  ^w^  =  +  13 

f  =  a0  -  2  J  .w,a.  +  U„w2  =  13  oc  +  16 
1  3  /  1  1  4  3 

f  -  a  /w?  -  180  JL 
o  o  1 

Then  f^(S)  becomes: 

S3  +  (o4  +  13)S2  +  (13^  +  16)S  +  180o4  =  0  (3-77) 

Equation  (3-77)  contains  only  one  variable  so  it  can  be  put  in  root  locus 
form  as  follows: 

S3  +  13S2  +  16S  +°£(S2  +  13S  +  180)  =  0 
3  2 

After  dividing  by  S  +  13S  +  16S  the  above  equation  becomes: 

(S2  +  13S  +  180)  =  (S  +  6.5  +  ill. 7)  =  -1  (3-78) 

S(S2  +  13S  +  16)  S(S  +  1.37)(S  +  11.6) 

By  plotting  the  root  locus  poles  and  zeros  of  equation  (3-78)  or 
by  inspection  if  one  wishes,  it  is  apparent  that  there  will  be  a  real 
root  of  f^(S).  Also  the  magnitude  of  the  real  part  of  the  complex  roots 
can  never  be  greater  than  6.5,  and  it  approaches  this  value  as  alpha  tends 
to  plus  infinity.  For  dominancy  therefore,  a  large  value  of  alpha  is  de¬ 
sired.  From  equation  (3-75),  the  maximum  value  of  alpha  is  for  Kt  =  0, 
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but  then  beta  is  too  small  since  Kt  =  0  implies  the  case  of  cascade  com¬ 
pensation  only.  It  is  obvious  that  should  be  as  small  as  possible, 
which  from  equation  (3-76)  implies  that  beta  should  be  as  small  as  possi¬ 
ble.  Applying  the  lower  limit  of  beta  =  .1  to  equation  (3-76)  one  can 
show  after  solving  the  quadratic  in  Kt  that  Kt  =  .023  is  the  necessary 
value. 

From  equations  (3-75)  and  (3-76)  one  can  solve  for  alpha  and  beta 
obtaining: 

alpha  =  .955  =  P  beta  =  .1  *  X 

The  cascade  compensator  is  then: 

.1(S  +  9.55) 

(S  +  .955) 

This  can  be  realized  by  an  R-C  lag  network.  Equation  (3-77)  for  f^(S) 
becomes : 

S3  +  13.955S2  +  28. 4S  +  172  =  0  (3-79) 

A  root  locus  of  equation  (3-77)  for  0  S  o4  Soo  is  shown  in  figure  (3- 
14).  The  root  locus  was  done  by  the  computer  program  presented  in  sec¬ 
tion  (5-2).  The  roots  of  equation  (3-79)  were  obtained  from  the  printed 
out  data  from  the  program  and  are  as  follows:  for  alpha  =  .955;  S=  - 
.585  +  j3.64  and  -12.8.  Since  .585  is  less  than  which  is  two,  the 

primary  or  specified  roots  are  dominant.  Now  the  maximum  value  of  alpha 
is  obtainable  when  K  =  =  0  and  even  for  this  extreme  case,  the  speci¬ 

fied  roots  will  not  be  dominant.  To  solve  the  problem,  therefore,  a 
compromise  on  the  specifications  must  be  made.  One  approach  would  be  to 
specify  omega  as  some  value  less  than  four  and  repeat  all  the  preceding 
calculations,  but  one  could  also  do  the  following: 

From  equation  (3-75)  observe  that  alpha  can  be  decreased  by  in¬ 
creasing  Kfc.  From  the  root  locus  data  it  is  observed  that  if  alpha  = 
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,1.17,  then  -  .5,  and  w^  =  1.34.  If  this  reduced  value  of  omega 
can  be  tolerated  then  the  problem  is  solved  since  the  secondary  roots 
located  at  ^  ^  and  are  now  dominant  and  only  the  settling  time  has  been 
modified  (due  to  the  decreased  value  of  omega).  On  the  other  hand  since 
the  roots  ^ ^  and  w^  are  still  located  at  zeta  =  .5  and  omega  =  4,  the 
overall  system  response  will  probably  be  acceptable  to  the  designer. 

The  final  results  are  therefore  as  follows:  . 

=  .173  obtained  from  equation  (3-75) 
alpha  =  P  =  . 117 

beta  =  gamma  =  .5405  obtained  from  equation  (3-76) 

Characteristic  roots  are  located  at  zeta  =  .5,  omega  =  4;  zeta  =  .5,  omega 
*  1.34;  and  a  real  root  at  S  =  -11.8.  The  cascade  compensator  is  of  the 
form: 

. 5405 (S  +  .216) 

(S  +  .117) 

This  can  be  synthesized  as  an  R-C  lag  network. 
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3-3  Some  sketching  techniques.  ^ 

The  graphical  solution  will  be  discussed  in  following  sections  and 
of  course  this  involves  plotting  curves.  Due  to  the  complexity  of  the 
parameter  plane  equations,  two  methods  of  curve  plotting  are  recommended. 
The  easier  and  faster  method  is  by  means  of  a  digital  computer.  A  com¬ 
puter  program  which  will  do  this  is  presented  in  section  (6).  The  second 
method  is  by  means  of  a  desk  calculator  or  slide  rule.  To  facilitate  the 
use  of  the  latter  method,  the  parameter  plane  equations  are  manipulated  to 
put  them  into  forms  more  suitable  to  non-computer  plotting. 

Due  to  the  time  limitations  usually  imposed  on  digital  computer 
operation  and  due  to  the  time  consumed  in  plotting  the  curves  by  hand,  it 
would  be  desirable  to  just  sketch,  say  the  zeta  equals  zero  and  the  zeta 
equals  one-half  curve  for  various  values  of  omega  to  see  if  the  type  of 
compensation  proposed  will  do  the  job.  This  type  of  rapid  sketching  is 
also  helpful  when  choosing  a  scale  for  the  digital  computer  graph  plot. 

For  this  reason,  in  the  following  sections,  special  characteristics  of  the 
zeta  equals  zero  and  the  zeta  equals  one-half  curve  have  been  derived  for 
characteristic  equations  of  order  two  through  five. 

3-3-1  Table  of  symbols. 

In  table  (3-2)  are  listed  various  symbols  which  apply  to  the  deriva¬ 
tions  made  in  section  (3-3-2).  The  change  of  variables  as  indicated  in 
table  (3-2)  is  useful  because  many  of  these  products  and  sums  appear  re¬ 
peatedly.  When  doing  hand  computations  of  the  parameter  plane  curves  it 
is  only  necessary  to  compute  these  quantities  once,  so  substantial  labor 
is  saved. 
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TABLE  (3-2) 


J 

K 

L 

M 

N 

P 

R 

T 

U 

A 

B 

C 

D 

E 

F 

G 

H 

I 

X 


-  Vi  +  bjc2 


b  c, 
o  1 

Cld2 

C  di 

o  1 

Vo 

Vi 


-  b.c 

1  o 

-  c2dl  '  cc 

-  C  d 

1  o 

-  b  d, 

o  1 

-  b  d  +  b 

1  2  c 


y 

z 


do°2  ‘  d2co 
b2co  -  boc2 


bod2 

b2c3 


b2do 

b3c2 


b3co  ‘  boc3 
c2d3  •  c3d2  '  C1 
C3do  •  cod3 
b3d2  -  b2d3  +  bl 


b  d. 
o  3 


M 

3  o 


B  +  J 
D  +  L  +  c 


P  -  b  +  F 
o 

c-d.  -  c.d.  +  c 
3  1  13  o 


b3cl 

bld3 


9( -  c4d2 

■£-  -  c2a4 

9li  =  e2d3 

^  =  cld4 

e  =  a0c4 

A  -  dlc3 

T  *  cod3 
t>  -  b4c3 

e  =  b2c4 

V  -  b3c2 

6  m  cob4 

S  -  b2d4 

A  s  b  d, 
o  4 

P  ■  b3d4 

+  =  d3«4 

v  •  b4di 


blc3 

b3dl 

c2d4 

d2c4 

c3d2 

c4dl 

c  d. 
o  4 

Cld3 

C3do 

b3C4 

c0b7 
2  4 

b2c3 

boc4 

b4d2 

dob4 

d3b4 

c3d4 

bXd4 
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3-3-2  Basic  derivations. 

Case  J.  (Second  order  characteristic  equation) 

The  characteristic  equation  is  of  the  form: 


S  +  c1^+  d^S  +  bQ0<.+  cq^+  dQ  =  0 

Employing  equations  (2-10)  one  can  obtain: 


B.  =  -b 
1  o 


C.  =  -c 
1  o 


B0  =  -b  w 
2  o 


°2  ‘  -cl" 


D.  =  -d  +  V2  D,  -  -d.  +  w2U, 

1  o  2  12 

Solving  for  alpha  and  beta  using  equations  (2-11)  results  in: 

alpha  =  (c^w2  -  cqU^w  +  M)/K  (3-80) 

beta  =  (-bjW2  •(-  b  U2w  4-  N)/K 

To  find  the  maximum  and  minimum  points,  the  first  derivatives  of  alpha 
and  beta  with  respect  to  omega  are  taken  and  set  equal  to  zero. 


d  ot-  /dw  =  (2c.w  -  c  U„)/K  =  0  or  w  =  c  U-/2c.  (3-81) 

1  O  Z  O  Z  1 

d  |3  /dw  =  (-2bjW  +  bQU2)/K  =  0  or  w  =  bQU2/2b^ 

In  equations  (3-80)  after  letting  w  =  0  one  obtains: 

oL  =  m/K  P  =  N/K  (3-82) 

Letting  w  tend  to  plus  infinity  one  obtains: 

alpha — >  +o o  beta — ^  (3-82)2 

Equations  (3-80)  through  (3-82)  are  valid  for  all  values  of  zeta 
between  zero  and  one.  Since  U2  =  0  when  zeta  =  0,  then  equations  (3-80) 
and  (3-81)  become: 

alpha  =  (c^w2  +  M)/K  beta  =  (-b^w2  +  N)/K  (3-80a) 

d°^/dw  =  2c^w/K  w  =  0  (3-81a) 

d  ^  /dw  =  -2b^w/K  w  =  0 

3 

The  relative  magnitudes  of  the  coefficients  determine  whether  alpha 
and  beta  approach  plus  or  minus  infinity. 


58 


Equations  (3-82)  remain  unchanged.  Also  since  =  1  when  zeta  =  .5 
one  can  obtain: 

alpha  =  (c^w2  -  cqw  +  M)/K  beta  =  (-b^w2  +  bQW  +  N)/K  (3-80b) 

d  o L  /dw  =  (2CjW  -  co)/K  w  =  cjlc^  (3-81b) 

d  £  /dw  =  (-2b  w  +  b  )/K  w  =  b  /2b 

r  1  o  o  1 

Equations  (3-82)  remain  unchanged. 


Case  II.  (Third  order  characteristic  equation) 
The  characteristic  equation  is  of  the  form: 


S3  +  (b2 o£.+  C2  P'*’  <^2^  +  d^)S  +  b^oi.  +  cc  ^  =  0 

proceeding  as  before  the  following  expressions  are  obtained: 


alpha  = 


~Vc2»3 


+  C2U2> 


-  c,U0w3  +  (L  +  c  +  c  Ujw2  +  U0Rw  +  M 

1  / _ O _ O  J _ / _ 

Jw2  +  U2Tw  +  K 


beta  -  »4<-»2D2  +  b2U3>  +  blU2”3  +  <P  *  bo  *  b„V"2  +  °2U"  +  N 


Jw  +  U2Tw  +  K 


For  w  =  0, 

alpha  =  M/K  beta  =  N/K  (Same  as  equation  (3-82)) 

For  w  +  c6 

alpha  =  +  o2>  beta  =  +  oo  (Same  as  equation  (3-82)) 

Due  to  the  increased  complexity  of  the  expressions  for  alpha  and  beta, 
the  derivatives  are  only  computed  for  the  zeta  equals  zero  and  the  zeta 
equals  one-half  curves. 

Letting  zeta  =  0,  equations  (3-83)  become: 


i  u  cow4  +  Lw2  +  M 
alpha:  2 _ 

Jw2  +  K 


4  2 

,  .  -bw  +  Pw  +  N 

beta  =  _2 _ 

Jw2  +  K 


(3-83a) 


For  w  =  0  and  w  =  °°  ,  equations  (3-82)  apply. 

Taking  derivatives  of  equations  (3-83a)  one  obtains: 
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,2,  ,A  2  w4  coJ  +  w2  2c0K  +  KL  -  JM 
d  <*-  /dw  =  2 _  2 _ 

2  2 
(Jw^  +  K)Z 

The  derivative  goes  to  zero  for: 

w2  =  -2c2K  +  J  (2c2K)2  -  4c2J(KL  -  JM) 

— 

d2  $  /dw2  -  w4<"b2J)  •  «2(2b2K)+KP  -  NJ 

'  -  g 

(JwZ  +  K) 

The  derivative  goes  to  zero  for: 

w2  =  2b2K  +  J(  2b2K)2  +  4b2J(KP  -  NJ) 

2b^J 

Letting  zeta  ■  .5,  equations  (3-83)  become: 

4  3  2 

alpha  «  c2w  -  c^w  +  (I,  +  cq)w  +  Rw  +  M 

Jw2  +  Tw  +  K 

beta  =  -b2w4  +  b^w3  +  (P  -  bQ)w2  +  Uw  +  N 

Jw2  +  Tw  +  K 


(3-84) 


(3-84) 


(3-84) 


(3-84) 


(3-83b) 

(3-83b) 


For  w  =  0  and  w  =  oO  ,  equations  (3-82)  apply.  Taking  derivatives  of 
equations  (3-83b)  one  obtains: 

d  oL  /dw  =  [2c2Jw5  +  (-3clJ  +  4c  T  -  c2T  +  2CjJ)w4  +  (4c2K  - 

2c1T)w3  +  (JR  +  LT  +cqT  -  3cxK  +  2JR)w2  +  (-2JM  + 

2  2 

2LK  +  2cqK)  +  KR  -  MT]/(Jw  +  Tw  +  K)  (3-85) 

d  §  /dw  =  [-2b2Jw3  +  (-2b 2T  +  bjJ)w4  +  (-4b2K  +  2b^T)w3  + 

(3b  ^  +  PT  -  bQT  -  UJ)w2  +  (2PK  -  2bQK  +  2JN)w 
+  (UK  -  NT) ] / ( Jw2  +  Tw  +  K)2 


In  the  general  case,  equations  (3-85)  involve  fifth  order  poly¬ 
nomials  in  omega.  In  specific  problems,  however,  some  quantities  in 


these  equations  will  be  zero,  enabling  one  to  more  readily  solve  for  the 
critical  values  of  omega. 

Case  III.  (Fourth  order  characteristic  equation) 

The  characteristic  equation  is  of  the  form: 

S4  +  (^3°^-  +  03  (3  +  +  (b2c4  +  c2  +  d2)S  + 

(bl  +  ci  ^  +  d1)S  +  bQOl+  cq  £  +  dQ  =  0 

Proceeding  as  before  the  following  expressions  are  obtained: 

alpha  =  [  w6c3(U3  -  U^)  +  w5  c2  (U4  -  U^)  +  w4(-C  -  ^ 

+  ClU3  +  U2  (C  +  c1))  +  w3(U2(X  -  cq)  -  cqU4)  + 

w2(-U-D  +  L  +  c  )  +  wU.R  +  M)/A 
3  o  / 

beta  =  [w^b3(U2U4  -  U2)  +  w^b2^U2U3  -  U^)  +  w  (U2(E  - 

bl)  -  bjlhj  -  U3(E  -  bx))  +  w3U2(Z  +  bQ)  +  boU4) 

+  w2(P  -  bQ  -  U3F)  +  wU2U  +  N]M 

A  =  w4A(U2  -  U3)  +  w3YU2  +  w2  (-BU3  +  J)  +wTU2  +  K 

For  w  =  0  and  w  =  °°  ,  equations  (3-82)  apply.  Letting  zeta  =  0, 

equations  (3-86)  become: 

6  4  2 

,  .  c„w  +  Cw  +  Hw  +  M 

alpha  =  __3 _ _ 

4  2 

Aw  +  Gw  +  K  (3-86a) 

5  4  2 

,  -b0w  +  Ew  +  Iw  +  N 

beta  =  3 _ 

4  2 

Aw  +  Gw  +  K 

For  w  =  0  and  w  =  ,  equations  (3-82)  apply.  Taking  derivatives  of 

equations  (3-86a)  one  obtains: 

d^  /dw^  —  [w  Ac^  +  w  3c^G  +  w  (AH  +  CG  +  3Kc^)  +  w  (2CK  + 

4  2  2 

GH  -  2A(H  +  M))  +  KH  -  G(H  +M)]/(Aw  +  Gw  +  K) 

(3-87) 
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d2  /3  /dw2  =  [-w3Ab3  -  wb2b3G  +  w4(EG  -  AI  +  3Kb3>  +  w2  (GI  + 

2EK  -  2AN)  +  KI  -  NG]/(Aw4  +  Gw2  +  K)2 

In  the  general  case,  equations  (3-87)  involve  eighth  order  poly¬ 
nomials  in  omega.  In  specific  problems,  however,  some  quantities  in 
these  equations  will  be  zero,  enabling  one  to  more  readily  solve  for  the 
critical  values  of  omega. 

Letting  zeta  «  .5,  equations  (3-86)  then  become: 


alpha  **  C3W  "  C2W^  +  (C  +  c^)w4  +  Xw3  +  (L  +  c^)w2  +  Rw  +  M 

4  3  2 

Aw  +  Yw  +  Jw  +  Tw  +  K 


(3-86b) 

.  .  -b.w^  +  b„w3  +  (E  -  b. )w4  +  Zw3  +  (P  -  b)w2  +  Uw  +  N 

beta  *  _ J  A  l _ 

Aw4  +  Yw3  +  Jw2  +  Tw  +  K 

For  w  =  0  and  w  -o°  ,  equations  (3-82)  apply.  Due  to  the  increased 
complexity  of  equations  (3-86b),  it  does  not  appear  practical  to  compute 
the  derivatives  for  the  general  case. 

Case  IV.  (Fifth  order  characteristic  equation). 

The  characteristic  equation  is  of  the  form: 

S5  +  (b^  +  c ^  (3  +  d4)S4  +  (bjO*-  +  c3  (3  +  +  c2  ^ 

+  d2)S2  +  (b^ot  +  c^  {3  +  d^JS  +  b'OL  +  cq  j3  +  dQ  =  0 

Proceeding  as  before  the  following  expressions  are  obtained: 

alpha  =  tw8c4(U5  +  U2)  -  w?c3(U2U5  +  U^)  +  w6(  &  (U^ 

-u2)  +  c2(u5  +  u2u4))  +  w5(u2u3  X  +  u4£-  Clu4) 

+  w4(U2>-  U31«+  U3^J  +  cq)  +  w3(U4^+  U )  + 
w2(U3T+  L  +  cq)  +  wU2R  +  M]/A  (3-88) 
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beta  =  t-w8b4(U l  -  U3U5)  +  vApu^  +  U^)  +  «6f(U2U4 

-U2  "  b2(U2U4  +  U5»  +  w5<(U2U3  "  +  blV  + 

w4(-U3(E  -  bp  +  U3T  +  U2(E  -  bp  +  boU5)  +  w3 
(U4/<+  U2(Z  +  bQ))  +  w2(-Up  +  p  -  bQ)  +  wU2U  +  N2/A 


/  2 

A  =  w60(U2U4  -  U3)  +  W5  0  (U4  -  U£U3)  +  w  ^  (U3  -  up  + 
w3(U4€  +  Up)  +  w2(J  -  U3B)  +  wUp  +  K 

For  w  =  0  and  w  =  °°  ,  equations  (3-82)  apply.  Letting  zeta  =  0,  equa¬ 

tions  (3-88)  become: 

alpha  -  V\  *  "6(C2  ill  1 1,4 ^  ''ft  +  *  "2<L  *  ~T)  +  M 

A  1  (3-88a) 

beta  =  w8b4  -  w6  p  (b2  +  1)  +  w4(E  -  bt  -  P  +  bQ)  w2(F  +  P-bQ  +  N 

— 

A  ^  =  -w80  -  w4X  +  w2(J  +  B)  +  K 
Letting  zeta  =  .5,  equations  (3-88)  become: 

alpha  -  [w7c3  -  w6('f  +  2cp  +  w5(c1  -  £)  +  w4<)T|  +  w3($  - 

0  )  +  w2(L  +  cq)  +  wR  +  M]/  &  2  (3-88b) 

beta  =  [-w8b4  -  w7b3  +  w8  P  (2b2  -  1)  +  w3  (S-bp  +w 

(E  -  b.  -  b  )  +  w3(Z  +  b  -/()  +  w2(p  -  b  )  +  wU  +  N]/  A  2 

a  6rt  -  w50  -  w4  +  w3(Y  -  6  )  +  w2J  +  wT  +  K 

A  2  -  -w  0 

Example  3-11  (Case  I  example) 

Problem: 

Sketch  the  zeta  =  0  and  the  zeta  **  .5  curves  for  the  following 
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characteristic  equations.  Let  the  abscissa  variable  be  alpha  and  the 
ordinate  variable  be  beta. 

S2  +  (  oC  +  2)S  +  (3+1  =  0 

Solution: 

From  equations  (3-80)  it  is  found  that: 

2 

alpha  =  -  2  beta  =  w  -  1 

From  equations  (3-81)  : 

d  /dw  =  1*2  d  jS  /dw  =  2w 

From  equations  (3-82): 

alpha  =  -2  and  beta  =  -1  for  omega  =  0. 

Therefore  for  zeta  =  0  the  following  relations  apply: 

2 

alpha  =-2  beta  =  w  -  1  and  d  (3  /dw  =  2w, 

which  implies  a  minimum  for  beta  at  omega  =  0.  The  zeta  =  0  curve  is 
therefore  a  vertical  line  in  the  alpha-beta  plane,  going  from  alpha  =  -2, 
beta  °  -1  to  plus  infinity  as  omega  increases  from  zero  to  infinity. 
Setting  zeta  to  .5,  one  finds  that: 

2 

alpha  =  w  -  2  beta  =  w  -  1 

or 

w  =  cL  +  2  beta  =  (o£+  2)2  -  1. 

The  alpha-beta  curve  is  therefore  a  parabola  with  vertex  at  alpha  =  -2 
and  beta  =  0;  symetric  about  the  line  alpha  =  -2,  and  opening  in  the  plus 
beta  direction.  The  curve  can  readily  be  plotted  after  calculating  a 
few  critical  points. 

The  above  curves  are  plotted  in  figure  (3-15). 

Example  3-12  (Case  II  example) 

Problem: 

Repeat  example  (3-11)  for  the  following  equation: 

S3  +  (  cL+  1)S2  +  ((3+2)S  +  3  =  0 
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Solution: 


From  equations  (3-83a)  for  zeta  =  0,  it  is  found  that: 

2  2 
alpha  =  - 1  +  3/w  beta  =  w  -  2 

For  omega  =  0,  alpha  =  plus  infinity  and  beta  =  -2. 

From  equations  (3-84): 

d^  o4  /dw^  =  -3/w^  d^  ^  /dw ^  =  1 

But  d  ^  /dw  =  2w,  so  alpha  has  a  minimum  at  omega  equals  infinity.  Beta 

has  a  minimum  at  omega  =  0  and  it  is  -2.  The  curve  is  asymptotic  to  the 

lines  alpha  =  -1  and  beta  =  -2. 

The  following  points  are  readily  obtained: 

omega  alpha  beta 

1.414  .5  0 

1.732  0  1 

From  equations  (3 -83b)  for  zeta  =  .5,  it  is  found  that: 

2  2 

alpha  =  w  -  1  +  3/w  beta  =  w  -2  +  3/w 

For  omega  =  0,  alpha  =  +  o°  ,  and  beta  =  +oO 
Also, 

d  oC  /dw  =  1  -  6/w^  d  /dw  =  2w  -  3/w^ 

It  is  therefore  obvious  that: 

3  3 

d  /dc4  =  w(2w  -  3)/(w  -  6),  from  which  the  critical  points  are 

1/3 

found  to  be:  w  =  0  and  w  =  (3/2)  =  1.145.  This  point  is  a  minimum 

since  the  second  derivative  is  negative. 

The  following  values  can  be  obtained: 


omega 

alpha 

beta 

1.145 

2.435 

1.93 

.1 

300 

28 

1 

3 

2 

3 

2.33 

8 

The  curves  are  sketched  in  figure  (3-16). 


Example  3-13  (Case  III  example) 


Problem: 


Repeat  example  (3-11)  for  the  following  equation: 

S  +  (  o£-  +  1)S  +  2S  +  oi.  S  +  jS  =  0 

Solution: 

From  equations  (3-86a)  for  zeta=0: 

alpha  =  -w2/(w2  -  1)  beta  =  (-w6  +  3w4  -  2s2) /(w2  -  1) 

As  omega  tends  to  plus  infinity,  both  alpha  and  beta  tend  to  minus 
infinity.  As  omega  tends  to  zero,  both  alpha  and  beta  tend  to  zero. 

Also: 

?  9  9  9  2  2 

d  e*-  /dw  =  [  - (w  -  1)  +  w  ]/(w  -  1) 

d^  ^  /dw2  =  [(w2  -  l)(-3w4  +  6w2  -2)  -  (-w^  +  3w4  -  2w2/|/(w2  -  l)2 

It  can  be  seen  from  the  above  derivatives  that  d  ^  /do1-  =  0,  implies  six 

critical  points.  It  appears  unpractical  to  compute  them. 

From  equations  (3-86b)  for  zeta  =  .5,  it  can  be  shown  that: 

3  2 

alpha  =  -w  +  2w  =  w(2  -  w  ) 

6  4  3  2  2  4  2 

beta  =  w  -  2w  -  w  +  2w  =  w  (w  -  2w  -  w  +  2) 

As  omega  tends  to  plus  infinity,  alpha  tends  to  minus  infinity  and  beta 

tends  to  plus  infinity.  As  omega  tends  to  zero,  alpha  and  beta  both  tend 

to  zero. 

Using  a  slide  rule  or  desk  calculator  along  with  the  above  informa¬ 
tion  the  curves  can  be  plotted  as  shown  in  figure  (3-17). 

Example  3-14  (Case  IV  example) 

Problem: 

Repeat  example  (3-11)  for  the  following  equation: 

S5  +  3S4  +  S3  +  2S2  +  c/  S  +j9  =  0 

Solution:  • 

From  equations  (3-88a)  for  zeta=0: 

2  4  2  4 

alpha  =  w  -  w  beta  =  2w  -  3w 
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exists  at  omega  =  .578.  The  following  points  can  readily  be  calculated: 


omega 

alpha 

beta 

.578 

.222 

.333 

.816 

.222 

0 

1  0  -1 

From  equations  (3-88b)  it  is  seen  that: 

3  5  3  2 

alpha  =  -3w  +  2w  beta  =  w  -  w  +  2w 

Therefore: 

d  /dw  =  -9w^  +  2  d  /dw  »  5w^  -  3w^  +  4w 

Hence  d  ft  /doL  -  0,  implies  that  w(5w  -  3w  +  4)  =  0.  It  is  not 
worthwhile  to  factor  the  above  cubic  in  omega  to  obtain  the  three  corres¬ 
ponding  critical  points.  The  critical  point  at  omega  equal  zero  indicates 
that  the  curve  starts  at  the  origin.  A  slide  rule  or  desk  calculator  can 
now  be  employed  along  with  the  above  information  to  plot  the  curves  as 
shown  in  figure  (3-18). 
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Example  3-15 


Problem: 

For  the  system  shown  in  figure  (3-19)  let  alpha  =  and  beta 

=  K^K^,  and  satisfy  the  following  requirements: 

1.  Assuming  that  a  computer  will  plot  the  curves  for  zeta  = 

0,  .5,  and  1,  determine  an  appropriate  graph  scale. 

2.  Only  first  quadrant  values  of  alpha  and  beta  are  of 
interest  and  the  graph  is  eight  inches  wide  and  fourteen 
inches  high.  Alpha  is  the  abscissa  and  beta  is  the 
ordinate. 

3.  Due  to  bandwidth  considerations,  omega  should  be  less  than 
1500. 

Solution: 

From  figure  (3-19)  the  characteristic  equation  is  determined  to  be: 

s3  +  1100S2  +  (lo5  +  o4  )s  +  £  =0 
From  this  it  is  seen  from  table  (3-2)  that: 


L  = 

1 

P  =-1100 

J  =  0 

M  = 

io5 

N  =  0 

K  =  -1 

R  = 

-1100 

T  =  0 

U  =  0 

Using  equations  (3-83a) : 

alpha  =  (Lw2  +  M)/(Jw2  +  K)  =  w2  -  105 
beta  =  (Pw2  +  N)/(Jw2  +  K)  =  1100w2 
For  omega  =  1500: 

alpha  =  21. 5x  105  beta  =  24. 7x  108 

From  equations  (3-84)  or  by  taking  derivatives  directly: 

d2  o*-  /dw2  =  1  d2  (9  /w2  =  1100 

Hence: 

d  £  /doL  =  1100 

Letting  zeta  =  .5  and  employing  equations  (3-83b) : 
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Figure  3-3-9 


alpha  =  (Rw  +  103) / (Tw  -  1) 

beta  =  (w3  -  1100w2  +  Uw)/(Tw  -  1) 

Therefore: 

alpha  =  llOOw  -  10"*  beta  =  -w2(w  -  1100) 

It  can  be  seen  that  the  zeta  =  .5  curve  is  not  in  the  first  quadrant  for 
omega  greater  than  550,  so  this  curve  does  not  influence  the  scaling  pro 
blem. 

The  origin  point  of  the  curves  for  omega  =  0,  is  obtained  from  equa 
tions  (3-82)  resulting  in: 

alpha  =  -10"*  beta  =  0 

From  the  above  data  it  is  concluded  that  it  would  be  best  to  scale 
using  the  values  of  alpha  and  beta  for  the  zeta  =  0  curve  corresponding 
to  w  =  1500. 

Therefore: 

alpha-scale  =  3xl03  units  per  inch 

g 

beta-scale  =  2x10  units  per  inch 

Example  3-16 

Problem: 

For  the  system  shown  in  figure  (3-20)  it  is  desired  to  have  roots 
with  zeta  greater  than  .5  and  a  maximum  error  coefficient.  Plot  para¬ 
meter  plane  curves  to  find  the  best  values  for  K  and  Z. 

Solution: 

The  characteristic  equation  is: 

s4  +  15S3  +  150s2  +  (100  +  K)S  +  KZ  =  0 
Let  alpha  =  K  and  beta  =  KZ  to  obtain: 

S4  +  15S3  +  150S2  +  (100  +  c4  )  S  +  (S  =0 
From  equations  (3-86a)  corresponding  to  zeta  =  0: 

alpha  =  (Cw4  +  Hw2  +  M)/(Aw4  +  Gw2  +  K) 

4  2  4  2 

beta  =  (Ew  +  Iw  +  N)/(Aw  +  Gw  +  K) 
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Using  table  (3-2)  one  can  obtain: 


C  =  0 
H  =  -15 
M  =  1000 
A  -  0 
G  =  0 


K  =  -1 
E  =  1 
I  =  -150 
N  »  0 
X  =  1 


R  =-150 
Z  =  15 
P  =  -150 
U  =  0 


Then: 

alpha  =  15w2  -  1000  beta  =  -w2(w2  -  150) 

Hence : 

2 

Beta  is  greater  than  zero  for  w  less  than  150  and  alpha  is  greater 
2 

than  zero  for  w  greater  than  66.5,  so  both  alpha  and  beta  are  greater 

2 

than  zero  for  66.5  <  w  <  150. 

Letting  zeta  »  .5  and  employing  equations  (3-86b)  it  is  found  that: 
alpha  =  -w(w2  -  15)  -  1000  beta  =  -15w2(w  -  10) 

From  these  it  is  seen  that  for  omega  greater  than  ten,  beta  is  less  than 
zero  and  alpha  is  less  than  zero.  For  omega  greater  than  zero  and  less 
than  ten,  beta  is  greater  than  zero  and  alpha  is  less  than  zero.  Hence 
the  zeta  =  .5  curve  is  not  in  the  first  quadrant.  Since  the  order  of 
the  characteristic  equation  is  only  four,  it  is  safe  to  assume  from  the 
above  data  that  all  constant  zeta  curves  for  zeta  greater  than  or  equal 
to  .5  are  not  in  the  first  quadrant.  Since  the  specifications  of  the 
problem  cannot  therefore  be  met,  it  is  unnecessary  to  plot  the  curves. 
Example  3-17 

Problem: 

Referring  to  figure  (3-21), 

1.  Determine  a  scale  for  alpha  and  beta  and  plot  the  curves. 

2.  Place  a  pair  of  roots  with  zeta  :  .5. 

3.  Make  the  roots  dominant. 

4.  Due  to  bandwidth  considerations,  omega  must  not  be  greater 
than  4000. 
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Figure  3-21 


Solution: 


From  figure  (3-21)  the  characteristic  equation  is  found  to  be: 

S4  +  1005S3  +  (5100  +  cL  )S2  +  (105  +  5  0^  )S  +  100  +  /S  =  0 

Here  alpha  =  and  beta  = 


Since  the  purpose  of  plotting  the  curves  is  to  determine  alpha  and 
beta  to  meet  specifications  2,3,  and  4,  the  sketching  techniques  are 
first  employed  to  determine  if  curve  plotting  is  necessary. 

From  table  (3-2)  one  can  find: 

c  =  0  L  =  -1  R  =  -5100 


o 

X  =  1 
Z  =  4925 
A  =  0 
T  =  1 


M  = 

io5 

E  = 

-1000 

U  - 

5. lxlO5 

N  - 

1 

I—* 

o 

Y  = 

0 

J  = 

0 

K  = 

-5 

P  = 

74600 

From  the  above  values  along  with  equations  (3-86b)  one  finds: 

alpha  =  [w2  (w  -  5100)  +  105]/(w  -  5) 

It  can  be  seen  that  for  omega  greater  than  5100,  alpha  is  positive, 
and  for  omega  less  than  5100,  alpha  is  negative. 

beta  [w4  (w  -  1005)  +  4925w3  +  74500w2  +  10  7  (.051w  -  1)]/ 
(w-5) 

Hence  for  omega  greater  than  or  equal  to  5100,  beta  is  positive. 

To  avoidKpositive  feedback,  it  is  concluded  that  zeta  =  .5  can  be 
obtained  only  for  values  of  omega  greater  than  5100,  which  violates 
specification  4.  It  is  therefore  not  necessary  to  plot  the  curves. 
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3-4  Graphical  solutions  on  the  parameter  plane. 

3-4-1  Advantages  of  the  graphical  solution. 

The  previously  discussed  algebraic  solutions  have  the  disadvantage 
that  a  fixed  value  of  zeta  and  omega  must  first  be  chosen.  In  some  cases, 
the  remainder  polynomial  can  then  be  modified  to  ensure  that  the  specified 
roots  are  dominant.  However,  it  is  not  always  possible  to  guarantee  that 
roots  placed  at  a  specified  location  can  be  made  dominant,  and  a  trial  and 
error  procedure  may  have  to  be  employed  to  achieve  the  best  values  for  the 
various  parameters.  Trial  and  error  may  also  have  to  be  used  in  the  de¬ 
sign  of  cascade  compensators  where  a  specified  root  location  may  require 
parameter  values  that  are  not  physically  realizable. 

In  these  instances,  the  calculations  have  to  be  redone  in  terms  of 
slightly  modified  specifications  or  a  different  means  of  compensation  may 
have  to  be  used.  In  general,  system  specifications  are  not  rigidly  fixed, 
but  can  be  met  by  a  given  range  of  values  or  by  some  upper  or  lower  limit. 

To  avoid  this  trial  and  error  analytical  procedure  one  can  employ  the 
graphical  solution.  If  a  net  of  curves  is  plotted  by  a  computer  or  other¬ 
wise,  one  can,  by  picking  an  M  point  or  operating  point  in  the  parameter 
plane,  read  from  the  curves  the  n  roots  of  the  n^  order  characteristic 
equation.  The  trial  and  error  procedure  can  then  be  done  visually  to  pick 
an  operating  point  which  best  meets  the  given  specifications. 

3-4-2  Some  examples  of  the  graphical  solution. 

In  this  section  only  first  quadrant  values  of  alpha  and  beta  are 
assumed  to  be  of  interest.  The  graphs  were  plotted  with  the  aid  of  the 
computer  program  presented  in  section  (6-1).  At  the  end  of  each  curve  on 
the  graph  is  a  letter  and  a  number.  The  letters  are  abbreviations  for  the 
following  quantities: 

Z(zeta),  S(sigma  or  the  real  part  of  the  complex  variable  S) ,  ZW 
(zeta-omega  or  the  real  part  of  the  complex  roots),  w(omega  or  the  undamped 
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natural  frequency  of  the  complex  roots). 


Example  3-18 

Problem: 

For  the  system  shown  in  figure  (3-22)  set  K  at  the  stability  limit. 
Place  a  dominant  root  pair  within  the  following  region  of  the  S-plane: 

.4  —  ~f  -  .7,  and  2  <  w  5  6.  Both  tachometer  and  acceleration  feed¬ 
back  can  be  used.  If  possible  use  only  one  or  the  other. 

Solution: 

From  figure  (3-22)  the  uncompensated  system's  loop  transfer  function 
is: 

GH  *  -1  =  K/[S(S  +  10)  (S2  +  5S  +  100)],  which  when  expanded  be¬ 
comes: 

S4  +  15S3  +  150S2  +  1000S  +  K  -  0 

To  determine  the  value  of  K  at  the  stability  limit  the  Routh  array  is 


formed: 

1  150  K 

15  1000  0 

1250  15K  0 

1.25xl06  -  225K  0  0 

15K  0  0 


From  the  Routh  array  the  stability  limit  is  seen  to  be: 

K  =  5555.5 

If  both  tachometer  and  acceleration  feedback  are  used  the  compensated 
system's  characteristic  equation  becomes: 

S4  +  15S3  +  (150  +  5555.5  cA  )S2  +  (1000  +  5555.5  jS  )S  +  5555.5  =  0 
where  alpha  =  Kfl,  beta  =  Kt  and  K  has  been  set  at  the  stability  limit. 
Parameter  plane  curves  for  this  characteristic  equation  are  shown  in  figure 
(3-23). 
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Figure  3-22 


From  these  curves,  the  following  analysis  can  be  made.  The  origin 
of  the  parameter  plane  corresponds  to  the  roots  of  the  uncompensated 
system.  Since  the  zeta  =  0  curve  passes  through  the  origin,  two  roots 
are  located  on  the  jw  axis  of  the  S-plane  as  was  to  be  expected  from  the 
Routh  array.  The  other  two  roots  are  also  complex  and  are  located  at 
zeta  =  .8,  and  omega  =  5. 

It  is  important  to  note  that  when  an  operating  point  involves  two 
different  pairs  of  complex  roots,  then  the  curves  for  two  different  values 
of  omega  and  two  different  values  of  zeta  will  have  to  pass  through  the 
point.  To  determine  which  value  of  omega  goes  with  which  value  of  zeta, 
it  is  necessary  to  refer  to  the  computer  printout  data.  Due  to  lack  of 
space,  this  data  is  not  provided  herewith. 

With  K  set  to  zero,  the  effect  of  tachometer  feedback  alone  corres- 
a 

ponds  to  movement  of  the  M  point  along  the  Y-axis.  In  figure  (3-23)  the 
unstable  region  is  determined  by  an  inspection  of  the  way  the  constant 
zeta  curves  tend  as  zeta  increases.  Since  the  Y-axis  is  always  in  the 
unstable  region,  it  is  concluded  that  tachometer  feedback  alone  cannot 
stabilize  the  system. 

The  effect  of  acceleration  feedback  alone  can  be  observed  by  moving 
along  the  X-axis  of  figure  (3-23).  If  K  is  varied  between  .01  and  .06, 
the  system  will  have  two  pairs  of  complex  roots  with  the  following  ranges 
of  values  for  zeta: 

.3  <  zeta  <  .5  and  .25  <  zeta  <  .32 

If  both  tachometer  and  acceleration  feedback  are  used,  it  is  concluded 
that  tachometer  feedback  will  in  general  cause  the  zeta  of  one  pair  of  roots 
to  increase  while  the  zeta  of  the  other  pair  decreases.  It  appears  that 
acceleration  feedback  alone  would  be  better. 

From  the  graph  it  is  determined  that  with  K  =0  and  K  =  .012,  that 

t  3 

complex  roots  are  located  at  zeta  =  .45,  omega  =  4,  and  zeta  =  .32,  omega 
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=  13.  Since  it  is  true  that  (.45)  (4)  =  1.8  <<  (.32)  (13)  =  4.15,  it 
is  apparent  that  the  roots  at  zeta  =  .45  and  omega  =  4  are  dominant. 

The  specifications  have  therefore  been  met  and  the  problem  is  solved. 

Example  3-19 

Problem: 

With  the  product  set  at  the  stability  limit  of  the  uncompensated 

system,  use  position  feedback  as  shown  in  figure  (3-24)  to  obtain  domin¬ 
ant  characteristic  roots  with  maximum  possible  zeta  and  omega,  and  with  a 
dominancy  factor  of  about  two  to  one  (i.e.,  the  ratio  of  the  real  part 
of  any  secondary  root  to  the  real  part  of  the  primary  roots  is  about  two 
or  greater) . 

Solution: 

From  figure  (3-24)  the  characteristic  equation  becomes: 

S3  +  2100S2  +  (  oC+  1.2x106)S  +  100  c/  +  (2  +  108  =  0 
where  alpha  =  and  beta  =  Setting  =  0  and  using  the  Routh 

criteria,  the  value  of  at  the  stability  limit  of  the  uncompensated 

9 

system  is  found  to  be:  2.42x10  . 

In  figure  (3-25)  are  plotted  parameter  plane  curves  for  the  above  sys¬ 
tem.  Constant  zeta-omega  curves  have  also  been  included  to  assist  in  the 
dominancy  considerations. 

The  analysis  proceeds  as  follows.  The  unstable  region  is  to  the 
left  of  the  zeta  =  0  curve.  Root  values  for  the  uncompensated  system 
are  obtained  for  M  point  values  along  the  Y-axis,  since  this  corresponds 
to  K3  =  0.  The  stability  limit  of  the  uncompensated  system  corresponds 
to  the  intersection  of  the  zeta  =  0  curve  and  the  Y-axis.  The  value  of 
beta  =  at  this  point  is  observed  by  the  Routh  check.  The  M  points  of 

interest  therefore  lie  along  a  horizontal  line  drawn  through  this  point  as 
shown  in  figure  (3-25). 

The  maximum  zeta  obtainable  for  any  value  of  alpha  along  this  line  is 
about  .33.  As  alpha  increases  then  zeta-omega  increases  and  the  magnitude 


of  the  real  root  decreases  and  the  root  becomes  more  dominant.  Omega 

also  increases  with  alpha.  Therefore  the  overriding  criteria  is  the 

dominancy  factor.  On  this  basis  an  M  point  located  at  alpha  = 

6  9 

=  2.5x10  and  beta  =  =  2.42  x  10  is  chosen.  Characteristic  roots 

are  read  from  the  graph  and  are  located  at:  zeta  =  .32,  omega  =  1680, 
zeta-omega  =  550,  and  a  real  root  at  -1000.  The  dominancy  factor  is 
computed  to  be  550/1000  or  1.82.  Since  this  is  about  two  and  the  maxi¬ 
mum  zeta  and  omega  have  been  obtained,  the  solution  is  complete. 

It  can  be  remarked  that  since  was  set  at  a  specific  value  for 

analysis  purposes,  the  problem  was  then  reduced  to  only  one  variable  and 
root  locus  techniques  could  have  been  employed.  This  example  illustrates 
the  flexibility  of  the  parameter  plane  and  it  also  points  out  the  interest¬ 
ing  fact  that  root  values  determined  along  either  a  horizontal  or  a  verti¬ 
cal  line  in  the  parameter  plane  constitute  a  root  locus  in  terms  of  the 
variable  parameter.  Root  values  along  a  sloped  straight  line  in  the  para¬ 
meter  plane  constitute  a  sort  of  hybrid  root  locus  where  the  two  root  locus 
parameters  are  linear  ly  related. 

Example  3-20 

Problem: 

1. 

2. 

3. 

4. 

Solution: 

It  is  known  that  steady  state  error  =  9^/Kg,  where  in  this  problem 
K  =  K/10,  and  9  =  1.0.  Error  **  2°  -  .0349  rad.  Therefore,  K  =  28.6 
and  K  =  286.  If  the  above  system  can  be  made  second  order  dominant,  the 
second  order  curves  can  be  used.  For  an  overshoot  of  307»,  a  zeta  of  .35 


Set  K  such  that  the  system  of  figure  (3-26)  will  follow  a 
ramp  input,  9  =  l.Ot. 

Steady  state  error  should  be  less  than  two  degrees. 

Step  response  overshoot  should  be  less  than  307o. 

Find  values  of  and  P  to  meet  the  above  specifications. 
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is  required. 

The  characteristic  equation  of  the  compensated  system  then  becomes: 

S4  +  (11  +  oC  )s3  +  (10  +  11  oL  )S2  +  (10  oL  +  286  @  )S  +  286 oC.  =  o 
In  the  above  equation,  alpha  =  p  and  beta  =  V  .  Fourth  order  parameter 
plane  curves  are  plotted  in  figure  (3-27).  It  is  seen  that  the  zeta  = 

.35  curve  has  a  minimum  value  of  beta  =  25,5.  Since  beta  is  greater  than 
ten,  a  multiple  lead  compensator  or  perhaps  a  lag-lead  compensator  is 
indicated.  To  enhance  the  chances  of  having  to  use  no  more  than  two 
sections,  the  minimum  value  of  beta  25.5  is  chosen  for  the  single  section. 
From  figure  (3-27)  the  following  root  locations  can  be  read  off  correspond¬ 
ing  to  beta  =25.5  and  alpha  =  87:  zeta  =  .35,  omega  =8.2,  and  real  roots 
at  -90  and  -4.5.  Zeta-omega  for  the  complex  roots  is  three.  Since  4.5  = 
1.5x3,  the  compensated  system  is  second  order  dominant  with  a  dominance 
factor  of  1.5. 

At  this  point  either  a  multiple  section  compensator  can  be  designed 
from  the  above  values  of  alpha  and  beta,  or  another  scheme  of  compensation 
can  be  used.  Both  alternatives  are  considered  in  section  (4). 
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Miscellaneous  aspects  of  the  parameter  plane. 


4-1  Some  general  comments. 

In  section  (3-2-2)  it  was  shown  that  constant  zeta  parameter  plane 
curves  of  order  two  through  five  originate  at  a  point  where  alpha  =  M/K 
and  beta  *  N/K,  where  M,  N,  and  K  are  determined  by  the  zero  and  first 
power  coefficients  only.  Inductive  reasoning  can  be  used  to  conclude  that 
constant  zeta  curves  of  any  order  originate  at  this  common  point  which  is 
determined  only  by  the  zero  and  first  power  coefficients.  An  exception 
is  when  K  =  0.  In  this  case  the  origin  point  depends  on  higher  order  co¬ 
efficients  and  its  location  will  be  obvious  given  a  specific  equation.  If 
K  is  not  zero,  the  origin  point  is  independent  of  the  order  of  the  charact¬ 
eristic  equation. 

Inspection  of  the  expressions  for  alpha  and  beta  in  section  (3-2-2) 
indicates  that  the  shape  of  the  constant  zeta  curves  as  omega  becomes 
larger  is  primarily  determined  by  the  coefficients  of  higher  power,  and 
in  general  the  curves  become  more  complex  and  less  well  behaved  as  the 
order  of  the  characteristic  equation  increases.  For  a  given  character¬ 
istic  equation,  an  increase  in  complexity  can  be  observed  as  alpha  and 
beta  appear  in  more  coefficients.  As  indicated  in  section  (3-2-2)  all 
constant  zeta  curves  tend  to  plus  or  minus  infinity.  The  relative  magni¬ 
tudes  of  the  coefficients  determine  whether  the  limit  is  plus  or  minus 
infinity.  It  is  therefore  necessary  to  choose  a  frequency range  of  inter¬ 
est  before  plotting  the  curves,  thus  limiting  the  analysis  to  one  "window" 
of  the  infinite  plane. 

Mitrovic  curves,  since  they  involve  equations  which  have  only  one 
parameter  appearing  in  each  of  two  coefficients,  are  in  general  simpler 
and  more  well  behaved  than  the  parameter  plane  curves.  It  is  also  interest¬ 
ing  to  note  that  when  parameter  plane  curves  are  plotted  for  characteris¬ 
tic  equations  of  the  Mitrovic  type,  the  resulting  curves  are  identical  to 
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those  plotted  from  the  Mitrovic  aquations.  This  can  be  seen  as  follows. 

From  reference  (8)  page  349,  Bq  and  are  given  by  the  following 
relations  where  notation  has  been  changed  to  conform  to  this  text. 


*. ' dkwVi  Bi =  vk'l0k 

where  0k  =  -(2  J  0k-1  +  0k_2)  for  k  ^  2 

and  0  =  0,  0.  =  1 

O  I 

Comparison  of  the  tabulated  0  functions  in  table  (10-1)  of  reference  (8) 
with  appendix  (IB)  of  this  text  leads  to  the  following  relation: 


K  - 


The  Bq  -B^  Mitrovic  characteristic  equation  is  of  the  form: 

Sn+...+d_S2  +  B.S  +  B  =0 
2  1  o 

If  Bq  =  alpha  and  B^  =  beta,  alpha  and  beta  can  be  computed  from  the 
parameter  plane  equations  (2-10)  and  (2-11)  and  are  as  follows: 


alpha  =  2?  (-l)kd,wkU, 

lc?T  k 


k-1 


beta  =  ('1)kdkwk'luk 


If  the  relationship  between  the  0k  and  Uk  functions  is  employed  in  the 
Mitrovic  expressions  for  Bq-B^,  one  obtains: 


B0%^vk<-»Vi  B1  - 12  Vk_1<-1>\ 

It  is  seen  that  Bq  =  alpha  and  B^  =  beta.  The  above  procedure  can  be 
repeated  for  variables  appearing  in  any  two  of  the  coefficients  of  the 
characteristic  equation.  This  duality  property  is  employed  in  section  (4-2) 
to  compensate  parameter  plane  type  characteristic  equations  employing  normal¬ 
ized  Mitrovic  Bq-B^  and  third  order  curves  which  have  been  plotted 

using  the  parameter  plane  computer  program  presented  in  section  (6-1). 

Since  in  the  parameter  plane,  a  negative  value  for  alpha  or  beta 
does  not  necessarily  imply  a  negative  coefficient  of  the  characteristic 
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equation,  one  is  not  restricted  to  first  quadrant  values  of  alpha  and 
beta.  In  most  applications  however  a  negative  value  for  alpha  or  beta 

corresponds  to  a  right  half  plane  pole  or  zero  in  a  cascade  compensator 
or  to  positive  feedback.  For  all  of  these  cases  only  the  first  quadrant 
of  the  parameter  plane  is  of  interest. 

Since  no  stability  criteria,  either  relative  or  absolute,  has  been 
established  for  the  parameter  plane,  it  is  necessary  to  base  the  stability 
analysis  on  observing  which  way  the  curves  tend  as  omega  and  zeta  are 
varied.  For  this  reason  it  is  desirable  to  plot  curves  for  as  many  values 
of  zeta,  omega,  sigma,  and  if  desired,  zeta-omega,  as  is  necessary  to  fix 
the  pattern. 

4-2  Normalized  third  order  curves. 

4-2-1.  Discussion  of  the  normalized  curves. 

Third  order  curves  are  available  for  finding  the  roots  of  polynomials 

where  the  coefficient  of  zero  power  is  one  variable  and  the  coefficient  of 

the  first  power  is  the  other  variable.  Such  a  set  of  curves  is  presented 

in  reference  (8),  and  the  curves  are  called  normalized  MJtrovic  B  -B,  curves. 

o  1 

In  this  section  a  method  is  presented  whereby  parameter  plane  type 
characteristic  equations  of  third  order  can  be  compensated  on  two  different 
types  of  normalized  curves.  One  type  is  the  Bo-B^  curves  of  reference  (8), 
and  the  other  type  is  the  BjB2  curves  which  are  presented  in  this  section. 

The  method  here  presented  is  divided  into  three  cases. 

4-2-2  Derivation  of  the  normalized  transformations. 

Case  I 

Characteristic  equations  of  the  type 

S3  +  d2S2  +  (b^  oL  +  c^  |3  +  d1)S  +  bo<54  +  CQ  |3  +  dQ  =  0  (4-1) 

are  considered,  where  alpha  and  beta  are  the  variable  parameters. 

In  equation  (4-1)  letting  S  *  d2s  one  obtains: 
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+  d0s^  +  (b  +c  (2  +  d  )d  s  +  b  q/L  +  c  fi-f-d  =0 
2  2  1  I'  12  o  o  I  o 

3 

The  £>bove  equation  is  divided  by  d^  resulting  in: 

S3  +  S3  +  sO^oC  +  c  £  +  d  )/d*  +  (bo0i.  +  co  (2+  dQ)/d3  "  0  (4'2> 


Let : 

Bj  -  (b1o«.  +  c1 12  +  d^/d2  (4-3) 

B  =(bo4+c  /3+d  )/d3 

o  o  o'  o  2 

Rearranging  equations  (4-3) ,  two  linear  equations  in  alpha  and  beta  are 
obtained: 


VL  +  Cifs 

b  cf~  +  c 


"  d2Bl  '  dl 


(4-3a) 

(2  ’  d2Bo  -  do 

Equations  (4-3a)  can  be  solved  by  Cramer's  rule  or  by  the  inversion  of 
the  matrix  of  coefficients  to  obtain: 


alpha  «  UdjBj  -  d[)co  -  <dBB0  -  d^l/fl^c,  -  b^J 

bata  =  t  (d|Bo  -  do)bl  -  (d^Bj  -  d^b,]/^  -  b^) 
Substituting  equations  (4-3)  into  equations  (4-2)  one  obtains: 


(4-4) 


s3  +  s2  +  B.s  +  B  =0  (4-5) 

1  o 

Equation  (4-5)  is  a  normalized  Bq-B^  type  equation. 

In  figure  (4-1)  parameter  plane  curves  have  been  plotted  for  equation 
(4-5)  where  B^  is  the  alpha  or  X-axis  variable  and  Bq  is  the  beta  or  Y- 
axis  variable. 

Systems  whose  characteristic  equation  is  of  the  case  I  type  can  there¬ 
fore  be  compensated  on  the  normalized  curves  of  figure  (4-1),  and  the 
necessary  values  of  alpha  and  beta  can  be  found  from  equations  (4-4). 

Case  II. 

Characteristic  equations  are  considered  of  the  type: 

S3  +  (b2  ot  +  c2  £  +  d2)S2  +  (b1<3d  +C]L  |3  +  dx)S  +  dQ  =  0  (4-6) 
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3  i - 1 

7, acting  S  =  \J  d^  s  in  equation  (4-6)  and  dividing  through  by  d^ 


one  obtains: 


sJ  +  (b2  S-  +  c2  |9  4-  d2)s2/d^  +  (bi°^  +  ci  p  *  d^)s/d*  +  1  =  0  (/+“* S * 7) 


7. 

3 


Let: 


(4-8) 


B2  =  (b2c^  +  c2p+  d2)/do 

Bi  -  <bi</  +  'lP  +  V/doVi 

Equations  (4-8)  when  rearranged,  reduce  to  two  linear  equations  in  alpha 
and  beta  as  follows: 

b2oL  +  c2  p  -  A2  -  d2 

»!-'  +  C1  £  *  doBl  -  dl 


(4-9) 


Solving  for  alpha  and  beta  as  in  case  I  results  in: 

z 

I  B,  -  d,  ) c. 

o  1  12  2  1  12 


_L  Z 

alpha  =  [(d^B2  -  d0)Cl  -  (d*Bn  -  d^cJ/O^c,  -  btc0) 


2  1 


(4-10) 


2  ± 

beta  =  [(d^  -  d1>b2  "  (dcB2  *  d2*bl^b2Cl  -  b^) 

As  a  result  of  equations  (4-8),  equation  (4-6)  becomes: 

3  2 

s  +  B2s  +  BjS  +  1  *  0 


(4-11) 


Equation  (4-11)  is  a  normalized  B^-B2  type  equation.  In  figure  (4-2), 
parameter  plane  curves  have  been  plotted  for  equation  (4-11)  where  B2 
is  the  alpha  or  X-axis  variable  and  B^  is  the  beta  or  Y-axis  variable. 

Systems  whose  characteristic  equation  is  of  the  case  II  type  can 
therefore  be  compensated  on  the  normalized  curves  of  figure  (4-2),  and 
the  necessary  values  of  alpha  and  beta  can  be  found  from  equations  (4-10). 
Case  III. 

Characteristic  equations  are  considered  of  the  type: 

S3  +  (b2c4  +  c2  p  +  d2)S^  +  (b^o£-  +  c^  +  d^)S  +  b^  d-  +  c^  |3>  4*  d  =  0 

(4-12) 

Letting  S  =  (b2c4  +  c^  (3  +  d2)s  and  dividing  through  by  (b^c^  +  c2  4-  d2) 
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equation  (4-12)  becomes: 


(4-13) 


Let: 


(4-14) 


In  equations  (4-14)  after  rearranging,  expanding,  and  collecting  terms 


one  obtains: 


and 


(4-15) 


When  the  substitution  of  equation  (4-14)  is  made  in  equation  (4-12),  a 
normalized  B  -B,  type  equation  results.  Systems  whose  characteristic 
equation  is  of  the  case  III  type  can  be  compensated  on  the  curves  given 
in  figure  (4-1),  to  obtain  a  value  of  Bj^  and  Bq  which  can  be  substituted 
into  equations  (4-15).  Equations  (4-15)  contains  the  solutions  for  alpha 
and  beta.  In  the  most  general  case  where  none  of  the  coefficients  in 
equation  (4-15)  are  zero,  a  graphical  solution  of  the  second  and  third 
order  polynomials  can  be  made.  That  is  the  polynomials  of  equation  (4-15) 
can  be  plotted  on  an  alpha-beta  plane  and  the  necessary  values  of  alpha 
and  beta  can  be  read  from  the  curve  intersections.  When  some  of  the  co¬ 
efficients  of  equations  (4-15)  are  zero,  an  analytical  solution  may  be 
easiest,  or  perhaps  a  combination  of  the  two  can  be  used. 

4-2-3  Application  of  the  method. 

Example  4-1  (Case  I  type) 


Problem: 


In  figure  (4-3)  find  the  values  of  K  and  Kc  to  obtain: 

1.  Characteristic  roots  at  zeta  =  .5,  and  omega  =  10. 

2.  Error  coefficient  should  be  greater  than  or  equal  to  6. 

Solution: 

The  appropriate  characteristic  equation  is: 

S3  +  11S2  +  (30  +c/  )S  +(6  =  0,  where  alpha  =  KKt,  beta  =  K. 

Specification  2.  can  be  satisfied  if: 

K  =  K/ (30  +  KK  )  =  /S  / (30  +cJ-  )  ±  .5 

e  t  1 

The  characteristic  equation  is  of  the  case  I  type,  and  the  appropriate 
frequency  transformation  is  therefore:  S  =  d^s  =  11s 

Letting  represent  the  undamped  natural  frequency  on  the  normalized  plane, 

one  can  obtain  via  the  frequency  transformation:  WN  =  *91 

The  value  of  zeta  is  unaffected  by  the  transformation.  From  the 

B  -B-  curves  of  figure  (4-1),  a  value  of  =  .91  and  zeta  =  .5  is  seen 
o  1  N 

to  correspond  to  Bq  =  .075  and  B^  =  .91.  Employing  equations  (4-4)  along 
with  the  appropriate  coefficient  values  one  can  solve  for  alpha  and  beta. 

The  result  is: 

alpha  =  121( . 91)  -  30  =  80  beta  =  (11)3(.075)  =  100 

Since  alpha  =  KK  =  100K  _  QA  . .  .  .  ^  ,,  _  ,  „  , .  _ 

r  t  =  80  it  is  seen  that:  K  =  .8  and  K  =  100  = 

u  t 

beta. 

The  error  coefficient  is  therefore: 

Kg  =  / (30  +  dL  )  =  .91,  which  is  greater  than  .5  so  the  specifi¬ 

cations  are  satisfied. 

Note:  If  the  error  coefficient  was  required  to  be  greater  than  .91,  the 
specified  value  of  zeta  and  omega  would  have  to  be  modified  so  as  to  in¬ 
crease  beta,  decrease  alpha,  or  both. 

Example  4-2  (Case  II  type) 
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Problem: 


In  the  system  shown  in  figure  (4-4)  find  K,  K  ,  and  K  to  obtain 

c  a 

the  following: 

1.  Characteristic  roots  at  zeta  =  .7,  omega  =  10. 

2.  The  error  coefficient  should  be  greater  than  or  equal  to  6. 

Solution: 

The  characteristic  equation  is: 

S3  +  (3  +o^  )S2  +  (2  +  (S  )S  +  K  =  0,  where  alpha  =  KKa  and 
beta  =  KKt*  The  equation  is  of  the  case  II  type  and  the  appropriate 

j. 

frequency  transformation  is:  S  -  d^s,  where  d^  is  as  yet  unknown. 
Therefore: 

l 

w„  =  10d"3  (4-16) 

No 

Employing  equations  (4-10)  with  the  appropriate  coefficients  one  obtains: 

(4-17) 

(4-18) 

(4-19) 

(4-20) 

In  figure  (4-2)  one  can  pick  of  various  values  of  w^  along  the  line 
zeta  =  .7.  The  desired  value  of  w^  is  one  that  satisfies  equation  (4-20). 
Referring  to  figure  (4-2),  the  following  values  can  be  obtained  when  zeta 
=  .7: 

w^  =  .8  and  =  2.2  and  using  equation  (4-20)  it  is  found 
that  =  2.1. 

Also  from  figure  (4-2)  when  zeta  =  .7  it  can  be  seen  that: 


£  I 

alpha  =  d  B  -  3  beta  =  d5B.  -2 

o  2  o  1 

The  error  coefficient  restriction  can  be  satisfied  if: 

d  =  12  +  6jS  or  (d  -  12)/6  =  beta  =  (^B,  -  2 
o  1  o  o  1 

From  equation  (4-16)  it  is  seen  that: 

d^  =  10/w  ,  d^  =  100/w2  ,  and  d  =  1000/w3 
o  N*  o  N  o  N 

Employing  equations  (4-19)  and  (4-18)  one  obtains: 

1000/6w3  -  2  =  1000B , /w2  -  2  or  B  =  1.67/wM 

N  IN  1  N 

The  Bj-B2  curves  of  figure  (4-2)  can  now  be  employed. 
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Figure 


WN  =  .76  and  B^  =  2.2  and  from  equation  (4-20)  :  B^  =  2.2. 

The  latter  set  of  values  are  the  desired  ones  and  the  corresponding  value 

of  B£  is  2.5.  From  equation  (4-19): 

d  =  1000/(. 76)3  =  2280  =  K 
o 

From  equations  (4-17) : 

alpha  =  29.9  beta  =  379 

Therefore:  beta  =  379  -  KKfc  =  2280Kt  or  =  .166  and  alpha  =  29.9  = 

KK  =  2280K  or  K  =  .0131. 
a  a  a 

The  error  coefficient  specification  should  automatically  be  satis¬ 
fied  due  to  the  way  the  analysis  was  performed,  but  as  a  check: 

K  =  K/(2  +  KK)  =  K/(2  +  B  )  ^6. 
e  t  1 

Now: 

K/(2  +  )  =  2280/381  -  6  which  confirms  the  result. 

Example  4-3  (Case  III  type) 

Problem: 

In  figure  (4-5)  find  values  of  gamma,  P,  and  K  to  obtain: 

1.  Characteristic  roots  at  zeta  =  .7,  omega  =  15. 

2.  The  error  coefficient  should  be  equal  to  20. 

Solution: 

The  characteristic  equation  becomes: 

S3  +  (o£  +  5)S2  +  (5<^  +  100  (3  )S  +  100 c4  =  0,  where  beta  = 
gamma,  alpha  =  P,  and  K  =  100  to  satisfy  specification  2.  This  is  recog¬ 
nized  as  a  case  III  type  and  the  appropriate  frequency  transformation  is: 

S  =  (o 4  +  5)s.  The  appropriate  curves  are  the  Bq-B^  curves  of  figure  (4-1). 

Employing  equations  (4-15)  after  dropping  the  terms  with  zero  co¬ 
efficients  one  obtains: 

B^  2  +  (10B1  -  l)o/  -  100  jQ  +  25B1  =  0  (4-21a) 

BqoC3  +  15Bqo4  2  +  (75Bq  -  100)o/  +  125Bq  =  0 
The  latter  equation  when  divided  by  Bq  becomes: 
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(4-21b) 


3  2 

o4  +  15  04  ^  +  (75  -  100/B  )oC  +  125  =  0 

o 

From  the  frequency  transformation  it  is  seen  that: 

wN  =  15/(oC  +  5)  (4-22) 

The  procedure  is  to  find  values  of  Bq  and  w^,  employing  the  zeta  =  .7 
curve  of  figure  (4-1)  such  that  equations  (4-21b)  and  (4-22)  are  simultan¬ 
eously  satisfied.  This  can  best  be  done  by  trial  and  error. 

Trial  1. 

From  the  curves,  Bq  =  .065  and  wN  =  .58.  From  equation  (4-22)  one 
finds  that  alpha  =  20.9.  These  values  are  used  in  equation  (4-21b)  to 
see  if  the  left  side  of  the  equation  becomes  zero.  Hence: 

(20. 9)3  +  15(20. 9)2  -  1465(20.9)  +  125  =  -14895*  0 

Trial  2. 

Try  Bq  =  .07,  then  from  the  curves,  wN  =  .55.  From  equation  (4-22) 
one  finds  that  alpha  =  26.  Using  equation  (4-21b) : 

(26)3  +  15(26)2  -  1355(26)  +  125  =  -17775  *  0 

Trial  3. 

Try  Bq  =  .075,  then  from  the  curves:  wN  =  .5.  From  equation  (4-22) 
one  finds:  alpha  =  25.  From  equation  (4-21b) : 

(25)3  +  15(25)2  -  1255(25)  +  5(25)  =  -250,  which  can  be  consider¬ 
ed  to  be  close  enough  to  zero  for  the  size  numb ers involved.  Therefore 
from  the  curves,  B^  =  .45. 

Equation  (4-21a)  can  now  be  employed  to  find  beta: 
beta  =  [  (B ,  oC.  2  +  (10B1  -  l)o4  +  256^/100 
or 

beta  =  [ ( .45) (25) 2  +  (10x.45  -  1)25  +  25(.45)]/100  =  3.8 
The  final  results  are: 

gamma  =3.8  P  =  25  K  =  100 


106 


4-3  Normalized  parameter  plane  curves  of  higher  order. 

In  the  preceding  section,  essentially  two  types  of  transformations 
were  made  to  normalize  the  third  order  characteristic  equations.  The 
first  being  the  magnitude  scaling,  was  tacitly  assumed  since  the  co¬ 
efficient  of  the  characteristic  equation  was  taken  to  be  one.  Dividing 
through  by  b^  constituted  the  magnitude  scaling.  The  second  transformation, 
which  was  frequency  scaling,was  used  to  make  one  of  the  remaining  coeffi¬ 
cients  unity.  This  left  only  two  variable  parameters,  and  normalized 
curves  could  readily  be  plotted.  If  the  characteristic  equation  is  higher 
than  third  order,  more  than  two  parameters  are  involved  since  only  two 
coefficients  can  be  made  unity,  so  general  normalized  curves  of  higher 
order  than  three  are  not  feasible  except  in  the  special  form  as  follows. 

In  reference  (5),  Choe  introduced  normalized  families  of  fourth  order 
curves  where  the  family  parameter  was  taken  as  the  coefficient.  Para¬ 
meter  plane  transformations  similar  to  those  given  in  section  (4-2)  could 
be  derived  for  the  fourth  order  case,  but  they  would  be  too  complex  to 
be  of  practical  use.  Obviously  normalized  curves  for  higher  than  fourth 
order  are  impractical. 
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4-4  Three  dimensional  parameter  plane  space. 

4-4-1  Discussion. 

Many  compensation  problems  involve  finding  values  for  more  than  just 
two  parameters.  Such  is  true  in  multiple  section  cascade  compensation  and 
combination  feedback  and  cascade  compensation.  This  type  of  problem  can 
be  solved  by  conventional  parameter  plane  methods  if  all  but  two  para¬ 
meters  are  set  at  some  arbitrary  value.  A  more  illuminating  approach, 
however,  involves  the  use  of  three  dimensional  parameter  plane  space. 

(When  the  problem  involves  or  can  be  reduced  to  three  parameters). 

Since  the  parameter  plane  equations  are  obtained  by  equating  the 
real  and  imaginary  parts  of  the  characteristic  equation  to  zero,  unique 
solutions  exist  for  only  two  parameters.  However,  a  third  parameter  can 
be  intorduced  by  plotting  families  of  alpha-beta  curves  in  two  dimensional 
space.  Theoretically  a  three  dimensional  parameter  space  surface  could 
be  plotted;but  interpreting  the  results  would  be  difficult. 

4-4-2  Example  problem. 

In  section  (6-1)  a  computer  program  is  presented  which  will  plot 
families  of  parameter  plane  curves  in  terms  of  a  third  parameter  as  a 
variable.  The  third  parameter  may  appear  linearly  or  non-linearly  in  any 
of  the  coefficients.  The  following  example  shows  one  application  of  the 
method. 

Example  4-4. 

Experience  has  shown  that  if  the  parameter  values  necessary  to  com¬ 
pensate  a  system  using  a  single  section  cascade  compensator  are  not  reali¬ 
zable,  then  the  application  of  tachometer  feedback  will  often  permit  the 
use  of  values  within  the  acceptable  limit.  The  system  shown  in  figure 
(4-6)  is  the  problem  of  example  (3-20)  but  with  tachometer  feedback  employed. 
Problem: 

Using  the  same  criteria  as  in  example  (3-20)  .  use  tachometer  feedback 
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to  find  more  reasonable  values  for  the  cascade  compensator  parameters. 
Solution: 

From  figure  (4-6)  it  can  be  seen  that: 

Kg  =  KP/(10P  +  KKtP)  =  K/ (10  +  KKfc) 

From  example  (3-20),  Kg  =  28.6. 

Hence:  K  =  286  +  28.6KKt. 

The  characteristic  equation  is  found  to  be: 

S4  +  (11  +</-  )S3  +  (10  +  11c £  +  KKfc  |3)S2  +  (10c/  +  KKfcc/  +  K(B)S  +  K^  =  0 

where  alpha  =  P  and  beta  =  gamma.  The  third  parameter  is  seen  to  be  KKt« 

The  above  characteristic  equation  is  seen  to  be  identical  to  that  of  ex¬ 
ample  (3-20)  when  KKt  is  set  to  zero.  The  error  coefficient  restriction 
can  be  incorporated  into  the  above  equation  by  letting  K  =  286  +  28.6KKt. 
This  results  in: 

S4  +  (11  +oi  )S3  +  (10  +  11  ct  +  +  KKtp)S2  +  (10o/  +  KKto£  +  286|S+  28.6KKt^ 
)S  +  286o£  +  28.6KKtcA  =  0 

Third  parameter  values  of  10,  50,  100,  200,  300,  400,  500,  and  1000 
are  investigated,  employing  the  computer  program  of  section  (6-1).  The 
resulting  constant  zeta  and  constant  omega  curves  are  shown  in  figure 
(4-7).  For  simplicity  the  only  constant  zeta  curves  that  are  plotted 
are  for  the  required  value  of  zeta  of  .35.  If  the  M-point  shown  in 
figure  (4-7)  is  chosen,  the  following  values  are  applicable: 

KK  =  300  gamma  =10  P  =  113 

t 

This  results  in  the  compensator  zero  being  located  at  S  =  -11.3,  which  is 
also  a  closed  loop  system  zero.  For  the  above  parameter  values  the  closed 
loop  system  poles  were  found  using  the  digital  computer  and  are  as  follows: 

-11.056  -90.835  -11.054  +  J29. 56 

The  residue  of  the  real  root  at  S  =  -11.056  is  approximately  zero  due  to 
the  closed  loop  zero  at  -11.3,  so  the  complex  roots  which  give  the  zeta 
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of  .35  are  dominant  and  the  problem  is  solved.  It  is  seen  that  tacho¬ 
meter  feedback  increases  the  undamped  natural  frequency  of  the  complex 
roots  by  about  a  factor  of  three  so  the  settling  time  is  decreased  by  a 
factor  of  one-third.  It  is  noted  however,  that  settling  time  is  not  of 
interest  in  this  example. 

At  this  point  one  should  investigate  to  see  if  tachometer  feedback 
alone  could  produce  the  desired  results.  The  Routh  check  is  first  employed. 
The  applicable  characteristic  equation  is: 

S3  +  11S2  +  (10  +  KKfc)S  +  K  -  0 

and 


K  =  28.6  =  K/ (10  +  KK  ),  hence  K  =  286  -  28.6KK  . 
e  t  *  t 

Using  this  value  for  K,  the  characteristic  equation  becomes: 

S3  +  11S2  +  (10  +  KKt)S  +  286  +  28.6KKt  =  0 
The  Routh  array  is: 


1 

11 

-176  -  17 . 6KKfc 
286  +  28.6KKt 


10  KKt 

286  +  28.6KKt 
0 
0 


Since  a  negative  value  of  KKfc  is  required  to  even  stabilize  the  system  , 
it  is  concluded  that  tachometer  feedback  alone  will  not  work. 

4-5  Characteristic  equations  involving  product  terms  of  alpha  and  beta. 
4-5-1  Basic  derivations. 

Assume  that  the  coefficients  of  the  characteristic  equation  are  of  the 


form: 

ak  =  bk°t  +  °k^+  ^k0^  +  \  (4-23) 

Employing  equations  (2-7)  one  can  obtain: 


ei  B1  +  (3  c1  +  cL  (S  A  +  =  0  (4-24a) 

o£  B2  +  £  c2  +  o4(3  A2  +  D2  =  0  (4 -24b) 
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Where : 


™  u  i,  m  U  1, 

ai  =  k?,  'An-i  A2  <?„  <-»v\  <4-25> 

Uk  1*  Uk’  Pl’  Cl*  Dl*  B2*  C2*  and  D2  are  as  defined  in  section  (2). 

Equations  (4-24)  contain  in  general,  two  unique  solutions  for  alpha 
and  beta.  Note  that  if  =  0,  equations  (4-24)  reduce  to  equations 

(2-9)  and  determinants  can  be  used  to  obtain  the  one  unique  solution. 

Solving  equations  (4-24a)  and  (4-24b)  for  alpha  and  equating  the 
results  one  obtains: 


Di  + f3  °i  =  *Di  -  P  °i 
B1  +  f3  A1  Bl  +  PAl 

Equation  (4-26)  can  be  solved  for  beta  resulting  in: 


(4-26) 


beta 


•<AAD  +  ABC>±  V 


*  A  AD  4  A  BC*  ~  4  A  AC  A  BD  (4-27) 


2  A 


AC 


The  deltas  in  equation  (4-27)  are  shorthand  notation  for  terms  of  the  form: 

ABC  ■  Bl°2  •  B2°l-  AAD  =  A1D2  •  A2D1-  etc' 

Proceeding  in  a  similar  manner,  alpha  is  found  to  be: 


alpha  = 


(  A  _  .  +  A  )  + 
DA  wrv  — 


BC' 


^  DA  +  ^  BC^ 


4  ^  BA  ^  DC  (4-28) 


2  A 


BA 


Equations  (4-27)  and  (4-28)  contain  four  solution  pairs  for  alpha  and  beta, 
only  two  of  which  satisfy  equations  (4-24).  To  find  the  two  correct  solu¬ 
tions,  equation  (4-27)  can  be  used  to  find  two  values  of  beta  which  can 
then  be  substituted  into  equation  (4-26)  to  find  the  corresponding  values 
of  alpha.  The  two  solutions  can  also  be  found  from  equations  (4-27)  and 
(4-28)  where  the  two  correct  solution  pairs  are  the  ones  that  make  equa¬ 
tions  (4-24a)  or  (4- 24b)  go  to  zero. 

Constant  zeta  and  constant  omega  curves  can  be  plotted  from  equations 
(4-27)  and  (4-28).  Proceeding  as  in  section  (2),  constant  zeta-omega 
curves  can  also  be  plotted  from  equations  (4-27)  and  (4-28)  where  B^,  C^, 
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D- ,  B  ,  C9,  and  are  defined  by  equations  (2-15)  and  where: 
12—  ^ 


m 


Kr  O 


hkpk-i 


m 

A2  =  £  hk\ 
^=0  k  R 


(4-29) 


Constant  sigma  curves  however  are  no  longer  straight  lines  when  an  alpha- 
beta  product  is  involved.  In  this  case  equation  (2-17)  becomes: 


c4 


yv\ 


i*i 


K-0 


(-l)kbk  k  +  fS  £.o  (-l)kck  k  +  |S  (-1)  \  £  k  <4-30) 


K=0 


rv\ 


+  £  <-Dkdk^k  =  o. 

fc-0  k 


If  one  assumes  the  following  notation: 


DDD  =  £  (-1)  d,  & 

W) 

ccc  =  £ 

(-Dkck^ 

K- 0  k 

K-O 

M  ,  , 

BBB  =  £  (-l)V  ^ 

BC  =  £ 

(-Dkh, 

K-0  k 

K=0 

k 

Then  equation  (4-30)  becomes: 

c-4  BBB  +  ^  CCC  +  @  BC  +  DDD  =  0  (4-31) 

Equation  (4-31)  is  a  special  form  of  a  conic  section  and  can  be  plotted 
by  solving  for  alpha  and  incrementing  beta  over  a  range  of  values  of  interest. 
The  above  equations  are  programmed  for  the  digital  computer  in  section  (6-2). 
Example  problems  involving  the  above  concepts  are  presented  in  the 


following  section. 


4-6  Design  of  double  section  cascade  compensators. 
4-6-1  Discussion. 


Double  section  compensators  can  be  designed  by  use  of  the  parameter 
plane  in  two  ways.  The  double  section  compensator  can  be  made  equivilant 
to  a  single  section  at  a  specific  complex  frequency  of  interest.  Here 
the  system  is  first  designed  using  a  single  section  compensator  but  the 
single  section  parameter  values  turn  out  to  be  physically  unreali: able. 
This  method  has  the  disadvantage  that  control  is  maintained  over  only  one 
pair  of  complex  roots  and  dominance  is  difficult  to  ensure. 

The  second  method  involves  writing  the  characteristic  equation  with 
a  double  section  compensator  inserted  and  drawing  the  parameter  plane 
curves.  Control  over  all  the  roots  is  obtainable  and  the  dominance  pro¬ 
blem  is  much  simplified. 

Both  the  above  methods  involve  characteristic  equations  with  alpha- 
beta  product  terms  appearing  in  the  coefficients,  and  new  parameter  plane 
equations  have  been  derived  to  handle  this  situation  in  section  (4-5-1). 
4-6-2  Design  of  a  double  section  compensator  on  the  basis  of  given  single 
section  parameter  values. 

This  method  was  proposed  by  Hyon  in  reference  (6)  but  Hyon  used  the 
Mitrovic  equations  to  obtain  a  solution.  Parameter  plane  techniques  will 
now  be  employed. 

Let  the  open  loop  transfer  function  of  a  double  section  cascade 
compensated  system  be: 


*1  *2  ^Pl4lKS  +  V(  2>  .  G  =  -1 

(s  +  ppcs  +  p2) 


(4-32) 


The  uncompensated  system's  forward  path  transfer  function  is  G, 
where  unity  feedback  is  assumed.  The  open  loop  transfer  function  of  a 
single  section  compensated  system  is: 


*(s  +  p/  y  )  .  g  =  -i 

(S  +  P) 


(4-33) 
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Assume  complex  roots  are  required  at: 


1  -/ 


2 

1 


(4-34) 


Equations  (4-32)  and  (4-33)  are  equated  and  G  is  divided  out: 


if  2  (s  +  y  ^i><s  +  p2/n2)  = 


(S  +  P1)(S  +  P2) 


(S  +  P/ft  ) 


(4-35) 


(S  +  P) 


After  rearranging,  collecting  terms  of  like  power,  and  dividing  by  S, 
equation  (4-35)  becomes: 


S2(  *  -  ft  x  ft  2>  Si (pi  +  P2) ft  +  P(1  -  ft  1  ft  2)  -  P1  ft  2 


(4-36) 


The  unknowns  in  equation  (4-36)  are  P^,  P2>  ft  ft  2>  an<* 

By  substituting  a  specific  frequency  for  S  and  equating  the  real  and  imagin¬ 
ary  parts  of  equation  (4-36)  to  zero  separately,  two  equations  in  two  un¬ 
knowns  can  be  obtained.  The  parameters  ft  ^  and  ft  2  can  be  pre-set  to  a 
fixed  value,  thus  determining  whether  the  compensator  will  be  double  lead, 
double  lag,  or  lag  lead.  Values  for  the  remaining  unknowns  P^  and  P2  can 
then  be  obtained  from  the  resulting  two  equations. 

It  can  be  noted  however,  that  when  ft  ^  and  ft  are  given  fixed  values, 
the  parameter  plane  techniques  of  section  (4-5)  can  be  applied  directly 
if  one  lets  P^  =  c/  ,  P2  =  jS  ,  and  P^P2  =  •  Equations  (4-27)  and 

(4-28)  give  the  desired  values  of  P^  and  P2  when  the  specified  value  of 
zeta  and  omega  for  the  complex  roots  is  substituted.  Only  one  or  at  most 
two  points  in  the  parameter  plane  are  of  interest  thus  obviating  the  need 
for  plotting  curves.  Hence  straight  analytical  techniques  can  be  used. 
Example  4-5. 

Problem: 

Apply  the  above  techniques  to  design  a  double  section  filter  equivalent 
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to  the  single  section  filter  of  example  (3-20)  at  the  specified  frequency 
of  zeta  =  .35  and  omega  =8.2.  Of  the  possible  solutions,  choose  the  one 
that  makes  the  specified  complex  roots  most  dominant. 

Solution: 

From  example  (3-20)  a  value  of  gamma  =  25.5  was  needed,  indicating  a 
double  section  lead  filter  might  work.  The  value  of  P  was  87.  Let  #  = 

"5  2  =  5,  which  is  a  reasonable  value  for  a  lead  filter.  When  the  above 
values  are  substituted,  equation  (4  36)  becomes: 


S2  +  (21</  +  21  p  -  2118)S  +  25  d-  (3  -  348(  oL  +  p  )  =0  (4-37) 

when  zeta  =  .35:  U  =  -1;  U  =  0;  U  =  1;  U  =  .7 

-1  o  1  2 

From  equation  (4-37)  the  coefficients  can  be  used  along  with  equations 
(2-10)  and  (4-25)  to  obtain  the  following  quantities: 


Al  =  -25 

A2  =  0 

Bx  =  348 

B2  =  -172 

Cx  =  348 

C2  =  -172 

Dx  =  67.2 

D2  =  17397 

then  found  to  be: 

A  A  =  -4.349  x  105 

AD 

A  =  -6.066  x  10^ 

DC 

> 

w 

o 

n 

o 

A  .  =  4.349  x  10^ 
DA 

Aac  -  4300 

Aba  =  -4300 

=  6.066  x  10 
BD 

Using  equations  (4-27)  and  (4-28)  the  solutions  are: 


j  =  16.7  ot  2  *  84,4  f3  i  =  84‘4  £  2  =  16.7 

From  equations  (4-24)  it  is  found  that  aL  ^  paired  with  ^  ^  and  c/.  ^ 
paired  with  2  are  the  consistant  solutions. 

Since  in  equation  (4-37),  the  coefficients  of  alpha  and  beta  are 
identical,  both  solutions  produce  the  same  result,  c/  ^  and  p  ^  are  arbi¬ 
trarily  chosen.  The  characteristic  equation  of  the  compensated  system 
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then  becomes: 


s5  +  112. 14s4  +  2533. 2S3  +  2.3678  x  104S2  +  1.587  x  105S  +  4.0344 
x  105  =  0 

Using  the  digital  computer  the  roots  are  found  to  be: 

-4.205  -16.79  -85.51  -2.82  +  j7 .674 

The  complex  roots  are  at  the  specified  value  of  zeta  and  omega. 

Comparing  the  above  roots  with  those  obtained  in  example  (4-4)  where 
the  single  section  plus  tachometer  feedback  was  employed,  one  concludes  that 
the  double  section  compensator  produces  dominant  complex  roots  whereas  the 
tachometer  feedback  and  single  section  scheme  does  not.  In  example  (4-4) 
however,  the  complex  roots  were  found  to  be  effectively  dominant  since  the 
residue  of  the  nearby  real  root  was  about  zero.  In  this  example,  closed 
loop  zeros  are  located  at  -3.34  and  -16.9.  Therefore,  the  magnitude  of 
the  residues  of  the  complex  poles  and  the  real  pole  at  -4.205  are  almost 
the  same  so  it  is  fortunate  that  the  complex  roots  are  dominant. 

4-6-3  Design  of  a  double  section  compensator  using  general  parameter 
plane  methods. 

This  method  involves  incorporating  the  double  section  cascade  compen¬ 
sator  equations  into  the  uncompensated  system's  characteristic  equation. 

This  technique  has  the  advantage  that  parameter  plane  curves  can  then  be 
drawn  and  control  ii.  maintained  over  all  the  characteristic  roots  rather 
than  only  the  two  specified  complex  roots  as  is  the  case  with  the  method 
of  section  (4-6-2). 

Here  the  technique  can  best  be  explained  by  an  example. 

Example  4-6 

Problem: 

Figure  (4-8)  shows  the  system  of  example  (3-20)  but  with  a  double 
section  compensator  employed.  Solve  example  (3-20)  using  a  double  section 
compensator  as  indicated  in  figure  (4-8) . 
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Solution: 


For  comparison  with  the  results  of  example  (4-5)  let: 


The  characteristic  equation  then  becomes: 

s5  +  (  +  [5»+  11)S4  +  (  dp  +  llo4  +  11(3,  +  10)S3  + 

(llc4p  +  10  oC  +  10  (3  +  7150)S2  +  (10o4  p  +  1430  o4  + 

1430  [2>  )S  +  286  cL  (3:  =0  (4-38) 

Here:  d  =  P^  (3  =  ?2,  and  d  (3  =  P^. 

Parameter  plane  curves  for  equation  (4-38)  have  been  plotted  in 
figures  (4-9)  and  (4-10)  utilizing  the  computer  program  of  section  (6-2). 
Since  the  curves  are  rather  complex,  the  constant  zeta  and  constant  omega 
curves  are  plotted  in  figure  (4-9)  and  the  constant  sigma  curves  are  plot¬ 
ted  in  figure  (4-10).  In  figure  (4-9)  it  is  noted  that  there  are  two 

sections  of  constant  zeta  curves  for  each  value  of  zeta  and  two  section  of 

constant  omega  curves  for  each  value  of  omega.  This  agrees  with  the  re¬ 
sults  proven  earlier  that  two  unique  solutions  exist  for  each  value  of 
zeta  and  omega. 

In  figure  (4-10),  the  discontinuities  of  the  constant  sigma  curves 
are  indicated  by  the  horizontal  straight  line.  The  straight  line  is  not 
part  of  the  curve  and  should  be  ignored.  A  study  of  the  curves  indicates 
that  for  a  given  curve,  if  one  chooses  a  point,  say  (  cd  (2  »  then 

there  exists  a  corresponding  point  (  ^  <d~  ^) •  This  is  due  to  the 

fact  that  the  coefficients  of  alpha  and  beta  are  identical  in  the 
characteristic  equation. 

If  the  M-point  indicated  in  figures  (4-9)  and  (4-10)  is  chosen  then 
the  following  characteristic  roots  are  indicated: 

-9.21  -26.14  -65.34  -2.1  +  j5.61 

The  accuracy  of  the  root  values  was  obtained  from  the  printed  out 
computer  data,  given  that  alpha  =  31.2  and  beta  =  62.5.  The  complex  roots 
are  located  at  zeta  =  .35  and  omega  =  6. 
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The  roots  obtained  in  example  (4-5)  are: 

-4.205  -16.79  -85.51  -2.82  +j7. 674 

Comparison  of  the  two  sets  of  roots  indicates  that  in  example 
(4-5)  the  ratio  of  the  nearest  real  root  to  the  real  part  of  the  complex 
roots  is  1.44  whereas  in  this  example  it  is  4.4.  Hence  the  dominancy 
factor  has  improved  by  a  factor  of  three  in  this  example. 

The  effectiveness  of  the  method  of  this  section  is  thus  apparent. 


5  Root  locus  digital  computer  programs. 

The  programs  of  this  section  were  written  using  Fortran  60  along  with 
subroutines  available  at  the  computer  facility  of  the  U.  S.  Naval  Post¬ 
graduate  School.  It  is  necessary  to  have  the  coefficients  of  the  poly¬ 
nomial  or  characteristic  equation  to  use  the  root  locus  program  of  section 
(5-2),  so  a  program  is  presented  in  section  (5-1)  to  compute  the  coeffi¬ 
cients  in  case  the  equation  is  in  factored  form. 

The  use  of  the  programs  is  explained  in  the  comment  cards  at  the 
beginning  of  the  programs. 


124 


x  a 


LlI 


O  CD 
2  X 
—  O 
h-  a 

h- 

x  x 
x  < 
*  a 
a  a 
•-H  O 
4-  or 

rH  Cl 


a 


00 

a 

CL 

a 

O 

UJ 

a 

•— 1 

O 

X 

a 

00 

X 

1— 

1— 

00 

0 

a 

a 

a 

h- 

0 

a 

X 

00 

0 

0 

< 

O 

!— 

0 

f— 

X 

a 

a 

•c* 

a 

X 

a 

a 

O 

a 

1— 

a 

X 

CD 

>— 

O 

X 

1— 

a 

X 

>— 

_i 

>- 

O 

a 

a 

< 

< 

CQ 

CD 

a 

—  a 

X 

2: 

1-  O 

X 

Ll! 

X 

h- 

< 

0 

oO 

X 

X 

a 

X  00 

X 

<=L 

I— 

< 

co 

a  x 

>- 

X 

0  a 

a 

a 

O 

X 

CD 

1— 

a  — 

0 

h- 

a 

a 

< 

00 

a 

a 

< 

1— 

X 

a 

00 

a 

a 

a 

•  a 

< 

a 

X 

a 

0 

0  x 

1 

a 

X 

00 

a 

r-H  }— 

Ll 

i— 

< 

X 

X 

a 

O 

a 

a 

a 

«— * 

< 

CO  }— 

X 

1— 

0 

-  O 

00 

s~ 

< 

i— 

a 

< 

0 

X 

L- 

1— 

00 

CD 

X 

rH 

00 

— * 

X 

« 

< 

< 

a 

a 

a  x 

0 

UJ 

a 

X 

t 

X 

c 

CO 

6  x 

0 

*— 1 

a 

a 

a 

i —  < 

rH 

a 

Ll 

a 

X 

a 

»*H 

X 

O 

X 

0 

CM 

00 

<  o0 

X 

• 

1— 

X 

1— 1 

a 

a  a 

a 

a 

00 

X 

a 

00 

O 

0 

UJ 

a 

l— 

X 

f— 

a  5“— 

— ■» 

O 

0 

0 

a 

a 

X 

X 

a 

x  a 

O 

a 

*— 

X 

00 

a 

a 

0 

< 

1-  < 

O 

a 

*— I 

3: 

< 

t— 1 

a 

a 

rH 

UJ  < 

> 

a 

1— 

a 

1— 

a 

X 

a 

0 

i— 

a 

< 

a 

O  a 

a 

h- 

a 

a 

CD 

VC 

X 

X 

X 

a 

O 

a 

a 

X 

0 

i— 

00  !— 

UJ 

X 

0 

a 

X 

< 

a 

1— 

r-K 

L- 

•— 1 

a 

t 

a 

a 

a 

a  a 

O 

X 

X 

b 

< 

CD 

a 

0 

<  0 

O 

0 

X 

a 

a 

a 

X 

CL 

rH 

->• 

0 

X 

O 

1— 

00 

00 

0 

a 

i— 

a 

VC 

■~c 

1— 

>  X 

♦— • 

a 

00 

< 

/ _ 

a 

a 

a  a 

a 

LU 

X 

X 

< 

• 

O 

< 

<  —• 

X 

a 

*— 1 

0 

a 

O0 

CL 

x  a 

X 

a 

X 

X 

a 

*— 1 

0 

»— 1 

0 

< 

i— 

a 

a 

a 

a  a 

0 

a 

}— 

►— * 

< 

< 

X 

< 

<  X 

rH 

< 

X 

a 

a 

a 

X  h- 

— • 

X 

a 

X 

X 

i— 

0 

a 

K-H 

a 

< 

X 

< 

a 

< 

• 

1— 

a 

a 

}— 

a 

a 

i— 

00 

a 

a 

a  < 

0 

X 

0 

o0 

<  X 

X 

nr 

x 

• 

0 

X 

i— 

X 

X 

a 

>— 

h- 

1—  >— 

X 

b 

a 

0 

»— « 

t— 1 

< 

a 

*— 1 

a 

a 

i— 

X 

a 

a 

rH 

CM 

co  a 

a 

a 

a 

*— 4 

X 

• 

O 

a 

X 

X 

uo 

X 

< 

VC 

< 

X 

X 

X  X 

a 

*—4 

X 

a 

a 

i— 

j— 

a 

a 

a  00 

a 

X 

X 

>- 

a 

a 

< 

< 

< 

<  X 

a 

i— 

00 

00 

co 

X 

00  X 

a 

a 

a  a 

< 

X 

\ 

\ 


a 


<  ~ 


< 

X 


S 


</) 

QC 

O 

I— 

a 

< 


o  X 
x  cl 


ii  it 


CM 


O  o  —  ~ 
X  rH  CM  CM 
UJ  *—■ 
co  co  m 
O  —  — 
(Ml —  I— 


o: 
—  a 

(M  — 


UJ 

*  X 

—  h- 
^  X 

—  r- 

*— •  rH 

a  * 


\ 

#k 

00 

a 

o 

f— 

VJ 

< 

a 

a 

o 


a 

a 

0 

X 

1— 

X 

0 

a 

a 

»— 1 

O 

< 

1— 

—  a 

a 

< 

a 

0 

O 

H 

X 

i— 

x  >- 

X 

#v 

0 

a 

~  a 

#> 

m 

a 

< 

«— 1  < 

rH 

rH 

Cl 

II  X 

II 

a 

»— 4  »— « 

>—4 

< 

O 

a 

^  0 

#V 

«— 1 

< 

—  < 

•v 

a 

a 

~  X 

CO 

a 

a 

— -  ►— 1 

-w. 

rH 

X 

X 

a  x 

»— 4 

a 

0  »— 1 

a  a 

a 

a 

CM 

0 

X  CM 

w  CM 

X 

rH 

X 

* 

#v 

»—■»—«*>  *  X 


X  < 


x 

< 

UJ 

a 


<  < 
X  X 
X  (X 

o  o 

Ll  Ll 


CO  \Q 
—  O 
h-  CM 
< 

X 

(V 

C 
Ll 


CL 


vO 

O 

CM 

X 

< 

Ll' 

QC 


Is-  rH 


CO 

m  o 

CM  CM 


CM  vO 

m  o 

CM  CM 


< 

JL 

DC 

o 

Ll 


i-  < 
X 
— '  a 
a  o 

Cl.  Ll 


I—  l—  <1—1— 


XXX 


a  a 

CL  CL 


QC 

O 

LL 


a  a 

CL  CL 


\ 


CO  ^0 

m  o 

CM  CM 


X  X 

CL  CL 
CL  CL 


II  •  •  1— l 

*-  o  o  ~ 

II  II 

o-  —  o 


+ 
o 

X  rH  !-< 

II  — 
U  O  CL 

X  x  a 


~  o 
a  o 


—1  cm  CO 
II  II  II 
->  ^ 


X 

CO  a  O  rH 

CM 

O 

rH 

CM 

CO 

r- 

a 

0 

0  On  0  cm  CM 

CM 

m 

in 

m 

in 

i—i 

rH 

X 

CM  CM 

CM 

CM 

CM 

CM 

uuu uuu uuuuuuu 


125 


hr 

,, 

“> 

or 

HH 

ez 

l-H 

or 

CL 

CL 

<— * 

— - 

«— 

or 

or 

— « - 

*— 

* 

* 

_J 

_j 

51 

s: 

*-H 

* 

NwT 

-j 

«— ■ 

— «• 

'W> 

•— > 

■w 

—r 

ora: 

H. 

— 

w- 

W 

hh 

•— i 

HH 

or 

l-H 

or 

HH 

or 

* 

HH 

HH 

HH 

or 

HH 

■ — • 

or 

CL 

CL 

or 

or 

or 

M  l-H 

- - - 

• - 

or 

or 

or 

•— i 

•— i 

* 

* 

* 

* 

* 

or  or 

HH 

HH 

* 

» 

5}C 

cl 

or 

hh 

M 

l-H 

l-H 

l-H 

l-H 

ao 

or 

or 

HH 

HH 

HH 

+ 

I 

or 

or 

or 

or 

or 

or 

i  + 

< 

+ 

or 

or 

or 

*— 

— - 

<  < 

CD 

CO 

u 

u 

*— ■  — . 

— 

— * 

<  < 

CO 

“•> 

l 

+ 

1 

+ 

i 

+ 

z:;z 

“•> 

\ 

+ 

I 

W 

w 

*— 

— » 

' < — - 

o 

T— < 

C\J 

CO 

m 

SO 

Is- 

00 

o 

^n— • 

- 

l-H 

CL 

_J 

-j 

51 

or  ►-* 

CO 

co 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

r-H 

or 

HH 

-J 

or 

CL 

w 

w 

w 

■ — ■ 

w 

w 

or  or 

o 

o  o 

o 

o 

CL 

or 

* 

* 

CL 

HH 

or 

l-H 

or 

l-H 

*  * 

z 

z  z: 

z 

zr 

ON  ON 

On 

ON 

ON 

ON 

ON 

ON 

ON 

ON 

* 

* 

or 

HH 

or 

— * 

.*— » 

CL 

or 

or 

or 

or 

or 

or  or 

or 

or 

or 

hh 

* 

* 

❖ 

* 

* 

or  or 

r-H 

r— I  r-H 

r— I 

r— < 

ON 

ON 

ON 

On 

ON 

ON 

ON 

ON 

ON 

O' 

HH 

HH 

* 

5}C 

V— 

or  or 

or 

or 

CL 

or 

00 

11 

II  II 

11 

11 

*— • 

— • 

— * 

— 

•w 

or 

or 

or 

or 

or 

CL 

CL 

or 

or 

CL 

or 

u  ii 

l-H 

N^ 

-J 

5: 

-J 

2: 

_i 

5! 

_j 

5:  or 

or 

or 

or 

or 

or 

or 

< 

< 

CO 

00 

u 

vj 

1 

i 

1 

i 

l 

l 

1 

1 

1 

i 

or 

or  <  < 

CO 

ii 

ii 

11 

11 

n 

ii 

11 

u 

r-  r- 

O 

ON  ON 

ON 

ON 

l-H 

l-H 

l-H 

HH 

“•> 

-0 

n*c 

_ j 

ii 

ii 

ii 

it 

11 

m 

O 

►— » 

CL 

or 

l-H 

or 

l-H 

or 

•H 

w  w 

w 

■ - ' 

w 

W 

■ — ’ 

• — ■ 

or 

HH 

or 

HH 

or 

ii 

ii 

II 

GT 

CL 

or  cl 

CL 

or 

or 

or 

or-H 

o 

o  o 

o 

O  Li. 

Li. 

Li. 

u_ 

u. 

Li- 

Li. 

u_ 

U- 

u_ 

or 

or  or 

or 

or 

M 

N 

A 

< 

CD 

CO 

o 

o 

uu  ooo 

Cl 

c 

I-H 

HH 

HH 

HH 

HH 

HH 

HH 

HH 

< 

<  CO  CO 

u 

in 

o 

H 

CM 

CO 

<*- 

m 

NO 

Is- 

00  o 

rH 

co 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO  r-H 

126 


CRI=8RR*RI (L)+BRI*RR(L) 
CR(6)=CRR*RR(M)-CRI*RI (M)+CR(6) 
Cl (6)=CRR*RI (M)+CRI*RR (M)+CI (6) 
PRINT  20,I,J.K*L*M 
CONTINUE 


I 


cl 

UJ 

o 

cl 

o 

o 


o 

Z 

UJ 

u 

UJ 

o 


a: 

UJ 

a 

a: 

o 

o 


o 

2: 

UJ 

u 

1/) 

UJ 

o 


cn  cn 

o'  1-4 

yj  u 


1/) 

►- 

z 

UJ 

4—4 

u 


Ul 

UJ 

o 

u 

IL 

o 


yf) 


UJ 

M 

u 

M 

u. 

LL 

UJ 

O 

u 


in 

in 

Sf 

st 

+ 

yf) 

W 

W 

w 

W 

yf) 

\— 

a: 

1—4 

CL 

*— I 

-*) 

b- 

CL 

r-* 

«— s 

u 

u 

yj 

U 

w 

W 

CL 

< 

• 

"> 

") 

+ 

+ 

"> 

+ 

+ 

»-4 

O' 

< 

a 

w 

w 

w 

— * 

O' 

O' 

a 

4—4 

O' 

-i 

-J 

M 

a: 

* 

♦ 

v 

— * 

a 

g: 

— 

— * 

W- 

w- 

O' 

O' 

W 

W 

— 

.j 

—  CL 

U 

* 

♦ 

M 

or 

* 

* 

»-4 

O' 

1-4 

1—4 

< 

U  < 

z 

— > 

W 

■ — 

O' 

a 

— 

— * 

O' 

O' 

w- 

W 

UJ 

z  z 

r 

1-4 

4—4 

4-4 

CL 

* 

* 

4-4 

M 

* 

»— 4 

1-4 

CL 

•*  ►— 4 

rH 

w- 

—  a:  or 

1—4 

4—4 

W 

W 

1—4 

1-4 

CkT 

CL 

X 

-s  O 

II 

1-4 

1— 1 

* 

* 

O' 

O' 

►—4 

1—4 

O' 

O' 

1 

+ 

0 

II  < 

"> 

o: 

or 

►— « 

4-4 

CD 

CD 

CL 

a: 

< 

< 

— * 

— 

— 

<t 

“">  X 

1 

+ 

a:  o' 

1 

+  U 

\ 

+ 

1 

+ 

") 

fVi 

9k 

•»  •— 1 

<  < 

9k 

a 

r*. 

W 

w 

w 

w 

V. 

-  X 

"> 

n 

1 

+ 

-J 

-J  ^ 

“> 

CL 

4-4 

O' 

1—4 

\ 

— • 

">  H 

m  <r 

0 

rH 

w 

w 

w 

W  9k 

CL 

a: 

u  yj 

V 

in 

w  in 

M-4 

CO 

o' 

M 

iZ  CL 

~  “"J 

in 

in 

0 

CL 

►—4 

O' 

* 

* 

+ 

*f 

•V 

• 

O'  •• 

VJ 

*  *  * 

» 

a: 

a: 

w 

w 

O' 

CL  ^ 

O  O 

O 

9k 

•* 

9k 

CL 

O' 

a: 

o:  • 

O 

O 

9k 

r- 

*** 

0  - 

r-» 

0k 

O 

U  X 

W. 

n  r-  n 

r- 

* 

♦ 

a 

M 

* 

Ti/L  1-4 

z  z 

z 

in 

in 

in 

* 

* 

* 

3^  4-4 

z 

z 

fO 

»-4 

1—4 

Z  “) 

-> 

rH 

C\J 

W  S 

*  *  * 

* 

— «. 

— 

o' 

O' 

O' 

CL  * 

— 

— 

CL 

a:  * 

<► 

w 

W 

•>>  w> 

w 

X 

sf 

UJ 

in  x 

go 

m 

r-  e- 

r- 

M 

M 

* 

* 

O' 

QC  *-*  UJ 

rH 

rH 

H 

in 

in 

m 

1-4 

1-4 

CL 

Of  CM  UJ 

rH 

rH 

fO 

CL 

CL 

UJ  *-H  O' 

4—4 

• 

• 

rH 

in 

in 

in  x 

in 

in  0 

^  ^ 

— 

w* 

w* 

CL 

O' 

CD 

co  cm  -5 

II 

II 

II 

«— •» 

w 

w 

< 

<  CM  ID 

II 

II 

CL 

CL 

ro 

II  O'  O' 

rH  O 

W. 

CM 

W 

eg  w 

eg 

eg  O' 

•J  ^  u 

-J 

cl 

o' 

CL 

O' 

II 

II 

z 

1—4 

") 

3^  O' 

CL 

II 

II 

z 

M 

") 

II 

II 

Z  “)  II 

II 

II 

II 

►- 

1- 

i— 

1  l  l 

1 

CL 

O' 

< 

< 

— * 

—  h- 

4—4 

1 

1 

1 

CL 

CL 

r-% 

—  h- 

*— • 

1 

»— 4 

< 

1- 

< 

1-  < 

i— 

1-  O 

*■«  “> 

*: 

II 

II 

II 

II 

in 

m  z 

h~ 

in 

in 

in 

1—4 

n 

1—4 

II 

II 

St 

sf  Z 

i— 

cn 

<n 

1—4 

<n 

cn 

1- 

rH  <\J 

eg 

rH  rH  2l 

Z 

X 

Z  X 

z 

z  1 — 

WWW 

w 

CL 

*— ■* 

O' 

M 

w 

W  M 

z 

W 

w 

w 

CL 

►—4 

w 

w  M 

z 

«— 

W 

W 

z 

w 

W 

W 

W 

o: 

r.  1—4 

CL 

—  g: 

*— 4 

4-4 

U.U.U. 

Ul 

CL 

o' 

CL 

O' 

CL 

-a  O  0 

O 

O 

Ul 

Ll 

U- 

CL 

CL 

CL 

~  C£  0  0 

O 

u. 

O' 

1—4 

0  0  cn 

*— 4 

O' 

1—4 

O 

V  O' 

O 

CL  O 

CL 

CH  O 

*H  *-l  1-4 

M 

< 

< 

CD 

CD 

U 

u  a  u  0 

O 

O 

•—4 

K-4 

*—4 

< 

< 

U 

u  a  u 

O 

O 

I-H 

yj 

u  u  0  yj 

UUUUl 

x  a 

Ul 

a  u. 

9k 

0l 

CL  O 

i 

rH 

rH 

►H  (M  fri 

st 

00 

r- 

cn 

0 

rH 

O 

in 

Cs| 

st 

<n  rH 

rH 

st 

in 

vO 

>t  * 

4- 

rH 

in 

in 

H 

rH 

in 

in 

in 

CM 

eg 

eg 

o  c 


127 


o 

re 

1 — 

UJ 

< 

o 

ce 

c0 

CO 

CL 

< 

X 

UJ 

CO 

< 

3 

i — 

CO 

UJ 

< 

cO 

1 - 4 

<c 

CO 

\— 

Z 

o 

3 

CJ 

< 

_j 

1 — t 

O 

re 

UJ 

GO 

z 

»— 4 

UJ 

ro 

a 

no 

h~ 

M 

UJ 

c0 

h~ 

x> 

< 

1 — 

u_ 

i— i 

i — « 

c0 

• 

. — 

3 

h- 

o 

a 

- — 

co 

a 

3 

O 

1— 

t — • 

• 

LU 

CO 

<1- 

ui 

CJ 

3 

ce 

o 

ce 

1— 

vj- 

1 

h- 

o 

LL 

UJ 

-J 

c0 

UJ 

1— 

1 

•—i  <c 

cj 

3 

(O 

CL 

1— 

V 

o 

f— 1 

»~f 

llJ 

i — H 

^  il 

-J 

c0 

ce 

h~ 

f- 

ce 

0 

X 

< 

a 

U) 

X 

o 

cO 

o 

< 

2 

CL 

z 

X 

X 

Lu 

• — « 

o 

< 

LiJ 

LlJ 

X 

3 

x 

IX 

cO 

Li  J 

ce 

UJ 

co 

no 

x> 

3 

CM 

UJ 

H- 

X 

CD 

7: 

_J 

O 

» — i 

h- 

* 

X 

1 — 

h— 

o 

o 

LO 

— * 

LJ 

X 

•— < 

liJ 

LlJ 

k- 

u 

< 

O 

ro 

x 

in 

o 

re 

X 

CL 

t— 

C0 

X 

i 

o 

< 

<c 

< 

co 

-j 

3 

i: 

H- 

• 

UJ 

o 

CJ 

— 

J— 

cj 

r 

o 

n 

i — t 

Ll. 

O 

< 

H- 

o 

i — • 

Li  J 

re 

-J 

UJ 

3 

< 

re 

< 

» — » 

< 

HJ 

LU 

-J 

CO 

**> 

CO 

3 

X 

h- 

-J 

h- 

UJ 

x 

UJ 

»— * 

h- 

o 

h- 

*— t 

a 

ro 

J— 

CO 

u 

^  - 

CO 

*— < 

1— 

u 

Lu 

i  i 

< 

3 

L  J 

X 

ce 

1 — 

I — 1 

sZ 

i — 

►— f 

u 

CL 

<r 

X 

1 — 

CJ 

oO 

r- 

rj 

o 

c0 

c) 

X 

X 

X 

c0 

M 

10 

o 

LU 

-j 

CO 

uux 

< 

»—» 

;  C 

CJ 

_j 

V 

LlJ 

UJ 

LiJ 

< 

< 

re 

oc 

UJ 

o 

a 

re 

h- 

no 

X 

h- 

tu  ce 

L0 

tu 

3 

3 

tu 

o 

1 — 

1— 

CO 

1 — 

UJ 

\— 

o 

l— 

*— * 

ce 

UJ 

CJ 

i 

i — 

01 

X 

ce 

U) 

cO 

o 

LiJ 

X 

< 

J 

o 

< 

► — < 

3 

< 

cn 

lu 

X 

i — 

ce 

o 

lL 

V 

ro 

1— 

< 

3 

re 

^0 

o 

o 

h- 

cO 

c0 

IU 

X 

* 

UJ 

>7 

X 

• 

CO 

Ll 

O 

u 

< 

< 

.j 

<C 

UJ 

• — i 

co  lu  ce 

o 

o 

:o 

< 

O 

CL 

ce 

< 

UJ 

LU 

CD 

h- 

fL 

•» 

-j 

L0 

a 

UJ 

< 

CJ 

LU 

X 

X 

z 

o 

UJ 

< 

1— 

7:7 

>—* 

J— 

►— 4 

3 

UJ 

re 

> 

UJ 

X 

CO  CO 

< 

i — * 

-J 

L~ 

a 

70 

«c 

X 

h- 

CO 

1— 

-J 

lu 

K- 

L  J 

:o 

CD 

1— 

o  oc 

a 

O 

3 

3 

O 

O 

"7 

UJ 

CO 

< 

o 

J— 

X 

3 

3 

3 

i  ) 

>— i 

>- 

H- 

cn 

u 

c0 

o 

e 

* — 4 

h- 

no 

lu 

ce 

CO 

LU 

U) 

1" 

> 

h- 

< 

ce 

> 

e  < 

o 

•—s 

UJ 

3 

3 

Cf 

o 

< 

n 

re 

h- 

LU 

CO 

ce 

ld  cj 

CO 

X 

i~ 

UJ 

IU 

< 

h- 

o 

z  o 

cO 

< 

OC 

L— 

c0 

re 

> 

uj  o 

iJ 

LU 

IU 

— •  3 

01 

re 

i — « 

2?* 

LU 

ro 

co 

X 

X 

LU 

1 —  LC 

UJ 

co 

• 

: 

UJ 

<i 

UJ 

h~ 

1 — 

X 

H- 

■  7 

0 

o 

X 

CL 

CD 

< 

X 

UJ 

1— 

3 

T  " 

re 

CO 

CO 

h- 

c0 

>  H- 

1 — 

X  <c 

< 

a 

LU 

i — t 

UJ 

1 — 1 

C\l 

co 

*  ce 

tf 

OC 

h- 

• 

LU 

h- 

_J 

_ l 

3  a 

O 

cO 

IU 

J— 

c0 

GO 

1 — 

< 

CL 

Q 

ro 

Q 

H  C) 

O 

►— t 

n 

CO 

:o 

CO 

2 

S’ 

ce 

ce 

re 

-T  X 

X 

X 

re 

-J 

-j 

O 

-J 

*— « 

o 

< 

< 

< 

f~i  Cl 

CL 

h- 

co 

CL 

Q_ 

1 — 

a 

U. 

CJ 

u 

CJ 

u 

O 

OJ 

O 

3 

• 

• 

CJ 

cj 

u 

cj 

CJ 

CJ 

o 

c 

c 

c 

u 

CJ 

CJ 

O 

H- 

ce 

< 

LU 

>r 

f\l 

ce 

co 

LU 

Lu 

'10 

— 

1 — ■ 

I— 

#» 

< 

in 

h- 

< 

:r 

• 

o 

X! 

OsJ 

re 

o 

X 

ce 

1 — 1 

o 

« — i 

o 

lu 

UJ 

J— 

Lu 

i — 

3 

co 

cO 

< 

3 

UO 

w* 

3 

CO 

X 

h- 

• 

*— • 

ce 

*— i 

o 

re 

o 

h- 

1 — ( 

LU 

LU 

UJ 

ro 

•» 

•* 

co 

re 

— 

m 

— 

— 

o 

i— 

— h 

• 

CO 

< 

O 

o 

CM 

o 

• 

co 

s; 

H 

i — i 

r — 1 

CM 

re 

X 

re 

1 

UJ 

i — 1 

LlJ 

• — * 

o 

H 

— ' 

i — « 

<w 

a 

CO 

X 

H- 

rsJ 

oc 

X 

nr 

• 

o 

< 

a 

LU 

m 

b 

LU  Q 

r« 

CJ 

• 

ce 

OC 

re 

— 

CD 

co 

o 

lu 

3 

* 

o 

X 

LU 

r— t 

— ' 

CO 

- — 

CM 

ro 

LU 

—t 

CD 

f-H 

CO 

• 

Lu 

r\j 

X 

X 

— 

a 

r— i 

X 

LU 

►— * 

UJ 

h- 

< 

CJ 

LU 

X 

X  OC 

#» 

<3 

tu 

3 

X 

< 

1 — 

— 

IU 

3 

ro 

< 

u 

o 

CO 

a 

ao 

< 

a 

< — * 

o 

Ovj 

< 

*— * 

co 

u. 

ce 

r-H 

X  ~  X  ~  \ 

*— •  X  <C  3  X  51  + 

<  >  (DO  O  h  h- 


U) 

> 

t-H 

i — i 

o 

X 

t — 

3 

o 

CO 

— ' 

o  — 

<  — 

LU 

X 

h— 

t— 

OC  h- 

CL  X 

3 

X 

1— 

3 

i — i 

o 

\ 

i — « 

H— 

cO 

X 

•» 

—  o 

>* 

CJ 

3 

UJ 

o 

00  OC 

X  h- 

w-» 

Jl 

o 

ro 

o 

CM  * 

<  X 

3 

CO 

< 

O  <~\i 

rH  OC 

X  < 

3 

IU 

CJ 

1 — 

rH 

—  1— 

—  X 

UJ 

c0 

3 

tu 

3  O 

co 

o 

H- 

3 

3 

3 

X 

u-  o 

<  3 

CJ 

X 

<r 

3 

O  X 

X  < 

UJ 

> 

3 

< 

— 

O  ^ 

►— »  3 

1 — 

>—* 

o 

CJ 

o 

X  O 

X  X 

X 

u 

3 

cO 

CM 

*  X 

—  3  X 

< 

1 — » 

OC 

i — 1 

^  •* 

c^ 

in  rH  co 

1— 

3 

»--4 

UJ 

X 

X  X 

CM 

•  * 

c0 

3 

1 - - 

CO 

X 

OC 

CM  < 

r — l 

m  \  x 

X 

UJ 

*-4 

7T 

< 

^  > 

• — 1  X  X 

o 

O 

X 

3 

OC 

X 

w  #V 

r-H 

3  X  X 

CJ 

CJ 

1 - 4 

X 

co 

o 

3  X 

ll 

o 

00  X  X 

1 — • 

#> 

X 

• 

w  ^  w 

3 

in 

vO 

r- 

CO 

cO 

—  X 

o 

X 

1 —  1 —  1 — 

X 

ON  O 

m 

u 

• — i 

-t 

<c  <  <c 

a 

co 

o 

o 

CO  3 

<m  :r 

H 

u 

ll 

x  :l 

re 

ce 

ce 

oc 

oc 

X 

rH  2* 

CO 

ro  x  ce  x 

< 

Cc 

< 

< 

< 

» - 4 

—  d 

o 

O 

< 

o  o  o 

u 

CJ 

CJ 

CJ 

cj  a 

<r  cj  o 

X 

X 

3 

3  3  3 

in  vo  co  o  cr\ 

■ — t  o  vo  r-  vo 

CM 

U  U  CJ  U  CJ 


N- 


Z 

o 


< 

O 

O 

UJ 


X 

u 

X 

HH 

X 

h- 

X 

c0 

+ 

• — * 

CO 

OC 

l 

►—1 

3 

ir 

H- 

< 

CJ 

1 — 

< 

< 

X 

CO 

< 

X 

X 

H- 

X 

- — 

CJ 

3 

X 

CM 

CM 

a 

* 

vD 

f-H 

3 

sO 

f-H 

X 

3 

X 

r» 

M 

3 

f-H 

r^- 

h~ 

rH 

h- 

1— 

II 

u 

ll 

il 

3 

*— « 

*— < 

►— » 

3 

» — t 

• — « 

X 

1— 

•* 

O 

1— 

X 

X 

f— f 

1— • 

X 

► — 1 

• 

r- 

CL 

-w 

3 

w 

• — i 

< 

31  3 

O 

3 

UJ 

OC 

3 

3 

X 

3 

3 

X 

CO 

H- 

i— 

O 

1— 

1— 

X 

X 

f— t 

*— • 

X 

»— • 

M 

X 

rH 

h- 

1 — 

vO 

m 

I— 

1— 

X 

.«-H 

*— 1 

*— • 

cn 

• 

•— « 

»— 1 

*» 

*» 

w> 

<w 

* 

C 

00 

- — 

o 

w 

rH 

X 

•* 

X 

Cc 

< 

CO 

• — i 

•* 

r» 

X 

X 

0 

0 

X 

* 

vO 

1 — * 

3 

O 

o 

rH 

X 

CM 

0 

0 

X 

00 

c< 

— 

a 

o 

CM 

cn 

CM 

CM 

w 

CM 

C 

1 — 

h— 

1 — 

CM 

CM 

H- 

h~ 

1— 

O 

<C 

< 

C 

< 

h-  < 

H- 

1— 

1— 

< 

h- 

u 

o 

CO 

X 

2; 

X 

X 

X 

2: 

X 

c 

X 

X 

oc 

< 

< 

X 

►h  ce 

1— t 

»— • 

t— » 

X 

1— 1 

< 

o  o 

O  3 

UJ  o 

X  0 

X 

X 

X 

0 

X 

u 

3 

3 

3 

X 

X 

3 

a 

Ll  a 

CL 

a 

3 

CL 

a 

O 

cn 

<r 

<r 

00 

O 

o 

o 

CM 

CM 

CM 

CM  CM  CNJ  CO 


128 


I 


N 

S 

♦* 

cz 

LU 

o 

cz 

o 

o 

2: 

a 

UJ 


CZ 

UJ 

o 

cz 

o 

o 
z 
►— « 
a 
2: 

LU 

U 

tO 

LU 

a 


o 

UJ 


LU 

_J 

< 

U 

to 

>- 

* 

UJ 

u 

< 

u 

to 


U 

LU 

X 

tO 

UJ 

-J 

9> 

UJ 

_j 

CO 

< 

UJ 

O 

00 

< 

CL 

-1 

< 

tO 

1— 

— * 

2! 

»— » 

CL 

« 

— * 

►— « 

X 

H— i 

CL 

< 

LU 

— 

m 

1— 

X 

< 

> 

GO 

X 

0 

r— i 

X 

to 

> 

X 

1 

* 

1— 

UJ 

O 

X 

UJ 

CO 

UJ 

2: 

LU 

N*x  X 

1— 

X 

• 

< 

CZ 

1 

UJ 

X 

I'¬ 

* 

rH 

_J 

< 

1 

►— « 

1 — 

tO 

UJ 

* 

•» 

u 

LL 

UJ 

CL 

►— « 

H 

to 

t-H 

Ll 

O 

0 

< 

1“ 

♦ 

0 

Ll 

O 

< 

O 

O 

cz 

Ll 

LU 

u 

tO 

O 

O 

LU 

UJ 

to 

=> 

UJ 

LU 

CL 

— * 

— >» 

X 

x| 

O 

1— 

U 

0 

u 

0 

O 

•> 

U 

2: 

< 

0 

CL 

2: 

M 

X 

UJ 

> 

Ll 

*— 

a 

\— 

1— 

UJ 

I— 

— * 

t—i 

O 

X 

O 

rH 

rH 

O 

1— 

2: 

u 

—  2T  -J 

\ 

X 

X 

O 

II 

II 

O 

to 

< 

2T 

►—» 

2:  ^  < 

CL 

X 

UJ 

cz 

X 

X 

CZ 

>- 

H— 

•* 

.-i 

Ll 

*  iH  >— < 

LU 

♦ 

»— 

♦ 

•» 

* 

*  — «* 

to 

tO 

rH 

11 

Ll 

H  ||  1 — 

CO 

LU 

tO 

0 

— 

CL  h- 

2: 

II 

X 

UJ 

ll  N-Z  •— * 

V 

U 

> 

2: 

X 

X 

h-  tO 

UJ 

O 

X 

•* 

0 

X  *  IX 

b 

< 

LU 

to 

•» 

w 

O  < 

.X 

U 

* 

u 

#>■  » — » 

2: 

u 

UJ  u 

X 

CL 

►— « 

O  _J 

I— 

m 

X 

X 

X 

—  X  X 

X 

to 

-i  < 

UJ 

#> 

1— 

1— 

CZ 

m 

• 

1— 1 

X 

w 

OO 

X  —  Ch 

cz 

rH 

X 

< 

U 

UJ 

X 

a 

O 

0 

•*  rH 

X 

0 

• 

(NJ 

>* 

< 

—  CQ  (NJ 

cz 

< 

CO 

O 

m 

u  to 

u 

I— 

< 

O 

0 

0  ^ 

rH 

2: 

0 

rH 

•* 

< 

•» 

CO  •» 

< 

> 

*s 

O  2: 

* 

tO  X 

< 

rH 

cz 

cz 

2:  CO 

XI 

•> 

rH  LU 

X 

«w 

•> 

N 

^  \ 

> 

* 

X 

2  * 

\ 

X 

•> 

u 

X 

•»  ^ 

1— 

w 

> — 

* 

•» 

CO 

UJ 

CO 

\ 

#> 

r- 

\ 

•*  r-  x 

X 

*  CO 

\ 

* 

<1- 

to  »-H 

rH 

rH  S 

CL 

*> 

is  r- 

X 

0 

00 

w 

N 

(NJ 

m 

0 

N 

com  0  x 

m 

0 

N 

v£) 

CO  0 

N 

r- 

0 

X  (NJ 

O 

II  — 

>- 

ON 

CO 

0 

CO 

<£  * 

X 

(NJ 

(\l 

0 

(NJ 

(NJO  (NJ  — 

(NJ 

O 

(NJ 

(NJ 

0  (NJ 

— 

(NJ 

0 

(NJ 

II 

(NJ 

X  < 

U 

vO 

NO 

r- 

vO 

CL  (NJ 

+  <  <  < 
o  2: 

20:0:0: 
II  o  o  o 
z:  ll-  u_  il_ 

m  r-  r\J 
OO  INJ 
(M  CM 


(NJ 

I—  1“ 

2:  a  e 

►— «  <f  «— t 
O'  LU  O' 
0.0:0. 


CZ 

o 


co 

r\i 


cm  l— 
h-  I —  <t  h- 
ZQ  Z  2:  2 

H<  >— t  /V  >— t 

czm  0:00 

0.0  Q.  Ll  CL 

m 

(nj 


<nj 


a 
<  •- 
LU  CC 
O'  CL 


2:  **: 


CL 

O  a: 

lu  a 

vO 

(NJ 


(NJ  h- 

I-  < 

a  2:  s: 

<  *—  CZ 

uj  a:  o 

o:  q.  ll 

r- 

(NJ 


(\j 


2:  a 

M  <t  ►— « 

O'  LU  CL 
Cl  CL  CL 


1 —  J  <  I- 


5:  2 

O' 

o 

Ll 


(nj 


—  —  CO 

II  O  O  *  I—  iZ 

~  rH  CL  I —  I —  I —  I —  C\J  <\J  <  —  r-H  O 

vo\:ijzzzzj^iis  co  + 

2:  -z  jmmmm  jvooq:h-^ 

iiocliicczczczczc  »  o  a  11  u.  11  o 

s:o<s:uaaaaurH5:u_'^^^o 


D 

(M 


r- 

vO 


nj 

o 

(NJ 


CO  (NJ 


129 


[ 


UJ 

-J 

< 

u 

tO 

>- 

UJ 

_l 

< 

u 

tO 

X 

•V 

UJ 


— 

m 

>— ♦ 

\ 

o 

h— 

Ci 

\ 

> 

l 

»— < 

Z 

1— 

UJ 

‘  * 

•> 

* 

M 

• 

CO 

«—■ . 

in 

? 

H 

< 

o 

• 

* 

-1 

<r 

o 

l-L. 

s 

* 

1 

* — i 

z 

1 — 

CO 

Ch 

UJ 

►— 4 

o 

* 

co 

•t 

o 

o 

* 

II 

f— 

(X 

O 

co 

o 

< 

* 

— * * 

cO 

x 

s: 

a 

a 

* 

* 

UJ 

UJ 

f- 

IX 

x 

Is- 

INJ 

x 

o 

o 

O 

i — 

co 

< 

o 

z 

o 

•v 

2: 

IX 

* 

* 

o 

vO 

UJ 

to 

* 

« — i 

i-H 

X 

co 

1 

h- 

o 

a 

ii 

II 

•v 

+> 

i 

to 

X 

IX 

v; 

>: 

in 

> 

UJ 

o 

* 

i — 

co 

• 

U) 

INI 

z 

* — 

' — 

O  h- 

r«. 

i 

* 

V 

O  to 

<r 

in 

^  UJ 

-1 

X 

X  < 

CO 

* 

—  X 

-1 

#1 

*  —1 

i 

GO  f- 

< 

fsj 

i — 

1 — 

a  * 

CO 

i 

x  x 

X 

< 

o 

o 

X  r-H 

CO 

i 

vO 

-vOO 

m 

z 

z 

» 

o 

o 

o  * 

#k 

i-H  »-H  X 

c\j 

z 

z 

f— 

X 

X 

Z  00 

C\J 

«-W 

+  *  UJ 

#V 

•k 

*■  — 

IX 

■ — ' 

' — ’ 

'  9» 

CO 

i— H 

X  x  INj 

\ 

i—H 

I-H  — *  C\j 

-H  S 

> 

•V 

* 

1 

^  \  ^ 

X 

+ 

II  ^  o 

II  — ' 

J  Choo  o 

00 

<  * 

UJ  H  vO 

i 

o 

CD  '  p-  ao 

o  — 

00  O 

-J  —  c\j 

^  ct 

O  o  vO  r- 

vO 

X  cm 

D  CO  h 

z 

X 

1  1 — 

i — 

a: 

m 

II 

CL 

o  * 

Z  —  O 

ii 

O 

ll  <  l —  o 

h-  < 

o  o 

O  ll  H  h 

NO  — 

i-H  1 — 

1 —  1 — 

►— 

(\J 

(\|M  OO 

QZOlZh2  2:hZH-  +  z  z 
x  — '  x  x  ►-«  m  ^  u  ._j  ^  ~  z 
OU_UJOXOXOOZO--'ll  x  II 
z«No.ae)Q.ii-ozQa^LL5: 

r-i  <t  r-  vO  O  in  O 


^  l  -J 


ON  II 
O  <  S 

vO 


<  x  x  x 
u  0-  a.  cl 


U  vO 
<  * 

O  r-H 


z 

o 

u 

co 


I —  •  H- 


HiA  H  (\JH  4 
rH  f\J  CO  CO 

O  O  O  O  O  O  O  O 
•  i —  •  h-  •  y—  ♦ 


H  CO  H  CO  rH  CO 

oo  <r  vo  <r  oo 

<t  m  so 

O  o  o  o  o  o 


•  I—  •  *—  •  f— 


O  ll  O  ll  O  II  o  II  O  ll  O  ll  O  ll  O  ll 


CO 


(M 

co 


co 

co 


co 


in 

CO 


vO 

CO 


r- 

co 


00 

CO 


‘'f -***•*>»* 


I 


LU 

-I 

CD 

< 

CL 

< 

> 


Ns 


cc 

< 

CL 


r- 

o 

• 

rH 

II 

o 

rj\ 

<0 


fH  O  H 

>t  oo  ro  w 

»— 

O  O  i-  < 
F-  •  Z  X 

i-H  *— i  CL 
O  II  CC  Q 
OOCUl 

O  «— i  o 

<*  CO 


O 

u_ 

O) 

< 


X 

o 

<\J  «-H 

<r  — 


< 

z  2; 

UJ  CL 
CL  O 
CL  .Li- 


CC 

< 


<  II 
IU  -) 
DC 

X  -H 

O'  O 

*  H 

X 

co  o 

r-4  O 


H  LlJ 

11  a 
x  1 


>  rH  + 

*  II  ~ 


X 

m 


o 

vO 


fM 


O  h- 

o  < 

f-H  :> 

CC 

O  O 
o  u. 

o 

vO 


X  Z  rO  ~ 


CL 


I- 

< 

T 

DC 

O 

u_ 


HO  H 
II  01  — 
3  —  3 
3  U-  3 

i0  — ' 
f~“*  iU 

rr  < 


rH  “■)  LU 
r-  3  3 


O  u. 
a  h 


o  o 

O  iTi 
CO 


3  O  3  *—•  F— 
I —  1 —  1 —  h—  Z 

o  z  -« 

OOO  CC 

KD  CL  U  CL 

(\|  r— 4 

vO  r*- 


O 

O 

a 


o 
o  «-» 
* 

CL  ~ 

<t  1—4 

>  I 
II  3 

CC  — 
<  lu 

>  >H 

o 

*5 


o 

CO 

c\J  |  CO 


-r- 
15  *  vO 
<(\J  o 
(L  *  cm 
Qrvj 
*  O 


II  3  II  -JvO  F— 
O  —  O  U  * 
Oil  O  <ch  o 

s:  ~  s:  u  <3 


00 

LU 


LU 

O. 

* 

r-H 

X 

<  j 

UJ 

CO 

1 — 

X 

|  -i 

s_~ 

| 

►— 4 

< 

LU 

LL 

X 

/ 

u 

IU 

0 

O' 

O) 

9 

* 

i  >- 

IS) 

X 

* 

00 

LU 

•> 

CO 

LU 

CM 

3 

1 - 

* 

-J 

H 

U 

CC 

1 — 

< 

< 

< 

(L 

U 

> 

CL 

< 

O) 

h- 

CL 

X 

0 

O 

O 

9s 

0 

LU 

<c 

-J 

LU 

LL 

►—1 

X 

< 

-J 

* 

Ll_ 

*— 1 

LU 

L— 

»— « 

X 

CC 

r--« 

no 

O 

O' 

X 

(M 

LU 

9 

O' 

*— < 

1 — 1 

CL 

X 

* 

CM 

-S 

* 

c/) 

CO 

X 

sO 

0 

no 

cc 

•k 

CO 

#> 

z 

< 

L— 

LU 

1 — 

> — 1 

rH 

•> 

-J 

O 

X 

DC 

» 

SO 

H 

* 

O 

1 — 

< 

1 — 

* 

II 

fM 

CL 

CL 

CL 

r~H 

*— « 

9s 

X 

< 

x> 

#> 

O 

— s 

0 

a 

— ■* 

O 

O'  *— • 

Li_ 

< 

— * 

• — 4 

X 

CM  h- 

IU 

0 

X 

<1- 

9s 

f-H  O 

c/) 

cl 

< 

0 

0 

3 

3 

—  O 

1 — 

X 

s> 

Vh 

1 

1 — 

1 — 

CC  iC 

O 

O' 

M 

#> 

LU 

O 

O 

9  ♦ 

O 

9 

X 

CO 

• 

O 

O 

—  CL 

CC  — 

X 

O' 

• 

m 

rL 

(L  — 

F- 

x  x 

<r 

* 

CM 

1 

— ’ * 

*h- 

CM  Q 

fH  \ 

•s 

X 

CL 

» — 1 

3 

— - 

O' 

Q0 

fLO) 

l 

—  O 

vO  Ns 

CL 

^o 

< 

LU 

3 

*—* 

O' 

O' 

j— < 

Q  CC 

*  Ns 

< 

#» 

> 

<\J 

3 

.  — 

* 

* 

O-J 

X 

9  9 

Ns  » 

> 

1— 

— 

* 

a 

3 

* — i 

CC 

f-H 

00 

O  ^ 

1 — 

-  O 

\  JO 

CL 

CO 

— >• 

»— 1 

1 — 

h~ 

O 

O' 

(LrH 

0 

rH  Z 

\  3 

X 

<  O 

O 

• 

-J 

O  H 

O 

0 

rH 

* 

♦>  9s 

0 

CO  •» 

\  O 

CL  CO 

rH 

CM  CL 

Z  ~ 

CO 

Z  1-  • 

O 

0 

* 

•-H 

Ooo 

CL 

—  CL 

+  _J 

9 

* 

#> 

f-H  <C 

GO 

♦  0  0 

CC 

(L 

f-H 

O 

Z  ^ 

U  < 

LU 


3 

o 

CL 
O  no 
Z  3 
LU  10 


> 

z  * 

O  CC 

•— « 

cO  Z 

z  o 

LU  X 
*2.  X 

'  o 
o  u 


Z  CO 


II 

< 

s: 

CM 

O' 


O' 

O' 


00  »-H 

O'  o 


•  • 

O  O 

II  II 

<  < 

x  x 

LU  LU 
IU  LL. 

CM 

O' 

<r 


«r.  a*# 


131 


20  BETAN  =R00T I ( M )  +EE ( M ) 
ALFAN=R00TR(M)  +FE(M) 

’  DO  7  1=1*100 
S=2.*ALFAN 

J  =  - ( ALFAN**2  +  6E  TAN**2 ) 


1 


oo 

CL 

+ 


< 

Li_ 

— I 

< 

* 

CM 

—  2: 
QL  < 
ijc  I — 

•  UJ 
CM  X 
* 


CM 

—  * 
~  * 


CM 

* 

* 


< 

O 

U) 

►H 

X 


< 

CO 

o 


< 

(X 

UJ 


+ 

V, 

CM 

2 

2 

CM  • 

X 

CM 

1 

O 

Q 

*  (N 

rH 

O 

+ 

+ 

—  * 

O 

#k 

CM 

X 

2 

CM 

CM 

O 

CM 

CM 

* 

* 

* 

* 

• 

rH 

rH 

1 

* 

• 

CM  O 

* 

* 

rH  • 

2 

«k 

-J 

2 

(M 

* 

* 

2 

+  rH 

M 

CM 

w 

< 

~  * 

2 

U 

U 

2  + 

rH 

1 

l 

U 

h- 

+ 

CM  2 

< 

<  2 

m 

X 

1 

1 

* 

X 

— 

1  < 

X 

\ 

\ 

X  < 

* 

u 

h— 

no 

—r 

CM 

2  Y— 

X 

X  1 — 

X 

2 

<r 

1 

< 

«— 

—  X 

< 

2 

2 

<  X 

•k 

X 

1 

+ 

O 

CM 

Y- 

X 

\ 

n  x 

1 

O 

U 

—  0 

vO 

0 

X 

2 

* 

UJ 

*  * 

— 

* 

* 

\ 

— 

CL 

• 

-  2 

— * 

CD 

* 

no 

+ 

h-  — 

rH 

2 

—  \ 

<r 

X 

in 

2  ^ 

rH 

U 

rH 

2 

* 

2 

♦ 

+  CM 

1 

rO 

X 

2  — 

1 

1 

> 

1 

w  *— » 

1 

* 

< 

2 

2  < 

—  1 

O 

+ 

1 

<  2 

X 

X 

2 

CL  h- 

-J 

2! 

* 

X 

< 

<  Y- 

i-i  O 

2 

2 

T3*’ 

X  < 

• 

• 

O 

Y—  O 

X 

— * 

< 

r- 

X 

X 

X  X 

— * 

1  2 

U 

O 

_ 1  1 — 

m 

m 

U 

2 

O  O 

rH 

j 

rH 

u 

X 

.-k 

1 — 1 

< 

.J 

.J  CD 

rH 

2  — 

Q 

❖ 

* 

<  X 

1 

1 

O  CL 

* 

w 

* 

-J  0 

< 

<  * 

— - 

—  O 

2 

2 

1  X 

r- » 

— * 

O 

*— » 

CL  | 

X 

X 

cO 

< 

2 

* 

* 

* 

*  — 

O 

O  * 

* 

< 

< 

<  1 

2 

< 

< 

h- 

rH 

* 

+ 

1 

— - 

t—i 

— * 

--k 

—  rH 

* 

*  • 

2 

X  < 

2 

X 

f-  O 

II  II 

•k 

CO 

u 

—k 

U 

CM 

1 — 1 

H 

rH  w 

<0 

3 

cO  CM  < 

1 

+ 

CL  H* 

X 

X  X  0 

X 

+ 

2 

X 

rH 

★ 

>— 

k— 

—  X 

rH 

+ 

2 

+  1 

Y~ 

2 

2 

X  X 

X 

X 

< 

< 

vO 

< 

X  CL 

—» 

rH 

•k 

•w 

+ 

2 

X 

X 

X  * 

w 

*- 

— »  — » 

X 

< 

< 

<  X 

X 

X 

X 

h- 

4- 

II 

II 

w 

rH  r-f 

rH 

CM 

CO 

X 

O 

< 

CO 

* 

* 

ON 

>!<  • 

Os 

u 

CM 

CO 

2  O 

X 

X 

Y— 

X 

X 

0 

— - 

CM 

X 

+  + 

r- 

O 

w 

rH 

II 

2! 

1— 

1 

• 

• 

rH 

•  CM 

rH 

rH 

II 

—  2 

* 

X 

X 

X 

X 

< 

X 

in 

rH 

2 

2 

no  2  2 

O 

CM 

X 

X 

+ 

-J 

II 

X 

0 

CO 

X 

CM 

II 

U 

1 

2 

<J  — 

• 

< 

X 

II  II 

nO 

nO 

11 

II 

1— 

w 

>— 

X 

—  >— 

2 

II 

II 

O 

u 

X 

2 

O 

II 

O 

II 

O 

II  U 

CM 

II 

II 

X 

X 

2 

2 

h-  < 

O 

CL 

r— 1 

< 

CL  ~ 

rH 

1  rH 

O  CL 

— k 

2 

CM 

—k 

II 

II 

H- 

II 

1 — 

— * 

2 

CO 

< 

< 

< 

< 

< 

< 

2  2 

I— 

I— 

1— 

h-  I— 

+ 

.  2  + 

h- 

3 

rH 

CM 

II 

-1 

II 

II 

rH 

CM 

II 

2  II 

II 

X 

1 — 

X 

t— 

l-H 

CL 

O  O 

O  O 

2 

-  2 

Y- 

— - 

u 

O 

<w 

2: 

2 

X 

2 

2  O 

2  2 

O 

■ — ■ 

D 

0 

—  2 

2 

X 

X 

X  X 

X 

X 

X 

X 

cc  0 

O 

O  O 

X 

O  O 

II 

x  iT 

O 

X 

u  u 

2 

Q 

u 

< 

X 

r— » 

u 

O 

O 

U 

O 

O  O 

2 

0 

O  U 

O 

< 

X 

X  X 

1 

— » 

< 

X 

CL 

X 

O 

CL 

CL 

r— * 

CL  CL 

2 

O 

CL. 

CM 

r- 

1 — i 

CO 

CO 

0 

m 

r- 

O 

vO 

rH 

X 

cm  m 

X 

rH 

CM 

rH 

rH 

m 

« — i 

rH  rH 

rH 

o 

2 


122 


6  Parameter  plane  digital  computer  programs. 

The  programs  of  this  section  were  written  using  Fortran  60  along 
with  subroutines  available  at  the  computer  facility  of  the  U.  S.  Naval 
Postgraduate  School. 

The  use  of  the  programs  is  explained  in  the  comment  cards  at  the 
beginning  of  the  programs. 
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224  FORMAT  (  14X * 6HALPHA+* 15X »5HBETA+ ♦ 14X *6HALPHA- * 15X* 5HBETA-,  1.5X,5HOME 
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7-  The  complementary  roles  of  the  parameter  plane  and  root  locus. 

As  was  mentioned  previously,  parameter  plane  curves  are  infinite  in 
extent.  The  curves  also  exhibit  discontinuities  and  become  less  well 
behaved  as  the  order  increases.  In  addition,  choosing  a  good  graph  scale 
involves  a  good  deal  of  trial  and  error  or  pre-plotting  calculations. 

The  root  locus  can  be  very  useful  in  overcoming  some  of  the  above  diffi¬ 
culties  . 

The  general  procedure  is  to  use  the  root  locus  to  limit  the  range 
of  parameter  values  of  interest,  and  to  gain  additional  insight  into  the 
problem.  This  is  particularly  true  if  a  computer  is  used  to  plot  the  root 
locus.  Since  a  correct  root  locus  graph  scale  can  be  fixed  in  one  or  two 
computer  runs,  several  root  loci  can  be  plotted  in  fairly  short  time. 

The  parameter  plane  can  then  be  used  to  complete  the  problem. 

The  technique  can  best  be  illustrated  by  the  following  example. 
Example  7-1 

Problem: 

For  the  system  shown  in  figure  (7-1)  find  values  for  K^,  K^,  K^, 

and  to  give  a  good  transient  response.  Low  settling  time  is  the 
primary  consideration.  Suggest  modifications  to  improve  the  system. 

The  following  parameters  have  fixed  values: 

P  =  6.28  x  106  P3  =  2.5  x  105  w  =  6.28  x  106 

zeta  =  0.5 
Solution: 

•It  is  known  that  the  block  diagram  elements  of  figure  (7-1)  that 
involve  and  were  added  on  the  basis  of  physical  reasoning  to  help 
reduce  the  settling  time.  The  basic  uncompensated  system  consists  of  the 
single  loop  unity  feedback  control  system  involving  K^,  K and  K^.  From 
the  block  diagram  reduction  of  figure  (7-1)  the  characteristic  equation 
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is  computed  and  the  given  parameter  values  are  substituted.  Letting 

alpha  =  K  K  and  beta  =  K,K-K0  one  then  obtains: 

14  1  Z  O 

S5  +  1.281  x  107S4  +  (8.194  x  1013  +  o^-)S3  +  (2.678  x  102° 

+  p+  2.5  x  105cyL  )s2  +  (6.18  x  1025  +  6.28  x  106jg  + 

K1K2K3K4K5)S  +  3.94  x  1013|S  =  0  (7-1) 

To  study  the  effect  of  the  compensator,  is  set  equal  to  zero 
in  equation  (7-1)  and  parameter  plane  curves  are  plotted  in  figure  (7-2) 
for  the  remaining  system. 

The  unstable  region  for  these  curves  is  above  and  to  the  right  of 
the  zeta  equals  zero  curve.  From  the  Routh  check,  the  stability  limit 

19 

gain  of  the  uncompensated  system,  i.e.,  is  equal  to  1.025  x  10 

This  value  corresponds  to  point  B  in  figure  (7-2).  The  effect  of  the 

compensator  is  to  make  the  system  less  stable  as  is  increased.  With 

19 

^1^2^3  Set  t0  x  10  ,  a  root  locus  is  plotted  in  figure  (7-3)  with 

as  the  variable.  This  plot  shows  the  system  to  be  unstable  for  all 
K1^4  8reater  ^an  zero.  Gain  changes  have  small  effect  on  the  right  half 
plane  root  locations. 

g 

Now  K,K.  is  arbitrarily  set  to  9 . 1  x  10  and  K-K^K0  is  left  unchanged. 
14  1  Z  J 

The  effect  of  the  compensator  is  seen  from  figure  (7-4)  to  have  a 
stabilizing  effect.  It  is  known  that  settling  time  is  inversely  propor¬ 
tional  to  the  undamped  natural  frequency  of  the  complex  roots.  In  figure 
(7-4),  point  A  represents  a  zeta  and  omega  that  would  be  satisfactory  for 
the  complex  roots,  but  the  location  of  the  real  root  shows  the  latter  to 
be  dominant. 

Two  avenues  of  approach  appear  to  be  open.  Try  a  different  means 
of  compensation,  or  modify  the  suggested  means.  The  former  is  investi¬ 
gated  first. 

The  effect  of  tachometer  feedback  around  the  entire  forward  path, 
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From  the  root  locations  as 


with  =  K,.  =  0,  is  shown  in  figure  (7-5). 
tachometer  gain  increases  from  point  A  to  B  to  C,  it  is  seen  that  the 
real  root  is  highly  dominant.  Tachometer  feedback  would  therefore  be  un¬ 
satisfactory.  Since  the  system's  operating  frequency  is  in  the  mega¬ 
cycle  range,  higher  forms  of  derivative  feedback  would  not  be  practical. 

The  fact  that  the  compensator,  which  is  the  feedback  path  around 
in  figure  (7-1),  makes  the  system  less  stable,  suggests  that  this 
path  could  be  eliminated.  When  this  is  done  the  resulting  system  is 
shown  in  figure  (7-6)  and  the  characteristic  equation  is  as  follows: 

S5  +  1.281  x  107S4  +  8.194  x  1013S3  +  (26.67  x  1019+  @  )S2  + 

(6.18  x  1025  +  6.28  x  106(3  +d  )S  3.94  x  1013^  =  0  (7-2) 

In  equation  (7-2),  alpha  =  TUTHTHTli^Yi  and  beta  =  K^Kj. 

The  advantage  of  this  scheme  over  the  original  system  is  that  the 
compensator  can  be  realized  by  an  R-L-C  circuit,  whereas  due  to  the 
physical  nature  of  the  problem,  the  original  scheme  would  require  electro¬ 
mechanical  implementation  with  the  inherent  disadvantages. 

In  figure  (7-6),  is  again  set  at  the  stability  limit  of  the 

uncompensated  system  and  a  root  locus  with  variable  is  plotted  in 

figure  (7-7).  A  study  of  figure  (7-7)  shows  that  a  slightly  better  domin¬ 
ance  factor  can  be  obtained  with  the  modified  system.  It  is  concluded 
that  the  latter  system  not  only  performs  as  good  or  better  than  the  former, 
but  it  is  simpler  and  cheaper  to  implement. 

To  see  if  system  performance  can  be  improved  by  increased  gain,  the 
forward  path  gain  is  increased  by  a  factor  of  ten.  The  resulting 

root  locus  in  figure  (7-8)  indicates  that  the  system  is  unstable  for  all 
K^K^.  This  illustrates  the  difficulty  in  choosing  the  best  values  for 
system  parameters  by  root  locus  techniques  when  more  than  one  parameter 
is  involved. 

The  parameter  plane  can  now  be  advantageously  employed.  The  resulting 
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curves  are  plotted  in  figure  (7-9)  with  variables  anc*  KlK2K3^4K5* 


A  close  study  of  the  curves  indicates  that  the  M-point  shown  is  perhaps 
the  best  one.  The  five  characteristic  roots  can  be  read  directly  from 
the  curves  and  are  as  follows: 


zeta  =  .55 
zeta  =  .71 


omega  =  3  x  10 

£ 

omega  -  6.5  x  10 


S  =  -2.2  X  106 

£ 

The  damping  factor  of  the  first  pair  of  roots  is  1.65  x  10  which 
when  compared  to  2.2  x  10  ,  indicates  that  these  roots  are  dominant. 

In  the  original  uncompensated  third  order  system,  for  a  zeta  of  .55, 
the  maximum  obtainable  omega  for  the  dominant  roots  is  only  about  .3  x 

g 

10  .  On  the  basis  of  the  dominant  roots  only,  one  can  conclude  that  the 
compensation  has  reduced  the  settling  time  by  a  factor  of  ten.  The  trans 
ient  response  of  the  compensated  system  to  a  unit  step  input  was  computed 
by  digital  computer  and  is  given  in  figure  (7-10).  The  settling  time  is 
about  four  microseconds  and  the  overshoot  is  50%. 
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8-  Conclusions. 


Parameter  plane  techniques  have  been  applied  to  the  compensation  of 
linear  feedback  control  systems.  In  particular^ general  equations  have 
been  derived  for  the  cases  of  feedback,  cascade,  and  combination  feedback- 
cascade  compensation,  to  enable  one  to  analytically  place  a  pair  of  com¬ 
plex  roots  at  a  specific  location  in  the  S-plane,  while  simultaneously 
satisfying  the  steady  state  accuracy  requirements. 

A  dominancy  technique  has  been  introduced  whereby  once  the  specified 
pair  of  complex  roots  is  fixed,  the  remaining  roots  of  the  characteristic 
equation  can  be  placed  to  ensure  that  the  specified  roots  are  dominant. 

Sketching  techniques  are  developed  enabling  one  to  quickly  sketch  the 
zeta  equals  zero  and  zeta  equals  one-half  curves.  This  can  be  useful  in 
determining  the  type  of  compensation  to  use  and  in  choosing  an  appropriate 
graph  scale  for  the  digital  computer  plots. 

Graphical  solutions  on  the  parameter  plane  are  discussed  in  terms 
of  engineering  example  problems. 

Miscellaneous  aspects  of  the  parameter  plane  are  discussed.  In  parti¬ 
cular,  transformations  are  derived  to  enable  one  to  compensate  parameter 
plane  type  characteristic  equations  on  normalized  Mitrovic  third  order 
Bq-B^  and  B^-B£  curves.  Three  dimensional  parameter  plane  space  is  dis¬ 
cussed  and  applied  to  an  example  problem. 

A  derivation  of  parameter  plane  equations  involving  product  terms  of 
alpha  and  beta  is  made  and  the  results  are  applied  to  the  design  of  double 
section  cascade  compensators.  Double  section  compensators  are  designed  on 
the  basis  of  unrealizable  single  section  parameters  and  by  incorporating  the 
double  section  compensator  into  the  characteristic  equation  and  plotting  the 
parameter  plane  curves. 

Digital  computer  programs  are  introduced  that  one  can  use  to  plot  root 
loci  and  parameter  plane  curves.  Parameter  plane  curves  can  be  plotted  for 
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characteristic  equations  involving  three  parameters  and  alpha-beta  pro¬ 
duct  terms. 

Finally  an  engineering  example  is  presented  that  points  out  the 
complementary  nature  of  the  root  locus  and  parameter  plane. 

A  basis  for  further  investigation  involves  plotting  constant  band¬ 
width  curves  on  the  parameter  plane  and  determining  the  nature  of  para¬ 
meter  plane  curves  resulting  from  characteristic  equations  containing 
squared  terms  of  alpha  and  beta. 
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APPENDIX  II 


A.  Synthesis  of  R-C  lead  network 


£ 


f?2  £< 


In  figure  (a),  — 

s-fp-v  r j^±el 

°  L  p. 


Fig.  a 

5 j~  */A?t  C  i 


] 


S  +  X 

s  +cfr 


where  (^,£,  and  cLz  Rz/(  /^,-f  ft ^  j  •  The  above  transfer  function 

has  a  D.  C.  gain  of  c£-  so  to  make  the  D.  C.  gain  unity  an  amplifier  of 

± 

gain  will  have  to  be  added.  The  pole  to  zero  ratio,  ^  is  greater  than 
one,  indicating  that  the  circuit  of  figure  (a)  is  a  lead  network. 

If  the  pole  to  zero  ratio  and  'X  become  very  large  the  filter  be¬ 
haves  like  a  pure  differentiator  and  noise  problems  arise.  If  the  pole  to  zero 
ratio  s  kept  less  than  ten,  the  noise  problem  is  reduced. 

B.  Synthesis  of  R-C  lag  network 
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K, 
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Fig.  b 


In  figure  (b)  the  transfer  function  iS  : 


Cz(R,+  Rz) 

S  + 


] 


S  4-  y 

s  +  ± 


Cz(Rt+Rz) 

where  +  Fj>2)  and  o£  =  Ri./(R,+  Rz)t 

The  D.  C.  gain  of  the  above  transfer  function  is  unity  and  the  pole  to 
zero  ratio  is  alpha.  Since  alpha  is  less  than  one,  the  circuit  of  figure 
(b)  is  a  lag  network. 
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